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EXECUTIVE  SUMMARY 


Objective 

The  objective  of  this  study  was  to  develop  guidelines  for  using  finite  element  analyses 
(FEA)  to  predict  the  life  of  new  and  untested  electronic  devices  used  in  military  equip¬ 
ment.  Previously,  there  had  been  an  inability  to  predict  the  life  when  empirically  derived 
failure  rates  were  not  available.  The  techniques  developed  under  this  study  are  directly 
applicable  to  any  advanced  electronics  acquisition.  These  procedures  have  been  used  on 
the  Air  Force  Advanced  Tactical  Fighter  (ATF)  and  the  Navy  Advanced  Tactical  Aircraft 
(ATA)  programs  to  predict  the  life  of  critical  electronic  components. 

Finite  element  analyses  are  computerized  mechanical  engineering  techniques  which 
make  it  possible  to  predict  material  response  when  a  modeleu  device  is  subjected  to  some 
internal  or  external  loading  or  environmental  disturbance.  FEAs  allow  physical  deflec¬ 
tions.  material  stresses  and  material  temperatures  of  complex  devices  to  be  predicted  be¬ 
fore  the  devices  are  fabricated  and  tested.  Although  this  computer  aided  engineering 
(CAE)  technique  can  successfully  predict  mechanical  performance,  a  need  exists  to  ex¬ 
trapolate  FEA  results  to  a  prediction  of  life,  time  to  failure  or  probability  of  failure.  The 
latest  developments  in  FEA  technology  relate  to  the  interfacing  of  the  results  of  FEA 
with  reliability,  or  life  prediction  methodologies.  Even  though  it  is  recognized  that  FEA 
simulations  cannot  address  all  possible  failure  mechanisms,  improper  fabrication  proce¬ 
dures,  etc.,  design  evaluations  and  reliability  assessments  of  electronics  from  a  mechanical 
integrity  and  strength  of  materials  perspective  is  achievable  given  the  proper  geometry, 
material,  boundary  conditions,  loading,  and  strength  information. 

Findings 

This  study  investigated  and  assessed  the  latest  developments  in  FEA  technology.  It  ex¬ 
panded  the  current  mechanical/structural  engineering  FEA  application  techniques  to 
electronic  device  and  equipment  applications,  and  resulted  in  a  practical  guide  for  elec¬ 
tronic  reliability  assessment  purposes. 
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Guidelines  were  developed  to  significantly  improve  the  accuracy  of  finite  element  anal¬ 
yses  of  electronic  devices.  A  method  has  been  outlined  which  allows  simplified  linear 
FEAs  to  be  used  instead  of  the  more  complex  elastic-plastic  nonlinear  FEA.  Guidelines 
for  mesh  generation  have  been  established  which  minimize  arithmetic  errors  caused  when 
materials  with  large  stiffness  differences  are  adjacent  to  each  other.  The  accuracy  of 
FEAs  when  dealing  with  very  small  dimensions  has  been  verified.  Procedures  for  com¬ 
bining  various  loadings  in  order  to  predict  life  have  been  established  for  materials  which 
exhibit  stress  relaxation  and  for  those  which  do  not. 

Existing  computer  codes  were  analyzed  for  applicability  to  electronic  equipment  along 
with  their  ease  of  use,  versatility,  and  computational  accuracy.  The  math  models  current¬ 
ly  used  to  predict  life  were  reviewed  and  documented.  Materials  used  to  manufacture 
electronic  equipment  and  their  respective  failure  mechanisms  (which  can  be  addressed  by 
FEA)  including  buckling,  deformation,  rupture,  fracture,  property  deterioration,  fatigue, 
and  creep  were  identified  and  documented.  Additionally,  the  available  physical  proper¬ 
ties  of  the  materials  which  are  necessary  as  inputs  to  FEA  were  collected  and  docu¬ 
mented.  A  number  of  step-by-step  examples  were  developed  to  illustrate  how  to  predict 
the  life  of  new  untested  electronic  devices  by  utilizing  FEA  outputs. 

Even  though  the  procedures  developed  under  this  study  can  be  directly  applied  to  existing 
electronic  programs,  there  are  some  limiting  factors  to  institutionalizing  the  results  of  this 
study.  These  limiting  factors  are  technology  areas  that  require  further  development: 

1)  Validation  of  the  prediction  techniques:  Life  predictions  using  FEA  outputs  differ 
by  as  much  as  an  order  of  magnitude  from  analytical  techniques  such  as  those  of 
Steinberg  and  Engelmaier.  None  of  the  techniques  have  been  verified  by  labora¬ 
tory  tests  to  confirm  the  accuracy  of  the  life  predictions.  Each  "expert"  comes  up 
with  different  answers  to  the  same  problem.  Program  managers  are  at  a  quandary 
because  each  "expert”  believes  their  approach  is  the  right  one.  However,  none 
have  been  verified.  This  is  probably  the  single  most  limiting  factor. 
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2)  Many  material  properties  necessary  for  input  to  finite  element  codes  are  nonexis¬ 
tent.  In  many  cases,  properties  for  bulk  materials  are  used.  Materials  sized  at  di¬ 
mensions  for  electronic  devices  do  not  behave  as  materials  sized  at  Duik  level 
dimensions.  Questionable  input  data  translates  to  questionable  analysis  results. 

3)  Virtually  every  major  electronic  house  and  aircraft  integrator  are  using  some  type 
of  FEA  procedures  and  reliability  prediction  tools.  However,  these  are  individual 
tools  which  require  tedious,  manual  operations  to  finally  arrive  at  a  life  prediction 
for  electronic  devices.  Electronic  design  engineers  require  expedient,  timely  an¬ 
swers  to  their  question  regarding  their  designs.  Fully  automated  computerized 
procedures  are  needed  to  perform  electronic  life  predictions. 

Study  Approach 

The  study  spanned  18  months  and  was  subdivided  into  four  tasks.  A  brief  description 
of  each  task  follows: 

1)  The  first  task  was  to  determine  the  state-of-the-art  of  methodologies  utilizing  fi¬ 
nite  element  output  results  for  the  prediction  of  reliability.  This  was  accomplished 
through  an  extensive  literature  search,  telephone  interviews,  and  personal  contact 
interviews  at  various  meetings  and  symposia.  The  techniques  utilized  for  predict¬ 
ing  life  of  electronics  by  virtually  every  prominent  electronics  house,  aircraft  inte¬ 
grator,  university,  and  government  agency  were  identified  and  described.  The 
algorithms  and  failure  models  of  each  methodology  were  identified.  The  reliability 
prediction  techniques  which  use  finite  element  analysis  were  reviewed  and  docu¬ 
mented.  The  above  was  accomplished  for  commercially  available  codes  which  ad¬ 
vertise  or  claim  to  be  able  to  predict  life  or  reliability,  and  for  techniques  used  in 
industry  which  are  not  marketed  as  public  domain  commercial  software. 

2)  The  second  task  was  to  select  a  wide  variety  of  materials  which  are  used  in  elec¬ 
tronic  applications  such  as  microelectronic  and  electronic  devices,  modules,  circuit 
boards  and  packages.  Failure  mechanisms  for  the  materials  which  could  be  ad¬ 
dressed  by  finite  element  analysis  were  identified.  The  failure  mechanisms  were 
viewed  from  a  perspective  of  material  failures  rather  than  device  failures  (ie., 
cracking  of  a  material  due  to  stress  rather  than  fatigue  failure  of  a  specific  device 
due  to  thermal  cycling).  Deformation,  buckling,  rupture  or  fracture,  change  or  de¬ 
terioration  of  material  properties,  fatigue  and  creep  were  among  the  failure  mech¬ 
anisms  which  were  considered.  Failure  theories  and  material  behavior  models 
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were  evaluated  and  appropriate  ones  selected  to  predict  life.  The  proper  inter¬ 
faces  of  the  failure  theories  and  material  behavior  models  with  the  appropriate 
FEA  results  were  identified.  The  prediction  techniques  selected  were  practical 
and  useful  for  reliability  assessments. 

3)  The  third  task  was  to  develop  a  description  of  the  above  and  present  in  a  format 
which  would  allow  easy  application  of  the  methodologies  selected  in  Task  2  to 
electronic  equipment.  The  result  is  the  format  of  this  final  report. 

4)  The  fourth  task  was  to  select  examples  and  to  demonstrate  the  application  of  mak¬ 
ing  reliability  assessments  using  FEA  results  coupled  with  material  behavior  mod¬ 
els.  Numerous  examples  were  generated  which  show  step-by-step  procedures  for 
predicting  the  life  of  electronic  devices  using  FEA  outputs.  This  task  was  accom¬ 
plished  by  actually  performing  FEA  using  various  finite  element  codes  and  com¬ 
paring  results.  Known  models  were  utilized  to  validate  accuracy  of  small 
dimension  modeling.  Finally,  examples  were  constructed  which  illustrate  predict¬ 
ing  life  for  materials  exhibiting  the  various  failure  modes  and  mechanisms  pre¬ 
viously  identified. 
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Chapter  1 

STATE-OF-THE-ART  ASSESSMENT 


1.0  Introduction 

This  chapter  describes  the  present  state-of-the-art  of  methodologies  utilizing  finite 
element  analysis  (FEA)  results  for  the  prediction  of  reliability.  The  algorithms  and  fail¬ 
ure  models  utilized  by  these  methodologies  are  identified  and  described.  Reliability  pre¬ 
diction  techniques  that  use  finite  element  analysis  are  reviewed  and  documented. 
Commercially  available  codes  that  advertise  or  claim  to  be  able  to  predict  life  or  reliabil¬ 
ity  are  assessed,  as  are  various  techniques  used  in  industry  which  are  not  marketed  as 
public  domain  commercial  software. 

The  assessment  of  the  state-of-the  art  was  accomplished  by  conducting  an  external 
and  internal  (to  McDonnell  Aircraft  Company,  hereinafter  referred  to  as  MCAIR)  litera¬ 
ture  search  and  information  survey  (conducted  over  the  telephone)  to  identity  current  fi¬ 
nite  element  reliability  analysis  procedures  (Figure  1-1).  Information  sources  included 
universities,  aerospace  companies  active  in  the  design  of  electronics  equipment.  FEA  code 
suppliers,  various  government  agencies  and  McDonnell  Douglas  Corporation  suppliers. 
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The  conclusion  reached  is  that  there  are  various  techniques  to  predict  electronics  reli¬ 
ability  (life)  utilizing  FEA  results.  However,  the  following  are  technology  areas  that  re¬ 
quire  further  development: 

1)  Procedures  have  not  been  accurately  validated  by  laboratory  testing. 

2)  Most  available  packages  use  separate  FEA  procedures  and  life  prediction  tools. 

3)  There  is  no  consensus  as  to  the  accuracy  of,  nor  the  exact  techniques  for,  utilizing 
finite  element  modeling  techniques  for  the  extremely  small  scales  required  to 
analyze  electronics  components. 

4)  Material  data  required  for  analysis  is  ve  '  limited. 

The  following  paragraphs  describe  the  information  obtained  in  researching  the 
siate-of-the-art  of  methodologies. 

1.1  Literature  Search 

In  recent  years,  publications  on  the  application  of  FEA  to  predict  microelectronic  reli¬ 
ability  have  been  appearing  at  an  explosive  rate.  This  report  will  discuss  some  of  the 
more  pertinent  works. 

The  study  of  finite  element  methods  (FEM)  to  predict  microelectronic  reliability  was 
initiated  by  Southland  (Ref.  1-1).  This  study  provided  documentation  and  guidelines  for 
the  application  of  FEA  to  perform  mechanical  and  thermal  analysis  of  microelectronic 
packages.  The  study  also  included  an  assessment  of  various  finite  element  codes  and  a  set 
of  guidelines  to  show  how  to  apply  FEA  to  electronic  devices.  The  STARDYNE  general 
purpose  finite  element  code  was  determined  to  be  the  most  applicable  for  analyzing 
electronic  packages. 

Computer  aided  engineering  (CAE)  software  which  predicts  the  vibration  response  of 
printed  circuit  boards  (PCB)  mounted  within  a  line  replaceable  unit  (LRU)  was  devel¬ 
oped  by  Dandawante,  Soovere  and  et  al.  (Ref.  1-2).  This  software  was  based  on  the  fi¬ 
nite  element  code  entitled  Numerically  Integrated  Elements  for  Structural  Analysis 
(NTS A).  This  report  demonstrates  finite  element  modeling  and  analysis  techniques  for 
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the  prediction  of  stresses  due  to  vibratory  loads.  NISA  was  previously  identified  by 
Southland  et  al.  as  having  the  capability  of  performing  both  vibration  and  thermal 
analysis.  The  advanced  thermal  analysis  capabilities  of  NISA  were  demonstrated  by 
Bivens  and  Bocchi  (Ref.  1-3).  Kallis  et  al.  (Ref.  1-4)  compared  the  capabilities  of  NISA 
to  NASTRAN  and  ANSYS  as  a  guideline  to  evaluate  analyses  that  are  done  by 
government  contractors  using  NASTRAN  and  ANSYS. 

Soovere,  Steinberg,  Dandawate,  and  et  al.  (Ref.  1-5),  developed  a  computer-aided 
design  process  to  predict  the  dynamic  stresses  induced  by  vibration  in  leaded  and  leadless 
chip  carriers  (LCC)  mounted  on  printed  circuit  boards.  They  demonstrated  how  these 
stresses  (calculated  using  NISA)  could  be  used  to  predict  the  fatigue  life  of  electronic 
components. 

Bivens  and  Bocchi  used  advanced  finite  element  simulations  to  see  how  various  physi¬ 
cal  factors  affect  surface  mounted  device  temperatures.  Detailed  finite  element  modeling 
techniques  for  PCBs  and  LCCs  were  demonstrated. 

Duncan  et  al.  used  finite  element  simulations  to  determine  the  stress  levels  in 
electronic  packaging  used  in  wafer  scale  integration  (WS1).  He  then  utilized  the  stress 
data  to  predict  reliability. 

The  transfer  of  finite  element  output  to  a  reliability  assessment  for  solder  connections 
was  demonstrated  by  Bivens  (Ref.  1-6).  Solomon.  Brzozowski  and  Thompson  (Ref.  1-7), 
using  3-D  FEA  showed  how  solder  joint  fatigue  life  could  be  predicted  using  joint 
strains  and  low  cyclic  fatigue  data. 

1'he  earliest  book  which  demonstrated  an  alternative  approach  to  FEA  for  predicting 
mechanical  reliability  of  microelectronic  packages  was  given  by  Steinberg  (Ref.  1-8). 
Steinberg  developed  a  "cook  book"  approach  f~r  the  'hock  and  vibration  analyses  of 
electronic  packages  based  upon  empirical  data  and  closed  form  solutions. 
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1.2  Universities 

Several  universities  are  widely  recognized  for  their  research  into  the  field  of  electronic 
life  predictions  and  electronic  material  characterizations.  The  following  two  are  the  most 
significant. 


1.2.1  The  University  of  Maryland 

The  University  of  Maryland  has  developed  a  program  called  Computer-Aided  Life 
Cycle  Engineering  (CALCE)  (Ref.  1-9).  This  software  package  performs  automatic  mod¬ 
eling  based  upon  inputs  by  the  user.  CALCE  provides  the  following  modeling  and  pre¬ 
diction  techniques: 

1)  Microelectronic  package  reliability  modeling  including  wire  and  wire  bond  failure 
modeling,  die  attach  failure  modeling,  corrosion  failure  modeling. 

2)  PCB  packaging  including  the  effect  of  components  on  PCB  natural  frequency,  de¬ 
sign  techniques  for  fatigue  of  surface  mount  technology,  coupled  thermal  and 
vibration  fatigue  analysis,  and  zonal  decomposition  techniques  for  thermal  analy¬ 
sis.  Currently,  the  vibration  analysis  is  based  on  Steinberg's  non-finite  element 
methods.  Finite  element  analysis  tools  for  future  versions  of  the  software  are  being 
developed  by  the  university.  Fatigue  life  predictions  are  based  on  S-N  data.  PCB 
assembly  packaging  analysis  includes  forced  convection  cooling  allocation  and  vi¬ 
bration  analysis  of  edge  conditions. 

3)  Interconnection  and  placement  theory  including  placement  for  producibility. 

4)  Experimental  measurement  techniques  including  infrared  experimental  setup  for 
PCB  analysis  and  infrared  studies  on  flaw  detection. 

5)  Design  matrix  for  system  evaluation  including  introducing  maintainability  into  the 
concurrent  design  process  and  a  user  definable  derating  system. 

6)  L  essons  learned  from  MIL-HDBK-217  studies. 

1.2.2  The  University  of  Wisconsin-Madison 

The  University  of  Wisconsin-Madison  has  an  established  Electronics  Packaging  and 
Interconnection  Research  Program  (Ref.  1-10).  The  following  three  areas  are  the  major 
thrusts  of  their  program: 
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1)  Isothermal  fatigue  testing  of  solder  joint  models  -  Evaluation  of  material  properties 
such  as  the  modulus  of  elasticity,  shear  modulus  and  Poisson’s  ratio  have  been  com¬ 
pleted  for  bulk  60/40  solder.  Current  research  employs  shear  specimens  to  develop 
life  prediction  methods  for  solder  joints.  The  methods  are  based  on  calculations  of 
damage  from  plastic  strain  and  creep.  In  order  to  accurately  account  for  damage,  life 
vs.  strain  amplitude  curves  have  been  obtained  for  several  temperatures. 

2)  SMT  Package  Test  Methodology,  Straddle  Board  Testing  -  This  area  concentrates 
on  measuring  the  stresses  and  strains  that  SMT  packages  experience  during  their 
lifetime.  Straddle  board  tests  permit  fatigue  life  comparison  between  package  de¬ 
signs.  Real  packages  are  tested  in  a  controlled  strain  environment.  Finite  element 
modeling  determines  the  stresses  and  strains  that  occur  within  the  solder  connec¬ 
tions  during  loading. 

3)  Expert  Systems  -  The  Expert  Systems  Software  combines  the  stress  and  strain  out¬ 
puts  from  finite  element  modeling  with  the  fatigue  life  data  from  the  tests  to  deter¬ 
mine  the  fatigue  life  of  actual  solder  joints. 

1.3  Industry 

The  following  paragraphs  briefly  describe  the  reliability  (life)  prediction  activities  by 
various  members  of  industry  which  are  involved  in  the  design  and  analysis  of  electronic 
equipment. 

1.3.1  AT&T 

Werner  Engelmaier  has  developed  many  of  the  tools  used  by  AT&T  Bell  Labs 
(Refs.  1-11,12,13).  These  tools  are  used  to  compute  the  thermal  expansion  mismatch 
between  the  board  and  the  chip  carrier  based  on  the  temperature  distribution  on  the 
printed  circuit  board.  This  mismatch  is  computed  with  the  following  non-finite  element 
relation: 


Ae  =  occ  (Tc  -  T0)  -  as  (Ts  --  T0) 


Eq.  1 
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where:  ac  is  the  coefficient  of  thermal  expansion  of  the  chip  carrier 
as  is  the  coefficient  of  thermal  expansion  of  the  substrate 
Tc  is  the  maximum  temperature  of  the  chip  carrier 
Ts  is  the  substrate  temperature 
T0  is  the  initial  temperature 
Ae  is  the  thermally  induced  expansion  mismatch 

AT&T  then  computes  a  shear  strain  from: 

Ay  =  Ae  Eq.  2 

'  Jl  h  M 

where:  L  is  the  length  of  the  chip  carrier 
h  is  the  height  of  the  solder  joint 
Ay  is  the  shear  strain  range 

Their  fatigue  life  computation  uses  the  Manson-Coffin  fatigue  life  relation  to  deter¬ 
mine  the  cycles  to  failure  based  on  the  strain  Ay.  Their  analysis  also  includes  models  of 
the  elastic-plastic  and  creep  behavior  of  solder. 

1.3.2  General  Dynamics 

General  Dynamics  primarily  uses  Steinberg’s  non-FEA  methods  to  predict  the  life  of 
electronics  (Ref.  1-14).  In  a  recent  study  conducted  for  the  US  Air  Force  Advanced  Thc- 
tical  Aircraft  (USAF-ATA),  they  applied  these  methods  in  the  analysis  of  a  circuit  board 
containing  surface  mounted  components  (Figure  1-2).  Two  edges  of  the  board  are  wedge 
clamped  to  the  sides  of  the  line  replaceable  unit  (LRU),  one  edge  is  unrestrained  and  the 
board  connection  is  on  the  other  edge.  They  focused  on  a  component  with  68  gull-wing 
leads  which  had  been  located  closest  to  the  unrestrained  edge  of  the  board.  Because  of 
its  location  and  size,  this  component  will  experience  large  vibration  displacements  and, 
therefore,  large  stresses  in  the  leads. 
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Printed  Circuit  Board 
(PCB) 

GP03-0262-2B 

Figure  1  -2.  Circuit  Board  Used  in  ATA  Analysis 

The  Steinberg  based  analysis  resulted  in  a  life  prediction  of  9,120  flight  hours  for  the 
component.  MCAIR  completed  a  parallel  effort  as  part  of  the  same  study,  but  used  finite 
element  methods  to  obtain  the  stresses  in  the  lead  and  the  solder  joints  of  the  same  criti¬ 
cal  component.  These  stresses  were  then  used  in  a  Miner’s  cumulative  damage  analysis 
to  predict  the  life  of  the  component.  MCAIR  predicted  a  life  of  117,905  flight  hours. 

This  comparison  indicates  that  the  Steinberg  method  results  in  a  more  conservative  life 
prediction  than  does  the  finite  element  techniques.  In  addition  to  this  work,  general 
Dynamics  has  developed  3-D  finite  element  models  of  leads  and  solder  connections. 

They  also  have  FEAP  and  the  University  of  Maryland's  CALCE  software. 
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1.3.3  Westinghouse 

Westinghouse  developed  a  FEA  model  for  analyzing  LCC  solder  joint  stresses  and 
strains  induced  by  thermal  or  mechanical  excursions  under  the  Air  Force  contract 
F33615-82-C-5047,  supported  by  the  Wright  Research  and  Development  Center 
(WRDC)  Materials  Laboratory.  The  Westinghouse  Electric  Computer  Analysis 
(WECAN)  calculates  the  strain  distribution  in  the  solder  joint  for  given  loading  condi¬ 
tions  allowing  the  user  to  study  the  levels  of  strain  and  their  location  (Ref.  1-15). 
Westinghouse,  a  member  of  the  CALCE  consortium,  is  incorporating  this  software  into 
its  tools. 

1.3.4  Martin  Marietta 

The  FEA  code  WECAN  (developed  by  Westinghouse)  was  used  by  Martin  Marietta  to 
obtain  Von  Mises  strain  in  solder  joints  (Ref.  1-7).  Solder  fatigue  data  was  obtained  in 
tests  documented  in  Reference  1-16,  work  sponsored  by  the  WRDC  Materials 
Laboratory,  Air  Force  contract  F33615-85-C-5065.  The  shear  strain  y  was  obtained 
from  Von  Mises  strain  e,™  by: 

y  =  /3  evm  Eq.  3 

These  strain  computations  were  then  used  to  predict  the  fatigue  life  of  LCC  solder 
joints.  Martin  Marietta  has  evaluated  the  effects  of  plasticity  and  crack  propagation  in 
the  solder  under  the  chip  carrier  and  across  the  fillet. 

Martin  Marietta  Electronic  Systems  in  Orlando,  FL  is  a  subcontractor  for  the 
USAF-ATA  As  part  of  this  work,  they  predicted  the  life  of  several  of  the  components  in 
their  designs  (Ref.  1-17).  Martin  Marietta  used  coarse  finite  element  models  of  circuit 
boards  to  determine  vibration  deflections  and  thermal  deformations.  This  information 
was  used  as  the  boundary  conditions  for  more  detailed  finite  element  models  of  the  chip 
carriers  mounted  on  the  PCB.  In  the  final  step,  detailed  models  of  the  solder  joints 
yielded  stresses  used  in  a  Miner’s  cumulative  damage  calculation  to  predict  the  creep / 
fatigue  life  of  the  components. 
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1.3.5  Texas  Instruments  (TI) 

The  TI  Materials  and  Control  Group  in  Attleboro,  MA,  has  used  the  ADINA  finite 
element  code  to  compute  3-D  strains  in  chip  mounts  (Ref.  1-18).  TI  in  Dallas  has  used 
ADINA  to  compute  stress  levels  inside  leaded  chip  carriers  (chip,  lead  frame,  chip  pad, 
mold  compound,  and  die  attach)  due  to  thermo-mechanical  loads  (Ref.  1-19). 

1.3.6  Hewlett-Packard  (HP) 

HP  has  completed  FEA  based  thermal  stress  analyses  of  gullwing  leads,  their  connec¬ 
tion  to  the  chip  carrier  and  the  solder  joint  connection  to  the  board  (Refs.  1-20,21).  This 
work  included  evaluations  of  solder  joint  quality.  HP  has  also  completed  FEA  of  J-leads, 
fatigue  tests  and  statistical  analyses  of  the  test  data.  The  tests  involved  boards  loaded  in 
cyclic  four  point  bending  until  solder  joint  failures  occurred. 

1.3.7  Litton 

The  current  method  used  at  Litton  is  Steinberg’s  non-FEA  technique  supplemented 
with  finite  element  analyses  of  the  board  temperature  distributions  and  board  vibration. 
Although  Litton  is  fully  capable  of  generating  detailed  finite  element  models  of  leads  and 
solder  connections,  their  studies  have  shown  that  inaccuracies  occur  when  using  finite  ele¬ 
ment  analyses  for  stress  and  strain  calculations  in  very  small  objects.  FEA  accuracy  will 
be  discussed  in  more  detail  later  in  the  report.  Because  of  the  inaccuracies  they  found, 
Litton  has  concentrated  on  hand  calculations  to  determine  stresses  in  materials  and  pre¬ 
dict  their  fatigue  lives  (Ref.  1-22).  Thermal  stresses  are  computed  with  strength  of  mate¬ 
rials  equations.  Vibration  analyses  are  based  on  resonant  frequency  equations  modified 
with  empirical  factors.  Fatigue  life  is  then  computed  by  adding  the  damage  caused  by  vi¬ 
bration  and  by  thermal  cycles. 

1.3.8  Pacific  Numerix 

Pacific  Numerix  Corporation  develops  and  markets  specialized  computer  aided  engi¬ 
neering  software  for  the  electronics  industry.  Pacific  Numerix  programs  are  used  in  the 
design  and  verification  process  to  analyze  the  effects  that  parts,  materials,  component 
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placement  and  environmental  conditions  have  on  the  performance,  manufacturability, 
reliability,  testability  and  maintainability  of  printed  circuit  boards. 

Pacific  Numerix  programs  comprise  the  PCB  Design  Expert  System™  (PDES)  and 
include  (Refs.  1-23,24): 

1)  PCB  Place  -  Manual  and  automatic  placement  program 

2)  PCB  Thermal  -  Finite  element  printed  circuit  board  thermal  analyzer 

3)  PCB  Vibration  -  Finite  element  printed  circuit  board  vibration  analyzer 

4)  PCB  Fatigue  -  Finite  Element  Printed  Circuit  Board  Stress/Fatigue 
Analysis  Program 

The  board  level  thermal  and  vibrational  results  obtained  from  running  PCB  Thermal 
and  PCB  Vibration  are  automatically  converted  into  appropriate  boundary  conditions  and 
are  applied  to  the  detailed  finite  element  models  for  stress  analysis  at  the  component  lev¬ 
el.  PCB  Fatigue  performs  thermal  and  vibrational  FEA  on  critical  detailed  regions  such 
as  die  bonds,  bonding  wires,  solder  posts,  leads,  plated  thru  holes  and  vias.  PCB  Fatigue 
also  has  the  capability  of  analyzing  the  effects  of  parts  tolerance  and  manufacturing  pro¬ 
cess  variability  on  the  life  of  the  equipment.  The  program  automatically  generates  de-- 
tailed  finite  element  models  for  a  variety  of  PCB  assembly  components. 

1,3.9  Engineering  Mechanics  Research  Corporation  (EMRC) 

EMRC  has  designed  a  program  called  Finite  Element  Analysis  of  PCBs  (FEAP)  which 
allows  a  structural  engineer  to  model  and  evaluate  PCBs  and  critical  components  without 
conducting  expensive  experiments.  FEAP  is  able  to  predict  temperature  distributions 
across  a  PCB,  calculate  stresses  in  components,  leads  and  PCBs  due  to  vibration  and 
thermal  loads,  and  estimate  the  life  of  critical  components. 

FEAP  (Figure  1-3)  has  a  sophisticated  built-in  graphics  capability  which  is  geared  to¬ 
wards  helping  the  user  define  the  PCB.  Once  the  PCB  and  the  components  have  been  de¬ 
fined,  FEAP  generates  a  coarse  finite  element  mesh  for  the  entire  board  which  includes 
modeling  its  critical  components  as  plates  and  beams.  After  the  finite  element  mesh  of  the 
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PCB  and  components  is  automatically  generated,  either  heat  transfer  or  dynamic  analysis 
can  be  performed.  Once  the  vibration  or  thermal  analysis  of  the  PCB  has  been  completed, 
the  forces,  moments,  displacements  and  rotations  of  the  PCB  are  available  at  a  preselected 
boundary  on  the  PCB  surrounding  the  electronic  component  selected  for  fatigue  analysis.  A 
much  finer  PCB  mesh  is  then  automatically  generated  for  this  component  and  a  cubic  inter¬ 
polation  is  used  to  obtain  a  more  detailed  distribution  of  nodal  forces  and  moments  along 
the  selected  boundary.  This  unique  sub-model  approach  improves  the  accuracy  of  the  solu¬ 
tion  within  a  defined  area  of  the  PCB  by  applying  the  displacements/loads  from  the  vibra¬ 
tion  or  thermal  analysis  as  statically  imposed  enforced  displacements  or  loads.  Further 
refinement  can  then  be  performed  by  generating  a  solid  finite  element  model  of  the  compo¬ 
nent  and  imposing  the  displacements/loads  from  the  refined  mesh  analysis.  Finally,  built-in 
S-N  fatigue  curves  are  used  to  relate  the  FEA  outputs  (stress/strain)  to  component  life 
(number  of  cycles  to  failure  or  lifetime)  for  common  materials  found  in  electronic  assem¬ 
blies.  The  fatigue  life  is  predicted  by  searching  the  finite  element  model  for  regions  of  criti¬ 
cal  stress/strain  and  then  using  the  appropriate  cycles  to  failure  curves  (Figure  1-4)  for  the 
material  with  the  high  stress  concentration. 

FEAP  Is  Comprised  of  and  Totally  Integrates  the  Following  Analysis: 


FEAP  Is  Able  to: 

•  Predict  Temperature  Profiles 

•  Calculate  Stresses  in  Components,  Leads,  PCBs 

•  Estimate  Life  of  Critical  Components  GP030007-275.D 

Figure  1  -3.  Finite  Element  Analysis  of  PCBs  (FEAP) 
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PCBs,  Chips 


Perform  Thermal  or 
Vibration  Analysis 


Figure  1  -4.  FEAP  Life  Prediction 


1.3.10  Lockheed 

As  part  of  a  USAF  -  WRDC  project  (Contract  Number  F33615-85-C-3403), 
Lockheed  has  developed  a  finite  element  based  reliability  analysis  in  combination  with 
EMRC  and  Steinberg  and  Associates  (Ref.  1-5).  The  analysis  uses  the  NISA  finite  ele¬ 
ment  code  with  a  pre-processor  to  model  the  geometries  of  the  electronics  and  generate 
the  finite  element  mesh.  A  post-processor  displays  the  output  stresses  as  contour  plots 
which  highlights  the  critical  regions. 

The  output  stress  is  then  used  to  predict  the  cycles  to  failure  by  using  stress  versus  life 
data  for  the  given  material.  The  life  analysis  model  holds  either  vibration  or  thermal 
stresses  constant  while  using  the  other  to  compute  the  life.  Only  linear  material  behavior 
is  considered  in  the  analysis.  (Non-linear  material  behavior  can  occur  if  the  yield 
.  trength  of  the  material  is  exceeded  and  plastic  deformation  occurs.) 
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1.3.11  McDonnell  Douglas  Corporation 

MCAIR  has  applied  FEA  and  fatigue  life  predictions  to  the  Air  Force  version  of  the 
Advanced  Tactical  Aircraft  (ATA)  Flight  Control  Computer  (Ref.  1-25).  The  fundamental 
frequency  of  the  PCB  and  dynamic  displacements  were  calculated  with  NASTRAN.  The 
displacements  were  then  imposed  on  detailed  models  of  the  leads  and  solder  joints 
(Figure  1-5)  to  compute  v  oration  stresses.  A  similar  analysis  was  completed  to  find  the 
thermal  stresses.  The  combined  thermal-vibration  stress  history  was  then  used  in  a 
Miner's  Rule  analysis  to  predict  fatigue  life.  Solder  joint  geometry  and  solder  quantity 
were  varied  to  evaluate  the  impact  of  initial  quality  on  stresses  and  life. 


MCAIR  is  also  developing  simplified  non-FE  procedures  to  predict  life  based  on 
stress  computations  using  strength  of  materials  equations.  Miner's  rule  is  used  for  life 
predictions.  These  procedures  are  being  automated  in  the  Fortran  program  “Electronics 
Life  Prediction”  (ELIFE).  This  code  will  serve  as  a  quick  analysis  tool  to  isolate  prob¬ 
lem  areas  which  may  require  subsequent  detailed  finite  element  analyses.  The  existing 
capability  of  ELIFE  includes  the  following  five  analyses: 

1)  Thermal  stress  analyses  of  LCCs,  leaded  surface  mount  devices,  dual-in-line 
packages  (DIP)  and  plated  through  holes  (PTH). 

2)  Thermal  fatigue  life  computations  for  components  in  item  (1). 

3)  PCB  vibration  analysis. 

4)  Vibration  stress  analysis  of  components  in  item  (1). 

5)  Life  predictions  under  combined  thermal  and  vibration  stresses. 


1-13 


Computer  Aided  Assessment  of  Reliability 

Using  Finite  Element  Methods - - - — — — — - 

Later  versions  of  ELIFE  will  include  a  heat  transfer  analysis  to  calculate  tempeiaiure 
distributions  and  an  expanded  component  and  material  library. 

McDonnell  Douglas  Electronics  Systems  Company,  in  combination  with  MCAIR,  has 
tested  printed  circuit  boards  to  obtain  solder  S-N  fatigue  data.  These  tests  subjected 
surface  mounted  LCCs  to  thermal  cycling. 

1.3.12  Hughes  Aircraft 

The  Electro-Optical  and  Data  Group  within  Hughes  Aircraft  is  currently  developing 
methods  for  analyzing  cracks  in  avionics  under  the  U.S.  Air  Force  contract  number 
F33615-87-C-3403  (Ref.  1-26,27).  As  part  of  this  project,  they  evaluated  failure  sources 
in  the  Hughes  APG-63  radar  used  in  the  F-15  aircraft.  They  have  also  fatigue  tested 
wires  in  vibration,  and  thermally  cycled  wire  bonds  and  plated  through  holes.  Further 
tests  will  involve  electronic  assemblies  subjected  to  thermal  cycles. 

Failure  Analysis  Associates  has  completed  FEAs  for  Hughes  to  determine  the  stresses 
and  strains  in  wire  bonds  exposed  to  mechanical  loads.  The  calculated  strains  were  simi¬ 
lar  to  strains  measured  using  stereoimaging  methods  developed  at  Southwest  Research 
Institute.  Southwest  Research  is  also  supporting  Hughes  with  fracture  mechanics  analysis 
of  cracks  in  electronics.  The  same  group  at  Hughes  has  completed  an  FEA  of  wafer  scale 
integrated  devices  as  part  of  the  USAF  contract  number  F30602-87-C-0118  (Ref.  1-24). 
They  compared  NISA.  ANSYS,  and  NASTRAN  finite  element  codes  and  developed  a 
computer  program  which  uses  finite  element  output  stresses  to  compute  fatigue  life.  This 
program  contains  a  material  properties  database  which  covers  low  cycle  fatigue  (strain  vs. 
cycles  to  failure),  high  cycle  fatigue  and  crack  growth  data.  The  life  computation  does 
not  account  for  mean  stress  effects  which  may  occur  when  vibration  and  thermal  stresses 
occur  simultaneously  or  when  the  temperature  cycle  does  not  result  in  fully  reversed 
strains. 
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i.4  U.S.  Government  Agencies 

1.4.1  National  Institute  of  Standards  and  Technology 

The  U.S.  Department  of  Commerce  National  Institute  of  Standards  and  Technology 
sponsors  electronic  packaging  workshops  to  establish  material  properties  and  measure¬ 
ment  techniques  (Ref.  1-28).  Areas  where  their  efforts  are  focused  include  test  methods, 
mechanical  properties  and  thermal  properties.  The  materials  of  interest  cover  most  elec¬ 
tronic  applications  and  include  polymers,  fiberglass,  ceramics  and  metals.  This  type  of 
data  is  essential  for  finite  element  based  reliability  assessments. 

1.4.2  Wright  Laboratory  (WL)  (Formally  Wright  Research  and 
Development  Center  (WRDC)) 

The  WL  Environmental  Control  Branch  compared  finite  element  vibration  analyses 
with  tests  and  reported  their  findings  in  WRDC-TR-89-3110  entitled  “Correlation/Vali¬ 
dation  of  Finite  Element  Code  Analyses  for  Vibration  Assessment  of  Avionic  Equip¬ 
ment”  (Ref.  1-29).  Plates,  representing  PCBs,  were  tested  dynamically  with  various 
boundary  conditions.  Capacitors,  diodes,  transistors,  relays  and  dips  were  mounted  on 
the  boards.  The  vibration  frequencies  of  the  boards  were  compared  with  calculated  fre¬ 
quencies  using  NISA  and  NASTRAN  finite  element  codes.  The  outputs  of  NISA  and 
NASTRAN  were  within  5%  of  the  measured  vibration  frequencies. 

WRDC  has  also  sponsored  various  contracted  efforts  in  this  area  including  the  pro¬ 
grams  documented  in  References  1-2,  1-5,  1-7,  1-16,  1-26  and  1-27. 

1.43  Rome  Laboratory  (RL)  (Formally  Rome  Air  Development  Center  (RADC)) 

Extensive  finite  element  analyses  for  electronics  reliability  assessments  have  been  con¬ 
ducted  at  RL,  Griffiss  AFB,  NY.  (Refs.  1-3,6).  Report  RADC-TR-87-177,  entitled  “Re¬ 
liability  Analyses  of  a  Surface  Mounted  Package  Using  Finite  Element  Simulation", 
documents  efforts  to  analyze  PCBs  and  electronic  package  assemblies  to  simulate  the 
effects  of  die  size,  heat  producing  areas,  voids  in  the  die  attach  and  thermal  undercoat  on 
the  package  thermal  resistance  using  the  NISA  finite  element  code.  The  thermal 
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computations  from  the  package/PCB  analysis  were  then  used  in  a  3-D  finite  element 
model  of  a  leadless  solder  connection  to  compute  thermally  induced  stresses  for  both 
models. 

In  another  study  (Ref.  1-6),  3-D  NISA  models  were  constructed  of  systems  containing 
the  chip  carrier,  circuit  board  and  connections.  The  models  included  leadless,  S-lead  and 
gull-wing  lead  systems.  Thermally  induced  deformations  obtained  with  these  models 
were  then  imposed  on  more  detailed  finite  element  models  of  the  leads  and  solder  con¬ 
nections.  The  output  stresses  were  then  used  to  estimate  the  strain  amplitudes  experi¬ 
enced  during  the  thermal  cycles.  Life  predictions  were  then  obtained  from  material 
curves  of  strain  versus  cycles  to  failure. 

RADC  also  sponsored  projects  documented  in  report  RADC-TR-82-133,  entitled 
"Finite  Element  Analysis  of  Microelectronic  Packages”  (Ref.  1-1)  and  efforts  by  Hughes 
Aircraft  under  contract  number  F30602-87-C-0018,  entitled  “Reliability  Assessment  of 
Wafer  Scale  Integration  Using  Finite  Element  Analysis”  (Ref.  1-4). 

1.5  Finite  Element  Code  Suppliers 

Thble  1-1  lists  the  suppliers,  the  name  of  the  finite  element  code  and  the  purpose  of  the 
code  for  state-of-the-art  finite  element  codes  presently  used  within  industry  for  reliability 
(life)  analysis.  All  of  these  can  be  used  to  predict  stress  and  strain  levels  in  electronic  mm. 
ponents:  however,  the  ability  to  do  automated  packaging  analysis  from  a  reliability  stand¬ 
point  is  limited  to  the  Pacific  Numerix  and  EMRC  codes.  FEAP  automates  and  integrates 
the  capability  to  do  stress  analysis  and  life  prediction  into  one  package. 

The  matrix  shown  in  Table  1-2  summarizes  the  applicability  of  various  finite  element 
and  non-finite  element  techniques  to  fundamental  electronic  reliability  analyses.  Addi¬ 
tionally,  Table  1-3  identifies  the  unique  characteristics  of  the  FEA  codes  as  they  apply  to 
structural  mechanics  applications,  heat  transfer  applications  and  to  fatigue  and  fracture 
mechanics  applications. 
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TABLE  1-1.  FINITE  ELEMENT  CODE  SUPPLIERS 


Coda 

Supplier 

Purpose 

ABAQUS 

Hibbit,  Karlsson,  and  Sorensen  Inc. 

Commercial/General 

ADINA 

ADINA  R&D  Inc. 

Commercial/Genera! 

ANSYS 

Swanson  Inc. 

Commercial/General 

IDEAS 

Structural  Dynamics  Research  Corp. 

Commercial/General 

MSC/NASTRAN 

MacNeal-Schwendler  Corp. 

Commercial/General 

FEAP'NISA  II 

Engineering  Mechanics  Research  Corp. 

Commercial/General 

PDES 

Pacific  Numenx  Corp. 

Commercial/Electromc 

Packaging 

PROBE 

MacNeal-Schwendler  Corp. 

(formally  of  Noetic  Technologies) 

Co  mmercial/Stress 
Analysis 

WECAN 

Westinghouse  R&D  Center 

Commercial/Electronic 
Component  Analysis 

GP03-0007-273-D 


TABLE  1-2.  CODES  AND  METHODS  MATRIX 


Complexity  of 
the  Analysis 

1  mile  Element 
Uased 

Strength  of 

Material 

Specific  to 
Electronic 
Application 

Representing  the  Failure 
Mechanism 


Input  Data 
Characteristics 


Computational  .  ..  .. 

Accuracy  Efliciency  AppUcat.ons 
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TABLE  1-3.  FINITE  ELEMENT  APPLICATION 


tmm/fim 


m/mMM/ 


Structural  Mechanics  Application 


Haal  Transtar  Application 


Types  of  Elements 

•  3-D  Rod 

•  3-0  Beam 

•  3iane  Stress 

•  Plane  Strain 

•  Membrane 

•  Skear  Pane-.. 

•  3la:es 

•  7h;n  Sheds 

■ "me*  Sheos 
•3-0  Sends 

•  J.screte  Si; “ness 

•  Jeer’d  ary 

•  3ap 

Range  of  Applications 

•  .  near  Statics 

•  E.genvaiue 

-  =ree  Vibration 

-  3_CKling 

•  \cnnnear  Statics 

-  Resoonse 

-  Postbcckling 

•  \ommear  Oynamics 

-  Nonlinear  Vibration 

-  T  nansient 

•  Nonnnear  interaction 

-  Fluid-Structural 

-  Tiermai-Vlecnanicai 

•  Plasticity 

-  Stat’C 

-  Dynamic 

Solution  Tecfiniques 

•  Nonnnear  Stat  cs 

-  incremental 

-  Newton-Tyoe 

•  Non  near  Dynamics 

-  Modei  Sooerocs.tion 

-  D:-ect  integration 

Explicit 

'mDliCit 

Types  of  Loading 

•  Concentrated 

•  Line 

•  ’normal 

•  Random 

•  Deformation  Deoendent 

•  'nitial  Stress 

Support  Conditions 

•  P-escroed  Disoacement 

•  Elastic  Founaat  on 


Space  Dimensionality 

•3-0 

•2-0 

Range  ot  Applications 

•  Linear  Steady  State 

•  Nonlinear  Steady  State 

•  Linear  Transient 

•  Nonlinear  Transient 

Material  Properties 

•  isotrcoic 

•  Anisotropic 
•Multilayered 

•  Temperature  Deoendent 
•Time  Dependent 

Boundary  Conditions 

•  Prescribed  Temperature 

•  Convection  From  a  Surface 

•  Radiation 

Other  Capabilities 

•  Thermal-Stress  Coupling 


XXX 

x  x 
xxx 
x  x 


x  x  x  x 

x  x  x  x 

X  X  X  X 


Fatigue  and  Fracture  Mechanics  Application 


Types  ol  Elements 

•  Plane  Stress 

•  Plane  Strain 

•  Plates 

•  Thin  shells 

•  Thicx  Shells 

•  3-0  Solids 

Range  of  Applications 

•  Linear  Statics 

•  Nonlinear  Statics 
•Nonlinear  Dynamics 

Loading 

•  Static 

•  FuncDon  of  Time 

•  Mode  l 

•  Mode  II 

•  Mode  ill 

•  Combines  Modes 

•  initial  Stress 

Surface  Crack  Geometry 

•Straight 

•Curve 


XXX 
XXX 
X  X 


Material  Properties 

•  Sbtrcoic 

•  Arsotrcpic 

•  Multilayered 
•’empe'ature  Dependent 

•  E  astic- Strain  hardening 
•Viscoelastic 

•  Nonnnear  Elastic 

•  Viscoplastic 

•  h.gh  Tempera!'jre  Creep 


GPC3  0624-72-0 


1-18 


Computer  Aided  Assessment  of  Reliability 
_ Using  Finite  Element  Methods 


References 

[1-1]  Southland,  J.R.  Beatty,  V.R.,  and  Vitaliano,  W.J.,  “Finite  Element  Analysis  of 
Microelectronic  Packages,”  RADC-TR-82-133,  May  1982. 

[1-2]  Soovere,  J.,  Dandawate,  B.V;  “Vibration  Stress  Analysis  of  Avionics”, 

AF WAL-TR-87-3 023 ;  April  1987. 

[1-3]  Bivens,  G.A,  and  Bocchi,  W.J.,  “Reliability  Analyses  of  a  Surface  Mounted 
Package  Using  Finite  Element  Simulation,”  RADC-TR-87-177,  October  1987. 

[1-4]  Kallis,  J.M.,  Duncan,  L.B.,  Van  Westerhuyzen,  D.H.,  Sandkulla,  D.C.,  MacFarlane, 
J.D.,  Naepflin,  H.P.,  Ingersoll,  J.G.,  Tierney,  B.D.,  McHorney,  P.E.,  Erickson,  J.J.,  and 
Terrill,  K.W.,  “Reliability  Assessment  of  Wafer  Scale  Integration  Using  Finite  Element 
Analysis",  Hughes  Aircraft  Company  under  RADC  contract  number  F30602-87-C-0118, 
February  1990. 

[1-5]  Soovere,  J.,  Dandawate,  B.V.,  Garfinkel,  G.A,  Isikbay,  N„  and  Steinberg,  D.S.. 
“Vibration  Reliability  Life  Model  for  Avionics,”  AFWAL-TR-87-3048. 

[1-6]  Bivens,  G.A,  “Reliability  Assessment  Using  Finite  Element  Technique,” 
RADC-TR-89-28 1 ,  Nov  1989. 

[1-7]  H.D.  Solomon,  V.  Brzozowski,  and  D.C.  Thompson.  “Prediction  of  Solder  Joint 
Fatigue  Life,”  AFWAL-TR-89-4002,  April  1988. 

[1-8]  Steinberg,  D.S.,  “Vibration  Anabasis  for  Electronic  Equipment,"  John  Wiley  and 
Sons,  1973. 

[  1  -9]  Conversation  with  Michael  Pecht,  Director  of  CALCE,  University  of  Maryland. 
1989. 

[1-10]  Conversation  with  Bela  Sandor,  Director  of  the  Electronics  Packaging  and 
Interconnection  Program,  University  of  Wisconsin-Madison,  1989. 


1-19 


Computer  Aided  Assessment  of  Reliability 

Using  Finite  Element  Methods  - 

[1-11]  Engelmaier,  W.,  “Reliability  of  Surface  Mounted  Assemblies:  Controlling  the 
Thermal  Expansion  Mismatch  Problem,”  Presented  at  the  SMART  IV  Conference, 

Los  Angeles,  CA,  January  11-14,  1988. 

[1-12]  Engelmaier,  W.,  “Fatigue  Life  of  Leadless  Chip  Carrier  Solder  Joints  During 
Power  Cycling,”  Proceed'ngs  of  the  Technical  Program  of  the  2nd  Annual  International 
Electronics  Packaging  Society  Conference,  San  Diego,  CA,  November  15-17,  1982. 

[1-12]  Engelmaier,  W„  “A  Reliability  Test,”  Circuits  Manufacturing,  June  1988. 

[1-14]  Conversation  with  Jeff  Intrieri,  General  Dynamics  Fort  Worth  Division,  1989. 

[1-15]  WECAN  -  Westinghouse  Electric  Computer  Analysis,  User’s  Manual, 
Westinghouse  R&D  Center,  1310  Beulah  Road,  Pittsburgh,  PA,  15235. 

[1—16]  Soloman.  H.D.,  “Low  Cycle  Fatigue  of  Surface  Mounted  Chip  Carrier/Printed 
Wiring  Board  Joints”,  AFWAL-TR-89-4Q04,  September  1987. 

[1-17]  Conversation  with  Bruce  C.  Steffens,  Staff  Stress  Engineer,  Martin  Marietta 
Electronic  Systems,  1990. 

[1-18]  Subramanyam,  S.  and  Rolph,  W.D.,  “Thermal  Analysis  of  Direct  Chip  Mount 
Electronic  Systems  Using  ADINA-T,”  Computers  &  Structures,  Vol.  32,  No.  3/4,  pp. 
853-859,  1989. 

[1-19]  Schroen,  W.H.,  Blanton,  P.S.,  and  Edwards,  D.R.,  “Finite  Element  Analysis  of 
Semiconductor  Devices,”  Presented  at  the  9th  Annual  IEE/AESS  Dayton  Chapter 
Symposium,  30  November  1988. 

[1-20]  Lau,  J.H.  and  Harkins.  C.G.,  “Thermal  Stress  Analysis  of  SOIC  Packages  and 
Interconnections,”  IEEE  Transactions  on  Components,  Hybrids,  and  Manufacturing 
Technology,  Vol.  11,  No.  4,  December  1988. 

[1-21]  Lau,  J.H.,  Harkins,  C.G.,  Rice,  D.,  Krai,  J.  and  Wells,  B.,  “Experimental  and 
Statistical  Analysis  of  Surface-Mount  Technology  PLCC  Solder-Joint  Reliability,"  IEEE 
Transactions  on  Reliability,  Vol.  37,  No.  5,  December  1988. 


1-20 


Computer  Aided  Assessment  of  Reliability 
- Using  Finite  Element  Methods 


[1-22]  Steinberg,  D.S.,  “Tools  Available  for  Implementing  AVIP,”  Presented  at  the  9th 
Annual  IEEE/AESS  Dayton  Chapter  Symposium,  30  November  1988. 

[1-23]  Company  Literature,  Pacific  Numerix,  1200  Prospect  Street,  Suite  300, 

La  Jolla,  CA,  92037,  (619)  587-0500. 

[1-24]  Conversation  with  Mark  Tinges,  Pacific  Numerix,  1989. 

[1-25]  "USAF  ATA  Initial  Avionics  Integrity  Assessment,”  McDonnell  Aircraft 
Company/General  Dynamics  Fort  Worth  Division,  Doc.  No.  86PR0579,  18  December 
1989. 

[1-26]  Kallis,  J.M.,  "Electronics  Reliability  Fracture  Mechanics”,  Presentation  at  the  2nd 
Workshop  on  the  Application  of  Fracture  Mechanics  to  Microscale  Structures  of  Modern 
Electronics,  El  Segundo,  CA,  May  1990. 

[1-27]  Burkhard,  A.H.,  Kallis,  J.M.,  Duncan,  L.B.,  Kanninen,  M.F.,  and  Harris,  D.O., 
"Application  of  Fracture  Mechanics  to  Microscale  Phenomena  in  Electronic  Assemblies." 
Proceedings  of  the  7th  International  Conference  on  Fracture  (ICF7),  Houston,  TX.  20-24 
March  1989. 

[1-28]  Conversation  with  Dr.  David  T.  Read,  US  Department  of  Commerce,  National 
Institute  of  Standards  and  Technology,  Materials  Reliability  Division,  Boulder,  CO. 

[1-29]  Bhungalia,  A.,  and  Kurylowich,  G„  "Correlation/Validation  of  Finite  Element 
Code  Analyses  for  Vibration  Assessment  of  Avionic  Equipment,”  WRDC-TR-89-31 10, 
December  1989. 


1-21 


Computer  Aided  Assessment  of  Reliability 
_ Using  Finite  Element  Methods 


Chapter  2 

ELECTRONIC  MATERIAL  PROPERTIES 
2.0  Introduction 

The  purpose  of  this  chapter  is  to  identify  the  physical  properties  of  materials  common¬ 
ly  used  in  the  production  of  electronic  hardware.  These  properties  are  necessary  inputs 
to  finite  element  analyses  used  to  calculate  thermal  distributions,  relative  deformations 
and  environmentally  induced  stresses  and  strains.  These  properties  fall  into  the  following 
two  classifications. 

1)  Thermal  -  The  primary  properties  in  this  classification  are  thermal  conductivity 
and  the  coefficient  of  thermal  expansion  (CTE).  Where  appropriate,  the  tempera¬ 
ture  at  which  a  material  changes  state  was  identified  if  these  temperatures  can  be 
expected  in  the  normal  operating  range  of  electronics. 

2)  Mechanical  -  This  classification  encompasses  numerous  properties  including  ten¬ 
sile  strength,  shear  strength,  modulus,  density,  Poisson’s  ratio,  stress  cycles  to  fail¬ 
ure  curves,  shear  stress  cycles  to  failure  curves,  strain /effective  strain/shear  strain 
cycles  to  failure  curves,  stress  vs.  strain  curves,  etc. 

For  the  purposes  of  this  Chapter,  a  failure  mechanism  is  defined  as  the  property  by 
which  a  material  will  fail.  These  mechanisms  include  ductile  fracture,  brittle  fracture, 
fatigue  and  buckling.  A  failure  mode  is  the  identifiable  result  of  the  failure  mechanism, 
ie.,  a  cracked  solder  joint.  Failure  mechanisms,  which  can  be  addressed  by  finite  element 
methods,  appropriate  for  each  material  were  identified.  The  failure  modes  of  electronics, 
based  upon  the  failure  mechanisms  of  the  material,  were  then  identified.  Based  upon  the 
preceding  definitions  and  the  scope  of  this  report,  failure  mechanisms  induced  by  chemi¬ 
cal  reactions,  processing  defects  and  electrical  overstress  were  not  identified. 

2.1  Failure  Mechanisms  and  Failure  Modes 

Figure  2-1  is  an  illustration  of  typical  line  replaceable  unit  (LRU).  The  failure  mecha¬ 
nisms  and  failure  modes  have  been  identified  for  the  printed  circuit  board  (PCB),  the 
housing  (or  chassis),  the  component  and  the  component  to  PCB  interface.  In  addition. 
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the  type  of  analyses  which  finite  element  methods  can  perform  on  each  part  of  the  LRU 
has  been  identified  under  the  heading  “Other”.  Table  2-1  identifies  the  environments 
and  forces  which  are  responsible  for  the  failure  mechanisms.  Table  2-2  identifies 
material  characteristics  and  the  type  of  FEA  analyses  that  should  be  performed  on  the 
material  based  on  their  application,  characteristics  and  failure  mechanisms. 


Printed  Circuit  Boards  (PCB) 

•  Failure  Mechanisms 

-  Buckling 

-  Fracture  -  Ductile 

-  Brittle 

-  Fatigue 

•  Failure  Modes 

-  Cracks  -  PCB 

-  Plated  Thru  Holes 

-  Heat  Sink 

-  Connector 

-  Traces 

-  Delaminations  -  Solder  Pad 

-  PCB 

-  PCB-PCB  Interference 


Component 

•  Failure  Mechanisms 

-  Fracture  -  Brittle 

-  Fracture  -  Ductile 

-  Fatigue 

-  Buckling 

•  Failure  Modes 

-  Buckled  Leads 

-  Cracked  Leads 

-  Cracked  Bond  Wire 

-  Cracked  Case 

-  Delaminated  Chip 

-  Cracked  Chip 


PCB  -  Component  Interface 

•  Failure  Mechanisms 

-  Fracture  -  Ductile 

-  Fatigue 

•  Failure  Mode 

-  Cracked  Solder 

•  Other 

-  Thermal  Conductance 


•  Other 

-  Natural  Frequency 

-  Displacements 

-  Dynamic  Modes 

-  Thermal  Conductance 


Housing  Structure 

•  Failure  Mechanisms 

-  Buckling 

-  Fatigue 

-  Fracture  -  Ductile 

•  Failure  Modes 

-  Cracked  Structure 

-  Warping,  Deformation 

-  Delamination 


Other 

-  Thermal  Conductance 

-  Dynamic  Transmissibility 

-  Natural  Frequency 


GP03-0624-70-dpt 


Figure  2-1.  “Black  Box”  Failure  Modes  &  Mechanisms 
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TABLE  2-1.  PREDOMINANT  FAILURE  MODES 


TABLE  2-2.  MATERIAL  APPLICATION  AND  FAILURE  MODES 


Material 

Failure  Mechanism 

Characteristics 

Application 

Failure  Mode 

Metals 

Rupture,  Fatigue 
Buckling 

Ductile 

Housing 

Cracks,  Warping, 
Permanent  Distortion 

PCB  Thermal  Plate 

Delamination  From 
Dielectrics,  Interference  With 
Adjacent  PCB,  Cracks 

Leads 

Cracking,  Distortion 

Plated  Thru 

Holes 

Cracking 

Solder 

Rupture,  Fatigue 

Ductile,  Creep 

Chip  to  PCB  Interface, 
Die  Attachment 

Cracking 

Organic  Resin 

Rupture,  Buckling 

Brittle 

PCB 

Delamination  From 
Conductors,  Thermal 
Plane;  Distortion; 
Cracking 

Ceramics 

Rupture 

Very  Brittle 

Chip  Housing,  PCB 

Cracking 

Semiconductor 

Rupture 

Very  Brittle 

Integrated  Circuits 

Cracking 
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2.2  General  Material  Properties 

This  section  contains  the  various  material  properties  necessary  to  perform  finite 
element  analyses.  Tables  2-3  and  2-4  are  a  collection  of  many  materials  and  properties. 
The  data  in  Table  2-3  was  collected  from  many  sources,  while  Table  2-4  is  an  excerpt 
from  Reference  2-1.  Data  in  these  tables  includes  both  thermal  and  mechanical 
properties.  Ranges  are  indicative  of  the  values  reported  by  different  sources,  a  result  of 
the  literature  search  approach  to  data  collection.  Table  2-5  is  a  listing  of  the  CTEs  of 
various  materials.  Figure  2-2  illustrates  CTE  information  in  a  manner  which  allows  rapid 
identification  of  materials  with  similar  values.  Table  2-6  (Ref.  2-2)  contains  data  on  the 
numerous  types  of  laminates  and  heat  sinks  used  in  PCBs.  Information  includes  the  glass 
transition  temperature,  the  CTE  in  all  planes,  thermal  conductivity,  the  elastic  modulus 
and  the  density. 
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TABLE  2-3.  GENERAL  MATERIAL  PROPERTIES 


Materiel 

Thermal  Conductivity 

TCE 

Tensile 

Yield 

Modulus  of 

Elongation  % 

Density 

8TU/hr*  ft'deg  f 

parts  per 

Strength 

Strength 

elasticity 

in/  in 

lb/in.*3 

X-T  Axis  2  Axis 

million 

(KPSI) 

(XPSI) 

(*10*6) 

Air 

0.017 

Acrylic  Conformal  Coat 

0.1 

50-90 

2.5 

2 

Alloy  42 

8.8 

4-5.3 

68 

35 

21 

43 

0.291 

A  l  uni  nun 

99-128 

10.4-23.8 

32 

28 

10 

15 

0.098 

Be 

87 

11.5 

42 

0.067 

BeO  (95%) 

120 

6. 4-7. 2 

13.5 

*24.18 

39 

0.036 

0.103 

8oron<  20%)  -  A  t  uni  nun 

100 

Carbon  Fiber  Epoxy 

-.5-  2 

Ceramic  (A1203  96%) 

9.6-10.4 

5. 5-7. 5 

25-30 

28.45 

38.5-40 

0.063 

0.135 

Copper 

227 

17-18 

32 

4.8-10 

17 

45-60 

0.323 

CuMn-Cu  <12.5*75- 12.5) 

62  8 

2. 8-3. 5 

64 

38 

19 

36 

0.3 

Cu- ln*Cu  (16-68- 16) 

81  9 

Cu-tn-Cu  (20-60-20) 

76-95  9 

5. 0-6. 3 

56 

32 

19 

36 

0.305 

Cu-Holy-Cu  (13-74-13) 

120  92 

5-6.5 

39 

0.357 

Epoxy 

0.13 

60 

8 

0.5 

0.066 

Epoxy  Glass 

.17-. 2 

12.4 

40 

2.5 

.065 -.092 

Epoxy  Kevlar 

0.13 

5.7-10 

3.5 

0.054 

fluorinated  Ethylene  Propy 

0.12 

26 

2 

0.078 

Gal l iun  Arsenide 

31 

5.7 

6.1 

12.3 

0.192 

German  iun 

38.5 

6.066 

18.8 

0.199 

Glass 

5.8 

8.6 

Glass  Fiber  (Si02) 

0.91 

.56-5 

60 

10.5 

0.079 

Gold 

173 

14.2 

18 

10.8 

30 

0.698 

Gold  Wire 

173 

14.2 

24 

10 

10 

1-12 

0.698 

Gold/Si l icon 

173 

9.8 

18 

10.95 

30 

0.698 

Graphite  Epoxy  Al 

90-100 

6.4-10.7 

Graphi te(75X)/Cu(25X) 

147 

5.6 

17.5 

0.131 

Indiun 

50 

32.1 

0.38 

1.57 

22 

0.263 

Kevlar  Fiber 

-4-(-2) 

Kovar 

14 

5. 0-5. 9 

75 

50 

20 

30 

0.302 

Holybdenun 

81 

5 

95 

80 

47 

10 

0.369 

Nickel 

53 

13.3 

50 

37 

0.322 

Parylene  Conf  Coat 

35-70 

8.2 

6 

Polyimide 

.006 

40 

0.4 

0.063 

Polyimide  Glass 

0.17-0.35 

11-14.2 

50 

2.5  -  3.0 

.066-. 090 

Polyimide  Kevlar 

.13-. 15 

3-7 

30 

3. 0-4.0 

.052-. 060 

Polyimide  Quartz 

.13-. 3 

5.0-10 

.070-. 061 

Polyurethane  Conf  Coat 

100-200 

Polyvinyl  Chloride 

.07-. 1 

15-18 

4.5 

0.057 

Quartz  Fiber 

0.54 

Sit  icon 

72.6 

2.3-4 

5 

23.7 

0.084 

Si  Iver 

243 

19.7 

18.13 

8 

10.3 

48 

0.379 

Solder 

25-30 

21.4-24.6 

2. 7-7. 8 

0.34 

4.49 

37 

0.3 

Steel  (low  C) 

27 

10-11 

30 

0.283 

Thermoplastic  Resin 

25-30 

7ungsten(75X)/Cu(25%) 

175 

38 

0.49 

*  Modulus  of  Rupture 
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TABLE  2-4.  GENERAL  MATERIAL  PROPERTIES  -  2 


C.ass 

Type 

Material 

Temp 

k 

rr 

CTE 

E 

G 

V 

P 

e 

Ftu 

Fsu 

Ft  • 

C 

W/in*C 

J/lD'C 

ppm/C 

TO  1 

TO  i 

pci 

X 

k$i 

ksi 

ksi 

Metal 

Pure 

Au 

25 

8.07 

58.05 

14.20 

11.31 

3.98 

0.42 

0.70 

30.00 

14.94 

10.01 

ni  l 

Metai 

Com* 

A- 

2J 

. 

5V.41 

»*.  20 

lU.ttu 

u./  U 

jO.UU 

16. uO 

«2.uc 

r  l 

Metal 

Wire 

Au 

20 

7.62 

59.41 

14.20 

10.00 

3.52 

0.42 

0.70 

1*12 

24.00 

16.08 

10.00 

Metal 

Wire/hard 

Au 

20 

7.62 

59.41 

14.20 

10.80 

3.80 

0.42 

0.70 

.5-2.5 

62.00 

41.54 

50.00 

Metal 

Alloy 

Au-si 

20 

7.62 

59.41 

9.80 

10.95 

4.21 

0.30 

0.70 

30.00 

18.00 

12.06 

ni  l 

Metal 

Alloy 

Au70-Pt30 

20 

8.07 

58.05 

14.20 

16.51 

5.81 

0.42 

0.72 

50.00 

92.70 

62.11 

Metal 

Pure 

Afl 

20 

10.87 

106.58 

19.00 

10.30 

3.76 

0.37 

0.36 

48.00 

18.13 

12.14 

8.00 

Metal 

Coeaa 

Ag 

0 

106.34 

106.12 

20.61 

10.30 

3.75 

0.38 

0.38 

48.00 

22.00 

14.74 

8.00 

Metal 

Pure 

Al  or  At  1199 

25 

6.27 

408.16 

23.60 

8.99 

3.63 

0.24 

0.10 

50.00 

50.00 

35.00 

20.00 

Metal 

Comm 

Al  1060 

20 

5.94 

408.16 

23.60 

10.00 

3.76 

0.33 

0.10 

43.00 

10.00 

7.00 

4.00 

Metal 

Alloy 

Al  IXSi  or  Mg 

20 

5.94 

408. 16 

23.60 

9.00 

3.76 

0.33 

0.10 

0.5-5 

34-62 

23-42 

27-50 

Metal 

Pure 

Cd 

20 

2.49 

104.31 

31.30 

7.98 

2.78 

0.33 

0.31 

50.00 

10.30 

6.82 

Metal 

Elec 

Cr 

20 

1.70 

2C8.33 

6.20 

0.04 

2.70 

0.33 

0.26 

0.00 

ic.04 

8.06 

Metal 

Pure 

Cu 

20 

10.11 

175.06 

16.50 

18.13 

6.73 

0.34 

0.32 

60.00 

30.30 

20.30 

4.83 

Metal 

AS-plated 

Cu/CuSOA  bath 

20 

10.11 

175.06 

9.50 

16.00 

5.96 

0.34 

0.32 

15-35 

25-35 

17-23 

25.00 

Metal 

Plate 

Cu/CuSOt  bath 

20 

10.11 

175.06 

9.50 

10.00 

3.72 

0.34 

0.32 

15-35 

25-35 

17-23 

25.00 

Metal 

AS-plated 

Cu/f luroborate 

20 

10.11 

175.06 

9.28 

12.00 

4.47 

0.34 

0.32 

15-32 

28-38 

19-25 

28.00 

Metal 

Plate 

Cu/f lurotoorata 

20 

10.11 

175.06 

9.28 

6.00 

2.23 

0.34 

0.32 

15-32 

28-38 

19-25 

28.00 

Metal 

Pure 

Pd 

20 

1.93 

111.11 

11.76 

16.30 

6.13 

0.33 

0.43 

30.00 

25.00 

16.75 

5.0u 

Metal 

Pure 

In 

20 

2.20 

105.67 

32.10 

f  .57 

0.59 

0.33 

0.26 

22.00 

0.38 

0.25 

Metal 

Coma 

Ti 

20 

0.52 

237.30 

8.82 

15.50 

6.50 

0.19 

0.16 

24.00 

35.00 

18.38 

25.00 

Metal 

Alloy 

KOVAR 

25 

0.44 

233.50 

5.86 

20.00 

7.59 

0.32 

0.30 

30.00 

75.00 

50.25 

50.00 

Metal 

Coma 

Pt 

20 

1.81 

59.86 

9.10 

24.80 

8.92 

0.39 

0.78 

30.00 

18.00 

12.06 

2.00 

Metal 

Pure 

Mi 

100 

2.11 

213.61 

13.30 

30.02 

11.02 

0.31 

0.32 

30.00 

45.97 

30.80 

8.60 

Metal 

AS-plated 

Mi-diff  baths 

20 

2.13*2.78 

213.61 

13.60 

21-31 

8-11.8 

0.31 

0.32 

5-35 

50-152 

33*102 

33.00 

Metal 

Alloy 

Mi7C-Cr30 

20 

0.36 

208.83 

12.20 

24.00 

9.00 

0.33 

0.29 

30.00 

128.00 

85.76 

Metal 

Al  toy 

Mi50-Fe50 

20 

218.36 

8.46 

24.00 

9.00 

0.33 

0.30 

0.30 

72.00 

48.24 

Metal 

Coma 

Sn-betta 

25 

1.54 

9.98 

23.10 

6.40 

2.41 

0.33 

0.26 

53.00 

1.30 

0.87 

Metal 

Pure 

MO 

20 

3.61 

125.17 

5.22 

47.00 

17.80 

0.32 

0.37 

10.00 

95.00 

63.65 

80.00 

Metal 

Pure 

P*3 

20 

0.89 

58.50 

29.30 

2.00 

0.69 

0.44 

0.41 

47.00 

2.00 

1.82 

Metal 

Alloy 

Sn63-PbJ7 

20 

1.29 

165.50 

21.40 

4.49 

1.60 

0.40 

0.30 

37.00 

7.80 

5.40 

0.34 

Metal 

Alloy 

Sn5-Pb93 

20 

0.89 

60.70 

28.40 

2.68 

0.94 

0.42 

0.39 

45.00 

4.00 

2.10 

0.20 

Ceramic 

96X 

A1203 

20 

7.87 

399.09 

7.10 

*0.00 

16.39 

0.22 

0.14 

0.06 

25.00 

16.76 

28.45 

Ceramic 

99X 

A 1203 

20 

9.84 

399.09 

6.50 

40.00 

16.39 

0.22 

0.14 

0.08 

30.00 

20.10 

28.45 

Ceramic 

Sapphire 

A 1203 

100 

0.64 

341.55 

8.00 

50.00 

21.50 

0.14 

* 

65.00 

Ceramic 

Diamond 

C 

20 

50.80 

0.80 

114.00 

21.00 

0.20 

Ceramic 

Coma 

AIN 

20 

55.12 

317.46 

4.10 

40.00 

16.00 

0.25 

0.12 

0.07 

28.00 

18.76 

56.90 

Ceramic 

Coma 

SiC 

20 

27.56 

362.81 

3.80 

59.00 

25.65 

0.15 

0.12 

0.04 

26.00 

17.42 

34.02 

Ceramic 

Comm 

BcO 

20 

102.36 

453.51 

7.20 

38.00 

14.18 

0.3*- 

0.10 

0.04 

13.50 

9.05 

24.18 

Ceramic 

HP  Grade 

BN 

20 

23.62 

294.78 

0.00 

6.20 

2.52 

0.23 

0.07 

O.U 

6.50 

4.36 

7.59 

Ceramic 

Fused  Silica 

Si02 

20 

0.04 

357.30 

0.56 

10.50 

4.53 

0.16 

0.08 

Ceramic 

Com  7052 

Boros i licate 

20 

0.03 

436.43 

5.00 

8.20 

3.36 

0.22 

0.08 

0.12 

10.00 

6.70 

6.00 

Ceramic 

Corn  9606 

Pryoceram 

20 

0.09 

442.33 

5.76 

17.30 

7.09 

0.22 

0.09 

0.08 

13.00 

8.71 

20.00 

Ceramic 

CERMET 

70Al2o3-30Cr 

25 

8.65 

52.30 

33.00 

0.21 

0.17 

0.07 

35.00 

23.45 

55.00 

Ceramic 

CERMET 

77Cr-23Al203 

25 

1.28 

303.40 

8.93 

37.50 

17.00 

0.21 

0.21 

0.06 

21.00 

40.00 

45.00 

Ceramic 

Forster i ta 

2Mg0-Si02 

20 

0.03 

9.80 

21.00 

8.54 

0.23 

0.10 

23.00 

Semiconductor 

Coma 

GaAs 

20 

1.37 

10.36 

5.70 

12.30 

4.77 

0.29 

0.19 

6.10 

Semiconductor 

Coma 

Si 

20 

3.19 

341.55 

2.33 

23.70 

9.63 

0.23 

0.08 

5.00 

Semiconductor 

Coma 

Ge 

20 

1.69 

140.48 

6.07 

18.80 

7.34 

0.28 

0.20 

Polymer 

Kapton  H 

Polyimide  film 

23 

0.00 

494.33 

20.00 

0.43 

0.16 

0.34 

0.05 

75.00 

25.00 

10.00 

Polymer 

Coma 

Polyimide 

23 

0.01 

494.33 

40.00 

0.40 

0.17 

0.33 

0.06 

Polymer 

Coma 

Epoxy 

23 

0.01 

800.00 

60.00 

0.50 

0.19 

0.35 

0.07 

8.00 

2.00 

Polymer 

Cast  Rigid 

Epoxy 

23 

0.01 

854.29 

59.40 

0.45 

0.17 

0.35 

0.04 

4.40 

9.50 

6.37 

Polymer 

Epoxy- conduct 

Ablebond  BA 

23 

0.05 

333.00 

55.00 

0.82 

0.31 

0.33 

0.09 

1.60 

Where: 

k  =  thermal  conductance  P  =  density 

Cp  =  Specific  Heat  e  =  elongation 

CTE  =  coefficient  of  thermal  expansion  Ftu  =  tensile  ultimate  strength 

E  =  modulus  of  elasticity  Fsu  =  shear  ultimate  strength 

G  =  modulus  of  rigidity  Fty  =  tensile  yield  strength 

V  =  poisson's  ratio 
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TABLE  2-5.  CTE  FOR  ELECTRONIC  MATERIALS 


Material 

CTE 

Material 

CTE 

ppm/C 

ppm/C 

•  1  /.  ^ 

5. * 

“■'lytdenur 

4  n 

Alunina  94X 

6.4 

Monet 

14 

Alenina  96% 

6.4 

Multi te 

2.3 

Alenina  99. 5% 

6.4 

Nickel 

13-15 

Aluminum  &  Alloys 

22-28 

Nickel  Silver 

16.2 

BeO  99. SX 

6.4 

Phosphor  Bronte 

17.8 

Berytifua  Copper 

17.8 

Platinuo 

8.9 

Brass  (66Cu,34Z.,; 

20.3 

Polycarbonates 

50-70 

Cadniun 

29.8 

Polyimides 

40-50 

Constanten  (45Ni,  55Cu) 

14.9 

Polyurethanes 

180-250 

Copper  &  Alloys 

16-18 

Porcelain  on  Steel 

11.4 

Oynamould  1Q3B  (epoxy) 

22 

RTV 

800 

Oynamould  1048  (epoxy) 

25 

Sapphire 

4.3 

Epoxies 

60-80 

Si l icon 

2.3 

Glass 

8.6 

Silicon  Nitride 

2.3 

Gold 

14.2 

Silver 

19.7 

Gold-Tin  Eutectic 

16 

Steel,  Low  Carbon 

12 

Invar 

2 

Steel,  SAE  1045 

15 

Iron 

11.7 

Sytgard 

300 

Kevlar 

-2 

Tin 

23 

Kovar 

6.2 

Ti  taniua 

10 

lead 

29 

Tungsten 

4.5 

Lead-Tin  Eutectic 

21 

Zinc 

39.7 

Magnesium 

25.2 

Material 


Solder 
Copper 
Qold 
Epoxy  Glass 
Polyimide  Glass 
Steel 
Aluminum 
Beryllium  Oxide 
Epoxy  Kevlar 
Alumina 
Polyimide  Quartz 
Cu-ln-Cu  20-60-20 
Cu  Clad  Molybdenum 
Kovar 
Glass  Fiber 
Alloy  42 
Polyimide  Kevlar 
Cu-ln-Cu  12-76-12 
Silicon 
Quartz 

Carbon  Fiber/Epoxy 
Kevlar 


-5  0  5  10  15  20  25  30 

Coefficient  of  Thermal  Expansion  (CTE) 


Figure  2-2.  CTE  for  Common  Materials 
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TABLE  2-6.  PCB  LAMINATE  PROPERTIES 


Printed  Circuit  Board  Laminates 


Material 

Clast  Transition 

CTE  x,y 

CTE  Z 

Thermal  Conductivity 

Modulus, El  attic 

Density 

Tg 

(C) 

ppm/“C 

ppm/-C 

(W/m*K> 

(10E6  psi) 

(g/ear-3) 

E-glass/epo*v 

125 

12-16 

80-90 

0.35 

2.5 

1.8 

E-gless/polyimide 

200-260 

11-14 

60 

0.35 

2.8 

1.83 

E-glass/PTFE 

75 

24 

261 

0.26 

0.14 

Kevlar  120U/0uartex 

185 

3-8 

105 

0.16 

3. 2-4.0 

1.5 

Kevlar  120U/potyimide 

180-200 

3-8 

33 

0.12 

4 

1.44 

Kevlar  120T/0uartex 

185 

0.75 

6.5 

Kevlar  108U/4093  resin 

125-135 

4-7 

50-110 

2. 3-2. 8 

Kevlar  108T/0uartex 

185 

.75-1.25 

66-82 

5.5 

Kevlar  2643MU/0uartex 

185 

5-8 

50-110 

2. 0-2. 3 

Ouartz/Polyimide 

260 

6-t2 

34 

0.13 

4 

1.68 

Ouartz/Ouartex 

185 

62 

2.7 

RO  2800 

16-19 

24-30 

0.44 

0.6 

6061  Aluairun 

23.6 

23.6 

190 

10 

2.71 

Copper 

17.3 

17.3 

131 

17 

8.95 

Cu-ln-Cu  20-60-20 

lateral  6.02 

6.02 

164 

19 

8.45 

normal  6.02 

6.02 

22 

19 

8.45 

Cu-In-Cu  16-68-16 

lateral  5.22 

5.22 

132 

IV 

8.59 

normal  5.22 

5.22 

19.6 

19 

8.39 

C,.  !"-Cu  12.5-75-12.5 

laterat  3.69 

3.69 

110 

19 

8.31 

normal  3.69 

3.69 

17.8 

19 

8.31 

Steel,  low  carbon 

10 

10 

27 

30 

7.84 

Alenina  94  X 

6.4 

6.4 

21 

37 

3.6 

Cu-Mo-Cu  13-74-13 

lateral  6.5 

6.5 

208 

39 

9.89 

normal  6.5 

6.5 

159 

J9 

9.89 

BeO  99. 5X 

6.4 

6.4 

208 

39 

’.77 

Alininua,  Boron  reinforced 

(20-50*1 

4.5- t2. 7  4 

.5-12.7 

120- 188 

12-30 

2.0-2.63 

Alunimm,  graphite  reinforced  (40-60X) 

-.32-3.6  - 

.32-3.6 

310-419 

2.41 

2.3  Solder  Properties 

This  section  contains  information  on  solder.  Information  includes  strength  properties 
at  various  temperatures  and  cyclic  frequencies,  cycles  to  failure  based  on  chip  lead  count, 
stress  vs.  strain  curves,  cycles  to  failure  based  on  stress  or  strain,  stress  relaxation  proper¬ 
ties,  creep  properties  and  modulus  of  elasticity  variations  with  temperature  and  cycle  rate. 


Table  2-7  (Refs.  2-2, 3, 4, 5)  is  a  summary  of  solder  mechanical  and  thermal  properties. 
Information  includes  tensile  and  shear  strength  at  various  temperatures,  creep- stresses  to 
obtain  a  1000  hour  life  and  thermal  conductivity. 
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TABLE  2-7.  SOLDER  PROPERTIES 


filename  a: CAARFEMXSOLDEft 
contains  solder  properties 

Creep  Stress  Thermal 


Alloy 

Tensile  Strength  (PSD 

Shear  Strength 

(PSD 

1000  hour 

life* 

Conductivity 

20  C 

100  C 

*130  C 

25  C 

150  C 

20  C 

100  c 

U/m/K 

63Sn/37Pb 

6120 

2700 

12700 

4130 

1165 

50.9 

6QSn/4QPb 

2700 

580 

4800 

4206 

653 

49.8 

50Sn/50Pb 

5945 

11100 

3515 

1100 

47.8 

96.5Sn/3.5Ag 

5260 

16600 

4650 

1510 

99Sn/1Sb 

15100 

2900 

1125 

95Sn/5Sb 

4410 

2900 

18150 

4625 

1880 

3046 

1305 

10Sn/90Sb 

2850 

1160 

7300 

2800 

1500 

97.5Pb/2.5Ag 

5300 

2590 

1440 

97.5Pb/1.5Ag/1Sn 

4980 

5900 

3040 

1520 

90Pb/5In/5Ag 

6220 

3470 

1755 

307 

*  Initial  creep  stress  (lb/in2)  for  life  of  1000  hours 

Sources:  1)  Surface  Mount  Technology;  Barbara  Roos*Ko*el;  International 

Society  for  Hybrid  Microelectronics.  1984. 

2)  Surface  Mount  Technology;  Peter  H.  Moy;  International  Society 
for  Hybrid  Microelectronics,  1984. 

3)  Soldering  in  Defense  Electronics;  K.  Nagesh;  Bharat 
Electronics  Limited 

Figure  2-3  (Ref.  2-6)  illustrates  cycles  to  failure  data  for  leadless  chip  carrier  solder 
joints  based  upon  the  number  of  leads  per  package.  This  data  was  based  on  a  one  hour 
thermal  cycle  of  -55C  to  125C.  Two  different  substrates  were  tested.  The  first  PCB  was 
polyimide  glass  with  an  aluminum  heat  sink,  resulting  in  a  combined  CTE  of  15  ppm/C. 
The  second  was  made  of  polyimide  glass  with  a  CTE  of  14  ppm/C.  The  components  had 
ceramic  packages.  Figure  2-4  (Ref.  2-7)  illustrates  a  method  of  predicting  cycles  to  fail¬ 
ure  based  upon  the  strain  factor  of  the  chip-PCB  combination.  The  strain  factor  consid¬ 
ers  factors  such  as  CTE,  temperature  range,  chip  size  and  solder  height  and  is 
represented  by  the  following  equation: 

SF  =  L  (  Q^-ai  )  AT  { 

2  J  2*  h 


where:  0*2  =  PCB  CTE 


c*i  =  component  CTE 
L  =  component  diagonal  length 
h  =  solder  joint  height 
AT  =  temperature  cycle  range 
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Figure  2-3.  Thermal  Fatigue  Data  for  LCC  Solder  Joints  by  I/O  Count 
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Figure  2-5  (Ref.  2-8)  illustrates  the  relationship  between  the  principal  stress  amplitude 
and  the  strain  amplitude  at  room  temperature.  Figure  2-6  (Ref.  2-9)  is  similar  except  it 
plots  shear  stress  as  a  function  of  shear  strain. 


Stress 

Amplitude 

ksi 


Strain  Amplitude  -  in./in. 

G  P03-0624- 1 5'D' Vs 


Figure  2-5.  Stress  vs.  Strain  for  Solder 


Figure  2-6.  Shear  Stress-Strain  Curve  for  Solder  63-37 
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The  next  group  of  figures  will  detail  cycles  to  failure  information  for  solder.  Figure 
2-7  (Ref.  2-8)  plots  cycles  to  failure  as  a  function  of  principal  stress  caused  by  bending 
for  solder  which  is  50%  tin  and  50%  lead.  The  figure  contains  separate  slopes  for  ran¬ 
dom  and  constant  amplitude  stress  cycles.  Figure  2-8  (Ref.  2-8)  is  identical  to  Figure  2-7 
except  the  solder  is  63%  tin  and  37%  lead.  Figures  2-9  and  2-10  (Ref.  2-8)  are  identical 
to  Figures  2-7  and  2-8,  respectively,  except  the  stress  is  a  fully  reversed  shear  stress  in¬ 
stead  of  a  principal  stress.  Figures  2-11  (Refs.  2-9,10)  and  2-12  (Ref.  2-8)  also  plot 
cycles  to  fai'ure  as  a  function  of  shear  stress.  Figure  2-12  shows  the  effect  cycle  rate  and 
temperature  has  on  the  fatigue  life  of  solder.  By  comparing  Figures  2-10,  2-11  and 
2-12,  variations  in  recorded  material  properties  becomes  obvious,  illustrating  the  impor¬ 
tance  of  carefully  choosing  data  as  an  input  to  finite  element  analyses.  The  life  variations 
as  a  function  of  temperature  and  stress  cycle  frequency  mandates  that  the  operational  en¬ 
vironment  be  correctly  quantified.  Figure  2-13  (Ref.  2-6)  plots  cycles  to  failure  as  func¬ 
tion  of  shear  strain  range  as  opposed  to  shear  stress.  This  data  was  collected  for 
polyimide  glass  PCBs  and  polyimide  glass  PCBs  with  an  aluminum  heat  sink.  Figures 
2-14  (Refs.  2-9,10)  and  2-15  (Ref.  2-7)  also  plot  cycles  to  failure  as  a  function  of  shear 
strain  range.  Figure  2-15  also  plots  strain  range  at  various  temperatures  and  cycles  rates. 


Figure  2-7.  S-N  Fatigue  Curve  for  Soft  Solder  (50%  Lead  -  50%  Tin),  Reversed  Bending 
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Figure  2-8.  S-N  Fatigue  Curve  for  Soft  Solder  (37%  Lead  -  63%  Tin),  Reversed  Bending 
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Figure  2-9.  S-N  Fatigue  Curve  for  Soft  Solder  (50%  Lead  -  50%  Tin),  Reversed  Shear 
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Figure  2-10.  S-N  Fatigue  Curve  for  Soft  Solder  (37%  Lead  -  63%  Tin),  Reversed  Shear 


Figure  2-11.  S-N  Fatigue  Curve  for  67-37  Solder  at  Room  Temperature 
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Figure  2-1 2.  Stress  vs.  Cycles  to  Failure 
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Figure  2-13.  Thermal  Fatigue  Data  for  LCC  Solder  Joints 
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F'-ure  2-14.  Strain-Life  Curve  for  67-37  Solder  at  Room  Temperature 


Figure  2-15.  Strain  Range-Fatigue  Life  Plots  for  63/37  Sn/Pb  Solder  at  Various  Conditions 

The  following  set  of  figures  will  illustrate  a  solder  phenomenon  which  raises  the  level 
of  difficulty  in  making  accurate  life  predictions  -  stress  relaxation  and  creep.  Stress  relax¬ 
ation  occurs  when  a  given  strain  is  imposed  on  solder.  If  this  strain  is  maintained,  the 
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stress  in  the  solder  will  begin  to  relax  or  disappear.  This  continues  until  the  solder  is 
virtually  stress  free.  Creep  occurs  when  a  constant  stress  is  applied  to  solder.  The  solder 
will  begin  to  elongate,  or  stretch,  while  the  stress  is  maintained.  When  the  stress  is  re¬ 
laxed,  the  elongation  does  not  return  to  zero.  Permanent  deformation  has  occurred. 
Figure  2-16  (Ref.  2-11)  plots  remaining  stress  versus  time  for  63/37  solder  at  various 
temperatures.  As  shown,  the  higher  the  temperature,  the  quicker  the  stress  relaxes. 

Figure  2-17  (Ref.  2-10),  plots  the  time  to  fail  for  solder  subjected  to  a  constant  stress. 
This  is  plotted  for  several  different  temperatures.  Figure  2-18  (Ref.  2-10)  illustrates  the 
relationship  between  the  speed  of  the  applied  stress  and  the  stress  to  fail.  The  faster  the 
applied  load,  the  stronger  the  solder  is.  Figure  2-19  (Ref.  2-7)  illustrates  the  correlation 
between  the  life  time  of  solder  and  the  percentage  of  creep  per  cycle. 
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Figure  2-17.  Solder  Creep  Properties 
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Figure  2-18.  Strength  vs.  Strain  Rate 
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PERCENT  CREEPICYCLE 


Figure  2-19.  Lifetime  Data  for  Eutectic  Solder  Versus  Percent  Creep  Per  Cycle 


The  last  bit  of  data  for  solder  is  shown  in  Figure  2-20  (Ref.  2-7).  It  illustrates  the  re¬ 
lationship  between  the  modulus  of  elasticity  of  solder  and  the  temperature  and  stress 
cycle  frequency.  It  clearly  shows  that  solder  is  less  elastic  at  the  lower  end  of  the  normal 
thermal  range  for  military  hardware.  Additionally,  it  is  less  elastic  at  higher  stress  cycle 
frequencies. 
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Figure  2-20.  Solder  Modulus  With  Temperature 


2.4  Component  Lead  S-N  Curves 

The  Luliov  i''£  group  of  figures  (Figures  2-21,22,23,24)  provides  stress  cycle  to  failure 
data  for  common  lead  materials.  All  stresses  are  principal  stresses.  This  data  was  col¬ 
lected  from  Reference  2-8. 
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Figure  2-21.  Typical  Random  and  Sinusoidal  S-N  Fatigue  Curves  for  Kovar,  Reverse  Bending 


Figure  2-22.  S-N  Fatigue  Curve  for  Electrical  Lead  Wire  (99.9%  Copper  Cold  Drawn), 

Reversed  Bending 
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2-23.  S-N  Fatigue  Curve  for  Electrical  Lead  Wire  (99.9%  Nickel),  Reversed  Bending 

O 


Figure  2-24.  S-N  Fatigue  Curve  for  Electrical  Lead  Wire  (99.9%  Silver  Hard  Rolled), 

Reversed  Bending 
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Chapter  3 

FAILURE  MECHANISMS 


3.0  Introduction 

This  chapter  describes  failure  mechanisms  caused  by  mechanical  stresses  which  can  be 
analyzed  with  finite  element  stress  analysis  methods.  These  stresses  are  a  result  of  the 
deformations  caused  by  temperature  changes,  vibration,  shock,  and  high  G  aircraft 
maneuvering  loads.  These  deformations  result  in  material  fatigue,  creep,  fracture,  buck¬ 
ling,  and  eventual  component  failure.  The  following  paragraphs  discuss  these  failure 
mechanisms  in  more  detail. 

3.1  Deformation 

Temperature  changes  occur  during  power-on  cycles,  changes  in  altitude,  coolant  tem¬ 
perature  fluctuations,  and  diurnal  cycles.  Thermal  cycles  cause  materials  to  expand  and 
contract  proportionately  to  their  coefficient  of  thermal  expansion  (CTE).  Since  printed 
circuit  boards  and  the  components  mounted  on  them  have  different  CTEs,  they  deform 
by  different  amounts  as  the  temperature  changes  (Figure  3-1).  This  results  in  stresses  in 
the  solder  joints  and  lead  wires  connecting  the  components  to  the  boards. 

Vibration  is  caused  by  ground  transportation,  engine  operation,  noise,  aerodynamic 
buffet,  or  gun  fire.  Vibration  results  in  out-of-plane  deformation  of  printed  circuit 
ooards  and  high  frequency  stress  cycles  in  the  leads  and  solder  joints  of  the  components. 
As  the  board  deforms,  the  leads  and  solder  joints  are  pulled  and  compressed  as 
illustrated  in  Figure  3-1. 

Other  sources  of  deformation  include  shock  and  high  G  aircraft  maneuvers.  Shock  is 
caused  by  installation  of  line  replaceable  units,  accidental  dropping  of  equipment, 
detonations  during  battle  situations,  hard  landings  and  collisions.  The  resulting  rapid 
deformations  are  similar  to  those  of  transient  vibration.  High  G  aircraft  maneuvers 
i..‘.uOse  steady  loads  (lasting  up  to  several  seconds)  on  circuit  boards  during  the  duration 
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of  the  maneuver.  The  out-of-plane  deformations  caused  by  maneuvering  loads  are 
similar  to  the  vibration  displacements,  but  at  much  lower  frequencies. 


Figure  3-1.  Circuit  Board  Deformation  During  Temperature  Cycles  and  Vibration 

3.2  Fatigue 

Material  fatigue  is  characterized  by  the  initiation  of  cracks  at  areas  of  high  stress  fol¬ 
lowing  a  period  of  repeated  loads.  These  cracks  grow  under  subsequent  load  cycles  until 
fracture  of  the  material  occurs.  The  repeated  loads  are  caused  by  vibration,  temperature 
cycles  and  the  other  deformation  sources  described  in  Section  3.1.  Material  fatigue  data 
is  typically  presented  as  plots  of  cyclic  stress  amplitude  versus  cycles  to  failure,  commonly 
referred  to  as  S-N  curves  (Figure  3-2).  As  the  figure  suggests,  higher  stresses  result  in 
shorter  fatigue  lives.  When  the  stresses  exceed  the  yield  strength  of  the  material,  failure 
can  occur  in  relatively  few  cycles.  Fatigue  cracks  are  frequently  observed  in  the  solder 
joints  of  leadless  chip  carriers  (Figure  3-3)  and  J-leaded  surface  mount  components 
(Figure  3-4). 
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Figure  3-2.  Typical  Format  of  Fatigue  Data 


Figure  3-3.  Fatigue  Crack  In  a  Leadless  Solder  Joint 


Figure  3-4.  Fatigue  Crack  in  a  J-Lead  Solder  Joint 
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3.3  Creep  and  Stress  Relaxation 


Under  a  constant  force,  materials  such  as  solder  will  deform  with  time,  or  creep.  Under 
a  constant  displacement,  the  stress  in  solder  relaxes  with  time.  These  conditions  are  illus¬ 
trated  in  Figure  3-5.  Models  of  combined  creep  and  stress  relaxation  have  been  developed 
by  Hall  (Reference  3-1)  and  creep  data  for  solder  is  documented  in  Reference  3-2. 

Material  Creeps  at  Material  Relaxes  Under 

a  Constant  Stress  a  Constant  Displacement 


Figure  3-5.  Solder  Creep  and  Stress  Relaxation  Behavior 


Stress  relazation  can  occur  under  the  following  circumstances.  When  electrical  power 
is  turned  on,  the  components  and  PCBs  heat  up  and  begin  to  expand.  With  time,  the 
temperature  will  reach  a  maximum.  At  this  point,  the  stress  induced  by  the  CTE  mis¬ 
match  will  be  a  maximum,  as  will  the  strain.  If  power  remains  on,  the  strain  will  remain 
constant,  but  the  stress  will  begin  to  relax.  This  relaxation  can  take  place  quite  rapidly 
(Reference  3-3).  Figure  3-6  (Reference  3-4)  illustrates  stress  relaxation  versus  time  for 
solder  at  various  temperatures.  This  Figure  shows  that  between  75-95%  of  stress  relaxes 
within  5  minutes  of  strain  stabilization.  This  phenomenon  will  be  used  to  greatly  simplify 
the  compilation  of  load  histories  (discussed  in  Chapter  5). 


From  a  finite  element  modellers  point  of  view,  failure  is  accelerated  for  the  following 
reasons.  As  mentioned  earlier,  when  a  material  which  exhibits  stress  relaxation  proper¬ 
ties  is  held  at  a  constant  displacement,  the  stress  will  dissipate  to  zero  with  time.  When 
the  displacement  returns  to  the  original  position,  the  material  will  be  compressively 
loaded  at  a  value  equal  in  magnitude  to  the  tensile  loading  it  was  under  prior  to  stress 
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Data:  L.  Fox  -  DEC, 

J.B.  Enns  and  C.J.  Aloisio, 
AT&T  Bell  Laboratories 
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psi 


Figure  3-6.  Solder  Joint  Stress  Relaxation  63/37  Sn/Pb 


relaxation.  Therefore,  the  part  has  been  subjected  to  a  stress  amplitude  twice  as  large  as 
a  material  which  does  not  exhibit  stress  relaxation  (Figure  3-7). 

For  loading  conditions  which  cause  creep  to  occur,  the  failure  mechanism  which  can  be 
modelled  is  the  strain  amplitude.  The  strain  amplitude,  in  this  case,  includes  the  initial 
strain  (before  creep  begins)  and  the  additional  strain  which  results  from  creep.  There¬ 
fore,  the  total  strain  is  much  greater  than  for  materials  which  do  not  creep  (Figure  3-8). 


Creep  can  occur  when  a  constant  load  source  is  confined  by  a  material  such  as  solder. 
There  are  few  constant  load  sources  in  avionics:  however,  one  situation  closely  approximates 
this  condition.  During  manual  soldering  operations,  the  leads  of  gull  wing  surface  mount 
devices  may  be  held  down  flush  on  the  pad  while  the  solder  is  flowing.  Once  the  solder  so¬ 
lidifies  and  the  lead  is  released,  a  constant  force  will  be  exerted  on  the  solder  by  the  lead  as 
it  attempts  to  return  to  its  original  position.  This  could  cause  the  solder  to  creep. 
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Figure  3-7.  Hysteresis  Curve  for  Materials  Exhibiting  Stress  Relaxation 
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Figure  3-8.  Hysteresis  Curve  for  Materials  Which  Exhibit  Creep 

3.4  Fracture 


3.4.1  Ductile  Fracture 


Pre-existing  flaws  and  cracks,  which  initiate  under  repetitive  fatigue  cycles,  propagate 
through  the  material  until  they  reach  a  critical  size.  Additional  cycles  result  in  fracture.  Sus¬ 
tained  stresses  in  solder,  for  example,  can  result  in  creep  deformation  and  eventual  rupture. 
Creep  can  also  act  in  combination  with  cyclic  fatigue  stresses  to  cause  fracture  after  a  given 
number  of  stress  cycles.  A  reliability  analysis  under  these  conditions  requires  a  fatigue  cal¬ 
culation  to  determine  the  number  of  cycles  to  failure,  or  a  creep  analysis  to  compute  the 
time  to  failure. 
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Fracture  can  also  occur  immediately  if  the  applied  stress  exceeds  the  ultimate  strength 
of  the  material.  If  high  stresses  are  limited  to  a  small  section  of  a  ductile  material  such  as 
solder,  the  material  *vdi  yield  in  that  region  and  stresses  will  be  redistributed  through  the 
rest  of  the  material  (Figure  3-9).  As  the  load  is  increased,  the  plastic  zone  will  increase 


Figure  3-9.  Stress  Behavior  In  a  Ductile  Material  as  Applied  Stress  Increases 
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in  size  until  the  entire  cross  section  has  reached  the  yield  strength  of  the  ductile  metal.  If 
the  load  continues  to  increase,  and  the  stress  reaches  the  ultimate  strength  of  the 
material,  rupture  occurs. 

3.4.2  Brittle  Fracture 

Brittle  materials  such  as  ceramics  generally  do  not  yield  or  undergo  plastic  deforma¬ 
tions.  If  stresses  exceed  the  ultimate  strength  of  the  material  in  any  section  of  the  part, 
the  material  will  fracture.  A  reliability  analysis  must,  therefore,  include  a  check  of  the 
stress  levels  in  the  part  to  determine  if  stresses  have  exceeded  the  yield  strength  in  ductile 
materials  or  the  ultimate  strength  in  brittle  materials. 

3.5  Buckling 

If  a  short  column  is  loaded  in  compression,  it  will  remain  straight  and  the  stresses  and 
strains  computed  in  a  static  finite  element  analysis  are  valid.  However,  a  different  type  of 
behavior  occurs  for  long,  slender  columns.  When  the  compressive  load  reaches  a  critical 
level,  the  slender  column  will  bow  out  of  the  plane  of  its  axis,  or  buckle.  The  column  is 
unable  to  carry  any  further  loads  because  large  lateral  deflection  occurs  with  little 
increase  in  load.  Buckling  is  illustrated  in  Figure  3-10. 


3.6  Failure  Modes  of  Electronics  due  to  Failure  Mechanisms 

Table  3-1  summarizes  the  typical  failure  modes  of  electronic  hardware.  Additionally, 
it  correlates  the  failure  mechanisms,  types  of  loading  and  probable  environment  with  the 
individual  failure  modes. 
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TABLE  3-1.  PREDOMINANT  FAILURE  MODES 
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Chapter  4 

FINITE  ELEMENT  TECHNIQUES 


4.0  Introduction 

Chapter  4  will  identify  concepts  which  increase  the  accuracy  of  FEAs  including  tech¬ 
niques  addressing  mesh  generation,  model  minimization  and  choice  of  analyses  (linear  vs. 
nonlinear).  Example  problems  will  demonstrate  FEA  accuracy  when  dealing  with  very 
small  structures  and  structures  made  of  materials  which  have  large  variations  in  their 
elastic  modulus.  Additionally,  example  problems  will  compare  FEAs  of  solder  joints  (for 
J-lead  and  leadless  devices)  using  various  FE  codes.  An  example  problem  will  also  be 
used  to  demonstrate  the  transfer  of  PCB  deflections  into  loading  on  component  leads. 
Before  the  example  problems  are  discussed,  a  brief  tutorial  on  finite  element  analyses  will 
be  given.  Discussions  will  cover  finite  element  methods,  finite  difference  methods  and 
hand  calculations. 

4.1  Numerical  Methods  in  Engineering  Analysis 

Engineers  in  the  electronics  industry  use  a  wide  variety  of  modeling  techniques  to 
investigate  the  thermal  and  structural  properties  of  electronics  systems  operating  in  many 
different  environments.  Mechanical  phenomena  typically  studied  by  engineers  can  be 
described  by  the  laws  of  physics  in  terms  of  algebraic,  differential,  or  integral  equations 
relating  various  quantities  of  interest.  While  the  derivation  of  the  governing  equations 
for  most  problems  is  not  unduly  difficult,  their  solution  by  exact  methods  of  analysis  is  a 
formidable  task.  The  methods  used  will  often  depend  upon  the  complexity  of  the  prob¬ 
lem,  the  time,  the  manpower,  the  funds  and  knowledge  of  the  engineers.  The  two  basic- 
methods  of  modeling  and  analysis  are  generally  described  as  either  closed-form  or 
numerical  analysis  techniques.  Closed-form  solutions  are  an  easy  and  efficient  form  of 
"'hand"  calculations.  Numerical  methods,  i.e.,  finite  difference,  finite  element,  boundary 
element,  and  statistical  energy,  enable  engineers  to  analyze  structures  too  complex  for 
closed-form  solutions.  By  far  the  most  common  of  these  numerical  techniques  is  the 
finim  element  method  (FEM). 
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4.1.1  Finite  Element  Methods 

In  the  finite  element  method,  mechanical  systems  and  structures  are  represented  by  a 
discrete  grid  of  node  points  interconnected  by  various  types  of  structural  elements  form¬ 
ing  a  finite  element.  The  complete  solution  is  obtained  by  comoining  the  individual 
elements  into  an  idealized  structure  for  which  the  conditions  of  equilibrium  and 
compatibility  are  satisfied  at  the  nodes  of  the  elements.  In  the  finite  element  method,  the 
assumed  displacement  fields  within  a  finite  element  are  assumed,  by  the  use  of  energy 
theorems,  to  derive  a  stiffness  matrix  relating  the  nodal  forces  to  the  nodal  displacements 
of  the  element.  If  the  equilibrium  conditions  are  applied  at  each  node,  then  a  set  of 
simultaneous  equations  can  be  assembled  and  solved  for  all  the  displacements  in  the 
structure. 

As  an  example,  the  stiffness  matrix  for  a  uniform  bar.  Figure  4-1,  can  be  derived. 
Given  the  governing  equation  derived  from  the  equilibrium  of  forces  within  the  bar: 

i(EAf)-W  =  0 

where  f(x)  is  the  internal  force  per  unit  length  and  EA  is  the  elastic  rigidity  (E  is 
Young’s  modulus  and  A  is  the  cross  sectional  area  of  the  bar).  The  boundary  equations 
for  the  problem  in  Figure  4-1  are: 

du 

P  =  EA—  at  x  =  L 
dx 

u  =  0  at  x  =  0  Eq.  2 

From  the  math  model  in  Figure  4-2,  the  displacement  along  the  length  of  the  bar  can 
be  expressed  by  the  unit  displacement  theorem: 

ux  =  ut  +  (u2  — ui)^  Eq.  3 
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The  strain  in  the  bar  is  given  by: 


= 


du 

dx 


Substituting  in  the  displacement  relationship  from  Eq.  3  yields: 


<Tx  =  —  (u2-ul) 


Expressing  the  strain  in  vector  form: 

Therefore,  the  strain-displacement  relationship  is: 


ui 

U2 


ib]  =  l[-i,i] 

The  stiffness  matrix  is  equivalent  to: 

[K]  =  ^  (b]T  [EJ  [b]  dv 


Eq.  4 


Eq.  5 


Eq.  6 


Eq.  7 


Eq.  8 


where  the  [E]  matrix  is  E. 


For  the  problem  at  hand,  the  stiffness  matrix  is: 

■L 

L 


[K] 


=  [LI 

Jo  L 


-  1 
1 


^[-1,1]  A  dx 


After  integrating  along  the  length: 


[K]  = 


EA 


1  -1 

-1  1 


Eq.  9 


Eq.  10 
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For  the  thermal  stiffness,  [H]: 


[H]  =  J  [b]T  [Et]  dv  Eq.  11 


Similarly, 


Eq.  12 


The  complete  force-displacement  relationship  is: 


Pi 

EA 

1 

-  1 

Ul 

p2 

L 

-  1 

1 

u2 

+  EA  cn  AT 


Eq.  13 


where  a  is  the  coefficient  of  thermal  expansion  and  AT  is  the  temperature  range. 

To  complete  the  problem,  discretize  the  domain  (Figure  4-3)  and  assemble  the  global 
stiffness  matrix. 


f(x) 


'Ai 


l  wMmmmmz a-»p-$ 

►  X 

- L - *H 


G  PO 3-0624-1  -D/ks 

Figure  4-1.  Uniform  Bar 


GP03-0624-2-D/ks 

Figure  4-2.  Bar  Math  Model 
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Figure  4-3.  Finite  Element  Idealization 


In  reality,  these  structures  are  continuous  systems  without  grid  boundaries.  The  finite 
element  method  is  in  the  form  of  a  mathematical  model  representing  a  continuous  struc¬ 
ture.  The  analyst  must  decide  upon  the  model  that  best  represents  the  structure  being 
examined  based  upon  the  type  of  analysis  desired,  the  time  available  and  the  experience 
of  the  analyst. 

One  very  important  advantage  of  the  finite  element  method  of  analysis  is  that  one 
model  can  often  be  used  to  perform  a  thermal  analysis  and  a  structural  analysis.  This 
dual  capability  makes  lais  technique  very  powerful  since  it  can  dramatically  increase  the 
productivity  of  the  design  engineer.  It  is  up  to  the  analyst  to  decide  how  to  properly  de¬ 
velop  the  finite  element  model  to  achieve  this  dual  capability.  Not  all  finite  element 
codes  are  capable  of  performing  both  thermal  analysis  and  structural  analysis  functions. 
Some  of  the  codes  that  have  dual  functions  are  ANSYS,  COSMOS  M,  STRUDL,  NISA, 
STARDYNE,  and  NASTRAN.  Programs  such  as  STARDYNE  and  STRUDL  only  have 
thermal  conduction  capability  while  others  also  have  convection  and  radiation  capability. 
However,  none  of  the  codes  have  the  capability  of  performing  transient  thermal  analyses 
where  physical  properties  are  mixed.  This  situation  can  occur  when  it  is  necessary  to 
change  radiation  and  convection  relations  as  a  function  of  altitude  and  temperature 
changes,  or  where  there  is  a  change  of  state  from  a  solid  to  a  liquid  that  requires  using 
the  latent  heat  of  fusion.  Problems  of  this  type  are  easily  solved  using  finite  difference 
methods. 
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4.1.2  Finite  Element  Computational  Process 

Most  computerized  finite  element  codes  break  up  the  computational  process  into  six 
steps  (Reference  4-1).  They  are: 

1)  Discretization  (or  representation)  of  the  given  domain  into  a  collection  of  prese¬ 
lected  finite  elements. 

a.  Construct  the  finite  element  mesh  of  pre-selected  elements. 

b.  Number  the  nodes  and  elements. 

c.  Generate  the  geometric  properties  (e.g.,  coordinates,  cross-sectional  areas,  etc.) 
needed  for  the  problem. 

2)  Derivation  of  element  equations  for  all  typical  elements  in  the  mesh. 

a.  Construct  the  variational  formulation  of  the  given  differential  equation  over  the 
typical  element. 

b.  Assume  that  a  typical  dependent  variable  u,  i.e.,  displacements,  is  of  the  form: 

n 

u  =  y  ui  Wi 

1=1  Eq.  14 

where  are  shape  functions  that  satisfy  the  equilibrium  and  boundary  condi¬ 
tions  across  the  element  and  substitute  it  into  step  2a  to  obtain  element  equa¬ 
tions  in  the  form: 

[K^]  {u<e>}  =  {F<e>}  Eq.  15 

c.  Derive  or  select,  if  already  available  in  the  literature,  element  interpolation 

functions  ip\  and  compute  the  element  matrices  where:  is  the  element  stif¬ 

fness  matrix;  F*e)  is  the  force  vector;  (e)  refers  to  the  element. 

3)  Assembly  of  element  equations  to  obtain  the  equations  of  the  whole  problem. 

a.  Identify  the  interelement  continuity  conditions  among  the  primary  variables  (re¬ 
lationship  between  the  local  degrees  of  freedom  and  the  global  degrees  of  free¬ 
dom  -  connectivity  of  elements)  by  relating  element  nodes  to  global  nodes. 

b.  Identify  the  equilibrium  conditions  among  the  secondary  variables  (relationship 
between  the  local  source  of  force  components  and  the  globally  specified  source 
components). 

c.  Assemble  element  equations  using  steps  3a  and  3b. 
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4)  Imposition  of  the  boundary  conditions  of  the  problem. 

a.  Identify  the  specified  global  primary  degrees  of  freedom. 

b.  Identify  the  specified  global  secondary  degrees  of  freedom  (if  not  already  done 
in  step  3b). 

5)  Solution  of  the  assembled  equations. 

6)  Postprocessing  of  the  results. 

a.  Compute  the  gradient  of  the  solution  or  other  desired  quantities  from  the  pri¬ 
mary  degrees  of  freedom  computed  in  step  5. 

b.  Represent  the  results  in  tabular  and/or  graphical  form. 

4.1.3  General  Considerations  and  Guidelines 

A  general  rule  of  thumb,  as  applied  to  finite  element  modeling,  is  that  the  more  e|p- 
ments  and  node  points  utilized,  the  more  accurate  the  model.  Many  text  books  show  sim¬ 
ple  problems  where  improved  accuracy  is  achieved  through  the  use  of  more  node  points 
and  some  text  books  show  that  the  stiffness  of  the  model  is  influenced  by  the  number  of 
nodes  and  elements,  where  the  fewer  the  number  of  nodes,  the  stiffer  the  model.  In 
general,  there  is  the  desire  to  use  more  node  points  to  improve  the  accuracy.  The  draw¬ 
back  to  using  more  node  points  is  the  rapid  increase  in  the  time  it  takes  to  solve  a  large 
finite  element  problem.  In  addition,  computer  memory  requirements  also  increase  rapid¬ 
ly  creating  problems  which  cannot  be  solved  on  a  small  machine.  A  common  practice  is 
to  use  a  coarse  mesh  in  areas  of  the  model  where  accuracy  is  not  important  and  to  use  a 
finer  mesh  where  accuracy  is  of  great  importance. 

An  understanding  of  the  structure  to  be  modeled  is  critical.  It  is  also  necessary  to  un¬ 
derstand  how  the  structure  is  expected  to  act  when  it  is  exposed  to  the  dynamic  environ¬ 
ment  ie.,  if  bolted  covers  are  used,  the  structure's  stiffness  will  be  altered.  Some  general 
considerations  in  the  FEM  techniques  are  outlined  below. 

1)  Understand  the  physical  characteristics  and  properties  of  the  individual  building 
block  elements  that  are  available  with  the  particular  FEM  program  being  used  to 
develop  the  model. 

2)  Visualize  the  geometric  shape  of  the  physical  system  and  the  general  behavior  un¬ 
der  the  action  of  applied  loads  and  restraints  that  may  distort  the  geometry. 
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3)  Determine  the  locations  of  areas  that  are  considered  to  be  critical  or  areas  where 
information  is  desired. 

4)  Select  the  geometry  that  properly  represents  the  physical  structure  and  its  behavior 
in  the  environment. 

5)  Select  the  type  of  element  (solid,  shell,  beam,  plane,  etc.)  that  best  represents  the 
characteristics  of  the  structure  being  analyzed. 

6)  Generate  the  mesh  density  that  will  obtain  the  desired  results  in  a  cost  effective 
manner  based  upon  the  capability  of  the  computer  being  used. 

7)  Apply  loading  and  boundary  constraints  that  are  consistent  with  the  geometry  and 
actions  of  the  physical  system. 

8)  Utilize  previous  testing  experience  or  modeling  experience  to  contribute  to  the  un¬ 
derstanding  of  the  load  path  through  the  structure.  Where  is  the  force  coming 
from?  Since  the  force  must  go  to  the  support  or  boundary,  how  does  it  get  there? 

9)  Make  use  of  symmetry  to  reduce  the  size  of  the  model. 

10)  Avoid  using  concentrated  loads  at  a  single  node  point.  This  can  lead  to  singulari¬ 
ties  where  stress  levels  are  much  higher  than  the  true  values. 

1 1 )  The  shape  of  the  model  should  approximate  the  shape  of  the  real  structure.  If 
stress  levels  in  a  fillet  are  desired,  then  a  fillet  should  be  included  in  the  model. 

12)  Nodes  shoaid  be  placed  at  load  points  and  at  support  poir:  s  and  anywhere  infor¬ 
mation  such  as  forces,  displacements,  and  stresses  are  desired. 

13)  A  large  mesh  pattern  is  desirable  in  open  sections  of  the  model  where  information 
is  not  critical.  A  gradual  change  should  be  made  to  smaller  mesh  sizes  using  a 
smooth  transition. 

14)  Good  stress  information,  in  general,  requires  a  finer  mesh  size  than  is  required  for 
displacement  or  resonant  frequencies. 

15  )  Avoid  meshing  elements  with  significantly  different  characteristics,  such  as  solids 
(which  do  not  bend)  with  beams  c  plates  that  can  bend  and  twist. 

lb)  Use  curved  elements  on  curved  surfaces  or  use  a  finer  mesh  with  flat  elements  at 
these  locations. 

17)  Avoid  high  aspect  ratio  (very  narrow  shell  or  plate  eien  ents).  Use  square  shell  or 
plate  elements. 
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18)  Gain  an  understanding  of  the  element  behavior  and  the  limitations  of  the  program 
being  used  for  FEM  analysis.  This  can  be  done  by  setting  up  classic  problems 
where  the  exact  answers  exist. 

4.2  Sources  of  Errors  in  a  Finite  Element  Analysis 

When  a  finite  element  solution  converges,  it  is  implied  that  the  exact  response  of  a  me¬ 
chanical  idealization  has  been  realized.  However,  the  finite  element  solution  is  only  an  ap¬ 
proximation  to  the  exact  response  and  different  sources  of  error  affect  the  finite  element 
solution  results.  Table  4-1  summarizes  the  various  sources  of  error  (Reference  4-2). 

TABLE  4-1.  FINITE  ELEMENT  SOLUTION  ERRORS 


Error 

Error  Occurrence  in 

Singularities 

Modeling 

Discretization 

Use  of  finite  element 

interpolation 

Numerical 

Evaluation  of  finite  element 

integration 

matrices  using  numerical 

in  Space 

integration 

Evaluation  of 

Use  of  nonlinear  material  models 

Constitutive 

Relations 

Solution  of  Dynamic 

Direct  time  integration, 

Equilibrium  Equations 

modal  superposition 

Solution  of  Finite 

Gausa-Seidel,  Newton-Raphson, 

Element  Equations 

quasi-Newton  methods, 

by  Iteration 

eigensolutions 

Round-off 

Setting-up  equations  and 

their  solution 

Singularities 

There  is  one  major  source  of  error  that  is  not  well  understood  by  most  analysts  which 
can  sharply  degrade  the  accuracy  of  any  FEM  analysis  -singularities.  Singularities  may 
produce  gross  errors  in  the  FEM  when  stress  calculations  are  requested.  Singularities  can 
be  found  at  a  point  source  of  heat,  a  point  load,  an  abrupt  change  in  the  boundary,  or  at 
a  sharp  corner  in  a  structure. 

Considering  a  point  load,  where  the  theoretical  area  under  the  point  is  zero,  the  result¬ 
ing  stress  level  is  infinite.  The  finite  elements  in  the  immediate  area  of  the  point  load 
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will  show  artificially  high  stress  levels  due  to  the  singularity  characteristics.  This  is  not  a 
real  condition  since  all  ductile  structures  will  deform  slightly  until  a  finite  area  is  formed 
that  will  support  the  load  or  the  structure  will  fail.  However,  the  computer  does  not  have 
this  information.  The  first  two  highest  stress  levels  are  often  discarded,  and  the  third 
highest  stresses  calculated  are  used  to  determine  the  pass  and  fail  criteria  for  the  struc¬ 
ture.  This  is  quite  arbitrary  and  not  an  acceptable  practice  for  evaluating  design  margins. 
The  accuracy  of  this  type  of  analysis  can  be  improved  by  using  a  nonlinear  model  where 
the  plastic  properties  of  the  material  are  included  in  the  comouter  model.  The  big  draw¬ 
back  to  nonlinear  models  is  the  great  amount  of  computer  time  that  is  required  to  obtain 
solutions.  The  analyst  must  decide  if  the  accuracy  improvement  is  worth  the  extra  com¬ 
puter  time  required  to  obtain  a  solution.  Another  way  to  improve  the  accuracy  of  the 
stresses  obtained  from  the  application  of  a  concentrated  load  is  to  alter  the  load  so  it  is 
not  applied  to  a  single  node  point.  Spreading  the  load  over  a  closely  clustered  group  of 
node  points  will  reduce  the  singularity  effects  and  reduce  the  peak  stress  levels. 

Geometric  shapes  can  also  result  in  singularities,  especially  at  holes,  notches,  angles, 
and  where  there  are  rapid  changes  in  the  cross  section  of  a  structure.  Consider  an  “L“ 
shaped  bracket  loaded  as  a  cantilevered  beam.  When  a  coarse  mesh  is  used  in  the  finite 
element  model,  the  inside  corner  of  the  “12’  will  be  sharp  with  zero  radius.  The  stress  lev¬ 
els  calculated  by  the  computer  at  this  inside  corner  will  be  much  higher  than  the  true  val¬ 
ue  for  the  same  reasons  outlined  above.  Normal  ductile  materials  will  simply  plastically 
deform  and  relieve  the  strain  so  the  part  will  not  fail.  However,  the  computer  does  not 
have  this  information  unless  a  nonlinear  analysis  is  performed.  One  way  around  this 
problem  is  to  add  a  small  radius  at  the  sharp  inside  corner,  which  will  reduce  the  magni¬ 
tude  of  the  maximum  stresses  to  a  more  realistic  value. 

Singularities  will  also  occur  at  the  boundaries  of  structural  members.  When  the 
boundaries  are  point  supports  or  clamped  edges,  the  conditions  can  be  treated  as  de¬ 
scribed  above  in  the  section  on  concentrated  loads  to  avoid  singularities. 
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4.2.2  Round-off  Errors 

Errors  critical  to  the  accurate  prediction  of  stress  in  any  general  engineering  structure 
such  as  finite  element  discretization,  have  been  discussed.  Errors  critical  to  the  finite  ele¬ 
ment  analysis  of  microelectronic  structures,  are  round-off  errors.  Round-off  errors  are  a 
result  of  the  finite  arithmetic  precision  of  the  computer  used.  Matrix  operations  are  typi¬ 
cally  the  culprit  in  introducing  solution  errors  due  to  round-off,  specifically  Guass  elimi¬ 
nation  (Reference  4-2).  Large  solution  errors  of  this  type  are  introduced  when  the 
diagonal  elements  in  the  stiffness  matrix  vary  by  a  large  amount,  or  when  very  smail  diag¬ 
onal  elements  are  used,  creating  a  large  multiplier.  The  reason  for  the  large  solution  er¬ 
ror  is  the  basic  mathematical  operation  in  Gauss  elimination  -  factorization. 

Factorization  is  a  subtraction  of  a  multiple  of  the  pivot  row  from  the  rows  below  it.  If 
numbers  of  widely  different  magnitudes  represented  to  a  fixed  number  of  digits  are  sub¬ 
tracted.  then  the  errors  in  this  operation  can  be  relatively  large.  In  other  words,  if  a 
structure  is  composed  of  many  different  materials,  such  as  a  component  lead/solder  con¬ 
nection,  then  the  soft  stiffness  of  the  solder  and  larger  stiffness  of  the  copper  lead  are 
represented  next  to  each  other  in  the  stiffness  matrix.  Hence,  large  and  small  stiffness 
values  are  present  on  the  main  diagonal.  Another  difficulty  is  analyzing  structures  where 
very  small  dimensions  are  being  modeled.  Since  refined  models  are  needed  to  ensure 
compatibility  and  completeness  of  the  solution,  truncation  errors  due  to  numerical  preci¬ 
sion  is  very  critical.  Therefore,  using  a  refined  mesh  can  cause  considerable  numerical 
errors  if  proper  precautions  are  not  taken  during  mesh  generation. 

4.3  Finite  Element  Accuracy  in  Stress  Analyses 

It  has  been  argued  that  small  feature  size  of  electronic  devices  makes  finite  element 
analyses  (FEA)  inappropriate  because  of  numerical  errors.  However,  these  errors  can  be 
minimized  by  following  some  general  guidelines  for  building  the  finite  element  model. 
The  following  examples  verify  numerical  errors  do  not  occur  when  small  dimensions  are 
used.  Additionally,  the  ability  of  FEA  to  simulate  predicted  stress  concentration  factors 
caused  by  geometric  changes  is  demonstrated.  Finally,  a  comparison  of  linear  and  nonlinear 
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analyses  demonstrates  the  variation  in  calculated  stresses  and  strains.  The  FEA  general 
purpose  code  ABAQUS  was  used  for  these  tests. 

When  utilizing  finite  element  analysis  techniques  for  small  scale  structures,  the  follow¬ 
ing  guidelines  must  be  adhered  to  in  order  to  achieve  realistic  results: 


1)  Quadrilaterial  elements  should  be  kept  as  square  as  possible.  The  aspect  ratios 
should  be  maintained  as  depicted  in  Figure  4-4.  In  cases  of  uniaxial  stress  fields, 
larger  aspect  ratios  are  acceptable. 


b 


a 

a/b  <  1 :4 


Figure  4-4.  Aspect  Ratio  for  Quadrilateral  Elements 


2)  When  analyzing  plates  and  shells,  be  sure  to  provide  a  sufficient  number  of  ele¬ 
ments  across  the  span  to  follow  the  deflection  surface  and  the  changes  in  shear. 
Remember  that  the  analysis  replaces  the  uniform  pressure  over  the  element  area 
by  a  set  of  equivalent  point  loads  applied  at  the  grids. 

3)  Finite  element  modeling  of  small  scale  structures  such  as  those  shown  in 
Figure  4-5,  requires  finite  elements  that  have  been  formulated  to  be  less  sensi¬ 
tive  to  curvature  or  skewed  geometries.  The  accuracy  of  typical  finite  elements 
based  on  simple  linear  shape  function  approximations  (typical  four  node,  2U 
degree-of-freedom,  quadrilateral  elements)  is  very  sensitive  to  modeling  geometries; 
for  instance: 

a.  It  is  essential  to  keep  the  relative  size  of  elements  in  areas  of  critical  stresses 
the  same. 

b.  Elements  must  also  be  relatively  square  in  shape  (aspect  ratio  equal  to  one)  in 
areas  of  critical  stress  prediction. 
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Avoid  Highly  Distorted  Elements 
0  Should  be  >  45° 


Elements  Around  a  Cutout 


1 5°  arc  or  Less 


Elements  Forming 
Cylindrical  Surfaces 


p<15°  4  Node  Quad 
P<30°  for  8  Node  Quad 


Figure  4-5.  General  Guidelines  for  Finite  Element  Modeling 
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4.3.1  Compatibility  and  Completeness  of  the  Finite  Element  Mesh 

As  a  general  rule  of  thumb,  the  more  elements  or  degrees  of  freedom  (DOF),  the 
more  accurate  the  solution.  However,  when  trying  to  decrease  model  size  but  still  retain 
accuracy,  a  coarse  mesh  can  be  used  in  areas  of  a  structure  where  stress  prediction  is  less 
critical.  For  the  finite  element  solution  to  converge  thru  mesh  refinement  (h-version 
FEA),  it  is  necessary  that  the  elements  that  make  up  the  assemblage  be  complete  and 
compatible.  The  requirement  for  completeness  of  an  element  means  that  the  displace¬ 
ment  functions  of  the  element  must  be  able  to  represent  the  rigid  body  displacements  and 
constant  strain  states.  The  refinement  of  a  finite  element  mesh  is  required  to  obtain  a 
constant  strain  state  within  an  element.  If  more  and  more  elements  are  used  in  the  as¬ 
semblage  to  represent  a  structure,  each  element  approaches  a  very  small  size  and  the 
strain  in  each  element  approaches  a  constant  value.  The  complex  variation  of  strain  with¬ 
in  the  structure  can  then  be  approximated.  The  requirement  of  compatibility  means  that 
the  displacements  within  the  elements  and  across  the  element  boundaries  must  be  contin¬ 
uous  since  the  stresses  in  an  element  are  calculated  using  derivatives  of  the  displacement. 
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When  calculated  in  adjacent  elements,  the  stress  may  vary  substantially  if  a  coarse  mesh  is 
used  or  if  adjacent  elements  are  not  the  same  size  (as  when  transitioning  from  a  coarse 
mesh  to  a  refined  mesh)  because  force  equilibrium  has  not  been  satisfied.  The  stress  differ¬ 
ence  at  the  element  boundaries  decreases  as  the  finite  element  idealization  is  refined,  and  in 
practice,  acceptable  results  are  usually  obtained  if  element  boundary  stresses  are  averaged. 

4.3.2  Small  Element  Size  Effects 


The  effect  of  micro-dimensioned  finite  element  models  on  stress  calculation  accuracy 
was  examined  for  a  plate  with  a  hole  in  tension.  The  model  in  Figure  4-6  shows  the 


ISOTROPIC  PHNEL  WITH  CIRCULRR  HOLE 


Figure  4-6.  Isotropic  Panel  With  Circular  Hole 

mesh  of  a  square  plate  with  a  distributed  load  along  the  top  and  bottom.  The  overall 
plate  size  was  varied  from  10  inches  to  0.0001  inches  in  length.  The  predicted  stresses 
were  equal  for  all  plate  dimensions  considered  with  a  gross  section  stress  concentration 
Kt,  of  3.370.  The  ratio  of  the  hole  diameter  (a)  to  the  plate  length  (r)  was  held  constant 
at  0.2  as  was  the  ratio  of  the  plate  length  {( )  to  plate  thickness  (t),  at  100.  The  loading 
was  proportionally  reduced  to  yield  the  same  nominal  stress.  anom,  which  is  defined  as: 
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Unom 


P 

i  t 


1.0 


Eq.  16 


where:  1  =  length 

t  =  thickness 
p  =  load 

Therefore,  each  finite  element  model  was  expected  to  predict  the  same  maximum 
stress  at  the  hole.  Figure  4-7  shows  the  typical  stress  distribution  around  the  hole.  The 
finite  element  results  are  shown  to  be  within  7%  of  the  3-D  elasticity  solution,  Kt  =  3.14 
(Ref.  4-3). 
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4.3.3  Beams  of  Two  Materials 

The  second  analysis  was  a  small  beam  made  of  two  different  materials.  The  beam  is 
illustrated  in  Figure  4-8.  The  concern  was  to  accurately  predict  the  stress  at  the  interface 
between  the  soft  and  stiff  materials.  A  stress  jump  (or  discontinuity)  will  always  exist  be¬ 
tween  two  adjoining  elements.  Flowever,  as  the  mesh  is  refined,  the  stress  jump  will  di¬ 
minish  (or  converge)  as  long  as  the  original  mesh  is  contained  within  the  new  mesh. 
Following  the  guidelines  mentioned  in  the  previous  section,  several  consecutive  refine¬ 
ments  in  the  mesh  are  made  to  determine  the  rate  of  convergence.  The  meshes  used  in 
this  study  are  depicted  in  Figure  4-9.  Each  mesh  is  a  double  refinement  of  the  previous 
mesh.  In  each  successive  mesh  refinement,  the  stress  jump  at  the  material  interface  was 
calculated.  Table  4-2  lists  the  stress  jumps  for  each  mesh.  As  a  rule  of  thumb,  a  stress 
jump  of  5  percent  across  the  element  boundary  is  considered  sufficient  refinement. 

The  mesh  used  did  not  have  any  abrupt  element  size  changes  and  it  did  not  have  any 
high  aspect  ratio  elements  in  the  area  of  high  predicted  stresses.  This  analysis  shows  FEA 
is  reliable  for  stress  prediction  of  objects  with  varying  elastic  modulus. 

Where 

•  L  *  0.09  in. 

1  5  •  h  «  0.01  in. 

•  w  =  0.01  in. 

•5  =  0.001  in. 

•  E1  =  17  E6  psi 

•  pi  ■>  0.31 58 

•  Eg  =  1.E6  psi 

•  p2  -  0.40 
GP03-0624-4-0/VS 

Figure  4-8.  Slender  Beam  With  TWo  Materials  and  Square  Cross-Section 

_ _ _  I  Mesh  1 


Figure  4-9.  Tip  Loaded  Cantilever  Beam  With  T\vo  Materials  -  Mesh 
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TABLE  4-2.  INTERFACE  STRESS  DISCONTINUITIES 


Mesh  Number 

Stress  Jump  (psi) 

Precent  Difference 

1 

1499 

10.8% 

2 

1313 

10.1% 

3 

722 

5.6% 

4.3.4  Comparison  of  Linear  and  Nonlinear  FEAs 

This  section  discusses  a  technique  which  will  allow  the  use  of  linear  FEAs  to  model  situa¬ 
tions  where  plastic  deformations  occur  in  the  structure.  Considerable  time  savings  can  be 
realized  with  this  technique  when  compared  to  using  nonlinear  (or  elastic-plastic)  FEAs. 

If  stresses  are  greater  than  the  yield  strength  (Oys)  of  the  material,  the  high  stress  sec¬ 
tion  of  the  component  is  in  a  plastic  state.  This  condition  is  characterized  by  nonlinear 
stress  (a)  versus  strain  (c)  behavior  is  illustrated  in  Figure  4-10a.  The  stresses  computed 
by  a  linear  FEA  analysis  are  not  valid  since  a  linear  relation  between  stress  and  strain  has 
been  assumed.  To  compute  the  actual  stress  and  strain,  a  nonlinear  elastic-plastic  FEA 
analysis  can  be  done,  but  this  is  usually  complex  and  time  consuming.  The  actual  stress 
and  strain  can  be  estimated  very  rapidly  by  using  the  stress  and  strain  computed  by  the 
linear  FEA  as  illustrated  in  Figure  4- 10b.  This  will  give  a  lower  (less  conservative)  esti¬ 
mate  of  the  actual  conditions.  Another  more  conservative  estimate  can  be  obtained  by 
using  Neuber’s  rule  for  notch  analysis  with  a  stress  concentration  of  1.0.  This  results  in: 

(ct€)actual  =  (oe)fEA  Eq.  17 

The  actual  stress  and  strain  combination  must  fall  on  the  curve  of  a  vs.  e  as  illustrated 
in  Figure  4-10c. 


1 


4-20 


Computer  Aided  Assessment  of  Reliability 
_ Using  Finite  Element  Methods 

A  simple  test  case  has  been  analyzed  for  the  purpose  to  verify  nonlinear  methods  for 
stress  prediction  in  microelectronic  structures.  A  two  dimensional  ABAQUS  model  has 
been  developed  for  elastic-plastic  analysis  and  was  compared  to  a  linear  ABAQUS  model 
for  stress  prediction.  The  linear  J-lead  analysis  in  Figure  4-11  demonstrates  how  high 
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Figure  4-10.  Stress  and  Strain  Approximations 
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Figure  4-1 1 .  Linear  J-Lead  FEA  Analysis 
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thermal  loads  can  cause  yielding  in  the  solder  connection  for  the  lead.  Stresses  in  the  sol¬ 
der  are  greater  than  the  yield  strength  (o^)  of  the  material,  as  shown  in  Figure  4-12.  The 
stresses  computed  by  the  linear  FEA  analysis  are  inaccurate,  since,  they  assume  a  linear 
relationship  between  stress  and  strain.  Flowever,  the  actual  stress  can  be  approximated 
by  dropping  vertically  down  to  the  material  stress  (o)  vs.  strain  (e)  curve  (Figure  4-12a). 
This  linear  correction  gives  a  lower  estimate  of  the  actual  condition.  A  more  conserva¬ 
tive  estimate  was  obtained  using  Neuber’s  rule  for  notch  analysis  with  a  stress  concentra¬ 
tion  of  one.  The  actual  stress  and  strain  predicted  by  the  ABAQUS  elastic-plastic 
analysis  falls  on  the  o  vs.  e  curve  between  the  linear  correction  and  Neuber's  rule.  For 
this  case  the  difference  in  the  linear  approximation  and  the  nonlinear  analysis  was  mini¬ 
mal.  In  general,  an  elastic-plastic  analysis  can  be  avoided  by  checking  the  difference  in 
strain  prediction  between  the  linear  correction  and  the  Neuber  approximation.  For  this 
case,  the  difference  in  strain  was  small  (  =  10%)  and  would  not  have  made  a  big  difference 
in  the  fatigue  life. 


•  Linear  FEA 
O  Linear  Approximation 
E  Elastic-Plastic 
X  Neuber  Rule 
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Figure  4-12.  Stress  -  Strain  Nonlinearity 
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For  the  case  when  the  degree  of  nonlinearity  is  great  (Figure  4- 12b),  then  an  elastic- 
plastic  analysis  is  critical  for  accurate  fatigue  life  predictions.  The  nonlinear  analysis  is 
necessary  in  this  case  due  to  the  large  horizontal  error  in  the  strain  prediction.  In  this 
case  the  strains  predicted  by  the  linear  approximations  disagree  by  more  than  50%.  This 
difference  can  relate  to  a  factor  of  5  to  10  difference  in  the  fatigue  life. 

For  both  test  cases,  the  vertical  error  in  stress  prediction  is  negligible  between  the  lin¬ 
ear  approximations  and  the  nonlinear  calculation.  This  shows  the  importance  of  using 
strain  versus  life  data  for  low  cycle  fatigue  (<  104  cycles)  as  opposed  to  stress  versus  life 
data  which  is  more  applicable  for  high  cycle  fatigue  calculations. 

4.4  Dynamic  Modeling/ Analysis  Technique  for  Electronic  Equipment 

Failures  that  occur  in  electronic  systems  during  exposure  to  vibration  environments  will 
usually  be  associated  with  PCBs  since  they  support  the  most  sensitive  elements  of  the  as¬ 
semblies.  These  sensitive  elements  include  the  electronic  components,  their  electrical 
leads,  solder  joints,  plated  through-holes,  and  interface  electrical  connectors.  When  the 
resonant  frequencies  of  the  PCB  are  excited  during  vibration,  the  PCB  will  bend  back 
and  forth  as  shown  in  Figure  4-13.  This  flexing  action  will  produce  relative  motion  be¬ 
tween  the  PCB  and  the  electronic  component  body.  When  the  component  has  leads  for 
electrical  interconnections,  the  load  path  will  pass  through  the  lead  and  into  the  solder 
joints,  the  plated  through-holes,  and  into  the  plastic  (usually  epoxy  fiberglass  or 
polyimide  glass)  PCB.  High  acceleration  vibration  levels  and  lightly  damped  PCBs  can 


Figure  4-13.  Printed  Circuit  Board  Out-of-Plane  Displacement 
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produce  large  dynamic  displacements  which  will  cause  rapid  fatigue  failures  in  the  leads 
and  solder  joints.  A  closer  examination  of  the  lead  bending  displacements  in  Figure  4-14 
shows  that  the  end  leads  have  the  greatest  movement  and  the  middle  leads  have  the  least. 


- Neutral  position 

— —  Deformed  shape  GP03-0624-5-D 


Figure  4-14.  Lead  Displacement 

There  are  four  main  components  of  the  electronic  assembly  which  must  be  modelled: 
the  chassis  (or  line  replaceable  unit,  LRU);  the  PCB;  the  component  and  leads;  and  the 
solder  joint  connection.  The  chassis  assembly  is  modeled  to  establish  the  magnitude  of 
the  dvnamic  coupling  response  between  the  PCci  and  LRU.  The  PCB  finite  element  model 
is  used  to  transfer  dynamic  displacements  and  strains  onto  refined  finite  element  models  of 
the  component,  lead  and  solder  joint.  Use  of  the  these  refined  models  is  restricted  to  deter¬ 
mining  the  stress/strain  levels  in  the  leads  and  solder  joints  for  fatigue  life  predictions. 

4.4.1  Determining  the  Chassis  to  PCB  Dynamic  Coupling 

Electronic  assemblies  typically  use  an  outer  housing  chassis  to  enclose,  support,  and 
protect  the  PCBs  (Figure  4-15).  These  enclosures  will  also  be  excited  during  the  vibra¬ 
tion  exposure.  If  care  is  not  exercised  during  the  design  phase  of  the  project  and  the 
resonant  frequency  of  the  outer  chassis  housing  is  close  to  the  resonant  frequency  of  one 
or  more  PCBs  supported  within  the  chassis,  dynamic  coupling  will  occur.  Under  these 
conditions,  the  response  accelerations  and  displacements  of  the  PCBs  can  be  further 
amplified  by  the  resonant  conditions  of  the  outer  chassis.  This  coupling  between  the 
chassis  and  the  PCB  can  sharply  reduce  the  fatigue  life  of  the  PCB. 


Computer  Aided  Assessment  of  Reliability 
Using  Finite  Element  Methods  _ 


Figure  4-15.  Electronic  Chassis  &  PCB  Assembly 


Therefore,  it  is  essential  to  establish  the  magnitude  of  the  dynamic  coupling  between 
the  LRU  and  PCB.  Hence,  the  finite  element  model  of  the  LRU/PCB  assembly  must 
consist  of  an  accurate  representation  of  the  LRU.  In  this  respect,  the  finite  element 
model  should  accommodate  cutouts,  concentrated  masses  (power  packs),  and  the 
compliance  between  the  PCB  and  the  side  wall  (edge  guides  and/or  wedge  clamps).  The 
PCB  finite  element  models  included  in  the  LRU  model  can  be  simplified  to  reflect  no 
mass  and  stiffness  effects  due  to  the  components.  The  modes  of  interest  in  determining 
the  dynamic  coupling  are  the  local  panel  modes  of  the  LRU  side  walls  and  the  primary 
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resonant  modes  (support  modes,  and  chassis  torsion  and  bending  modes).  If  the  LRU 
chassis  is  too  stiff  (which  could  be  the  situation  when  not  all  mass  is  represented  or  when 
the  support  mounts  on  the  chassis  are  not  accurately  represented),  lower  than  expected 
PCB  vibration  levels  will  result  by  missing  the  coupling  effect  in  the  analysis.  Once  the 
FEM  model  is  complete  (Figure  4-16),  the  computer  can  be  used  to  determine  the 
frequency  response  of  the  assembly.  This  will  show  the  magnitude  of  the  acceleration 
levels  at  every  node  point  in  the  model  across  a  broad  frequency  band,  typically  from  10 
Hz  to  2000  Hz.  The  computer  can  show  direct  axis  response,  where  a  1.0  G  dynamic 
stimulus  is  input  along  the  X  axis  and  the  response  is  measured  along  the  X-axis.  The 
computer  can  also  show  cross  axis  response,  where  a  1.0  G  dynamic  stimulus  is  input 
along  the  X  axis  and  the  response  is  measured  along  the  Y  or  Z  axis. 


GP03-0624-6/gms 

Figure  4-16.  Chassis  Assembly  FEM 

Steinberg  (Ref  4-4)  has  developed  a  simple  method  to  approximate  the  transmissibility 
factors  applicable  to  the  dynamic  response  of  PCBs.  If  the  chassis  (LRU)  frequency  (tc)  and 
the  PCB  frequency  (fp)  are  known,  the  uncoupled  PCB  transmissibility  is  given  by  equation: 
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Q  =  1  /  (1  -  R2) 


Eq.  18 


\>here:  R  =  ff/tc 

ff  =  forcing  frequency 

When  the  forcing  frequency  is  equal  to  the  natural  frequency  of  the  PCB(fp  =  ff),  the 
coupled  transmissibility  is: 

Qp  =  Q*SQRT(fp)  Eq.  19 

Knowing  the  damping  of  the  PCB  and  the  acceleration  levels  of  the  chassis,  the 
coupled  transmissibility  can  be  determined.  Qp  accounts  for  the  additional  energy  from 
the  chassis  due  to  the  coupling  with  the  chassis  at  the  ff  and  is  the  amplification  factor 
applied  to  the  dynamic  loads  that  excite  the  PCB  at  the  frequency  fp. 

4.5  Techniques  for  Modeling  the  PCB 

Once  the  LRU  to  PCB  dynamic  load  transfer  is  established,  the  PCB  finite  element 
model  can  be  analyzed  in  more  detail  by  itself.  The  individual  PCB  finite  element  model 
is  used  to  establish  the  mode  shapes  and  to  determine  which  components  are  critical  in 
terms  of  relative  displacement  to  the  PCB.  Determining  which  components  are  most  crit¬ 
ical  is  not  a  straight  forward  task.  Experience  and  empirical  techniques  will  minimize  the 
critical  components  which  need  to  be  analyzed  in  further  detail.  The  PCB  is  not  a  critical 
component  by  itself  since  very  few  failures  will  occur  in  the  board  before  the  components 
on  the  board  fail  (Ref  4-5.). 


The  PCB  should  be  analyzed  in  detail  for  possible  failure  of  the  etched  copper  traces 
and  for  failures  of  plated  through  holes  (PTH).  It  is  also  necessaiy  to  predict  the  individ¬ 
ual  layer  stress  distributions.  This  requires  that  each  layer  be  modeled  with  solid  finite 
elements  or  composite  plate-type  elements  where  a  composite  laminate  failure  theory  can 
be  applied.  More  commonly,  PCBs  are  modeled  with  plate-type  elements  (Figure  4-17) 
to  determine  dynamic  displacements.  Boundary  conditions  of  the  PCB  can  generally  be 
2'  med  to  be  simply  supported  at  the  mother  board  connection,  unsupported  at  the  top, 
simply  supported  on  the  -;dcr.  "’hen  edge  guide.,  (r;rwher)  are  used.  nr  simply  supom  > ;  J 
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SOLID-TYPE 


Figure  4-17.  Finite  Elements  for  Structural  Analysis 


with  rotation  flexibility  when  wedge  guides  are  used.  Steinberg  (Reference  4-4)  discusses 
the  rotational  flexibility  of  wedge  clamp  devices.  He  states  that  wedge  clamps  act  more 
like  a  fixed  boundary  condition  when  a  PCB’s  fundamental  natural  frequency  is  less  than  100 
Hz  and  more  like  a  simply-supported  edge  when  fundamental  natural  frequency  is  above 
600  Hz.  Therefore,  if  the  fundamental  frequency  is  between  100  and  600  Hz,  the  rotational 
flexibility  or  spring  stiffness  of  the  wedge  guide  can  be  determined  by  using  the  relationship: 


f,  +  1-10  (ft  -  F,) 
1  +  0.001  (F,  -  f.) 


Eq.  20 


where  fn  =  expected  natural  frequency  of  PCB 

ff  =  natural  frequency  of  PCB  with  fixed  (clamped)  sides 
fs  =  natural  frequency  of  PCB  with  simply-supported  sides 

The  frequencies  ff  and  fs  are  determined  from  the  F.E.A.  for  a  PCB  with  both  fixed 
sides  and  simply-supported  sides.  Then,  the  expected  natural  frequency,  fn,  is  calculated 
using  the  above  equation.  Using  the  finite  element  model,  rotational  springs  are  used  to 
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model  the  wedge  guides  rotational  flexibility.  Therefore,  several  F.E.A.s  are  performed 
at  different  spring  constants  until  the  F.E.A.  natural  frequency  matches  the  expected 
natural  frequency  calculated  using  the  above  equation.  In  other  words,  the  finite  element 
model  must  be  “tweaked”  to  more  closely  match  a  real  structure’s  response. 

Component  modeling  is  usually  restricted  to  concentrated  mass  representation  in  PC B 
finite  element  models.  Bivens  (Ref  4-6.)  showed  that  the  leaH  design  had  a  negligible  ef¬ 
fect  on  the  expansion  of  the  component  package  and  the  PCB  as  a  result  of  thermal  de¬ 
flections.  A  similar  behavior  can  be  assumed  for  dynamic  deflections. 


When  performing  a  dynamic  response  analysis  of  the  individual  PCB,  the  symmetry  of 
the  structure  can  be  used  to  reduce  analysis  cost  and  effort.  For  plate  type  structures, 
analyzing  a  quarter  symmetry  model  is  limited  to  a  plate  with  similar  boundary  conditions 
on  all  four  sides.  Typically,  PCBs  exhibit  only  half  symmetry  (symmetric  in  only  one  direc¬ 
tion).  For  a  one  half  symmetry  model,  it  is  necessary  to  perform  the  eigenvalue  extrac¬ 
tion  two  times  in  order  to  extract  both  the  symmetric  and  anti-symmetric  modes. 


Symmetric  boundary  conditions  about  the  X  plane  are  specified  with  ux  =  ey  =  ez  =  0- 
Anti-symmetrie  boundary  conditions  can  be  specified  with  uy  =  ez  =  0x  =  0  (Figure  4-18). 
As  is  typical  with  PCB  structures,  the  fundamental  resonance  contributes 


PCB  is  symmetric  about  the  y-z  plane. 


Boundary  conditions  along  line  of  symmetry: 
SYMMETRIC- 


Ux=  0Y=0Z=O 

ANT1-SYMMETRIC- 

UY=  9X=9Z=0 


where: 

U  -  displacement  degree  of  freedom 
6  •  Rotational  degree  of  freedom 


Figure  4-18.  Printed  Circuit  Board  Boundary  Conditions 
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Half  Finite  Element  ModelMode1 


Symmetric  Modes 


598.6  Hz 


Anti-symmetric  Modes 


Figure  4-19.  Half  Model  Finite  Element  Model  Resonant  Modes 


the  most  damage  to  the  components  and  can  be  predicted  with  only  the  symmetric 
boundary  condition.  Note  that  only  symmetry  to  the  stiffness  matrix  is  assumed.  When 
making  refined  PCB  models  which  include  the  component  mass,  the  finite  element  mass 
matrix  may  no  longer  be  symmetric  and  will  cause  errors  in  the  frequency  prediction. 

As  an  example,  the  natural  modes  of  vibration  in  a  half  symmetric  PCB  model  are 
shown  in  Figure  4-19.  The  frequencies  and  modes  shown  are  for  both  the  symmetric  and 
antisymmetric  boundary  conditions.  These  are  in  complete  agreement  with  the 
frequencies  and  modes  extracted  in  the  full  model  analysis  shown  in  Figure  4-20. 
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Full  Finite  Element  Model 


598.6  Hz  689.3  Hz 


The  Natural  Modes  of  Vibration  for  a 
Clamped-Clamped-Pinned-Free  Plate 

Figure  4-20.  Full  Finite  Element  Model  Resonant  Modes 

To  accurately  predict  stresses  in  lead  geometries  from  a  dynamic  model  of  the  PCB,  it 
is  necessary  to  model  the  components  individually.  Detailed  solid-type  finite  elements 
are  used  throughout  the  board,  ceramic  case,  and  leads  to  account  for  the  small  lead 
geometries  (Figure  4-21).  This  requires  extensive  modeling  and  an  exorbitant  number  of 
degrees  of  freedom.  To  reduce  the  size  of  the  finite  element  model,  simplified  models 
can  be  assumed  using  plate-type  elements  for  the  chip  carrier  and  beam-type  elements 
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Figure  4-21.  Small  Lead  Geometries 


tor  the  copper  leads  as  described  by  Soovere  et.  al,  (Ref.  4-5).  The  modeling  described 
lumps  groups  of  four  to  five  of  the  side  leads  into  a  single  beam  placed  at  the  center  of 
the  group  pin  location  (Figure  4-22).  The  corner  leads,  which  are  the  most  critical,  are 
individually  modeled.  The  loads  predicted  by  fin..e  element  analysis  in  the  beam  elements 
can  then  be  applied  to  material  failure  models. 


GP03-0624-7-0/qg 

Figure  4-22.  Lead  Grouping 
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4.6  Transferring  Vibration  Loads  and  Displacements  onto  Detailed  Finite 
Element  Model  of  Component  Leads 

The  following  example  demonstrates  how  loads  from  a  PCB  level  analysis  are  t  rns- 
terred  to  a  detailed  lead  analysis.  The  PCB  configuration  is  depicted  in  Figuie  4-23. 

The  out-of-plane  displacements  due  to  random  vibration  were  determined  from  a  pre¬ 
vious  analysis  to  be  C.0028  inches  rms  at  the  chip  center.  Figure  4-24.  Loads  and  dis¬ 
placements  in  the  leads  were  found  using  the  three-dimensional  PCB 'chip  model  in 
Figure  4-25.  It  consists  of  a  PCB  and  chip  carrier  modeled  with  plate-type  elements  and 
leads  individually  modeled  with  bar-type  elements.  The  PCB  is  clamped  on  the  sides, 
simply  supported  on  the  top  and  free  on  the  bottom.  The  model  was  only  used  to  predict 
defections;  hence,  the  meshing  lackc  the  necessary  refinement  to  accurately  predict 
stresses.  The  deflections  in  the  lead  wires,  extensional  and  rotational,  were  determined 
for  an  applied  displacement  of  0.0028  inches  at  the  chip  center. 


Hinge  Support 
Board 

Dimensions:  5.65  in.  x 


Chip 

•  Dimensions:  0.95  in.  x 
0.95  in.  x  0. 10  in. 

•  Material:  Ceramic 
•Weight:  0.007827  1b 

•  Description:  68  Lead, 
Chip  Carrier 

•  CTE:  6.4  ppm/deg  C 

•  Modulus:  45E6  lb/in.**2 


6.14  in.  x  0.064in. 

Material:  Sandwich  of 
Poly/Glass  and  Copper 
Weight:  0.327  lb 
Copper 

-Modulus:  16.0L6  lb/in**2 

-p:  0.34 

Poly/Glass 

-  Modulus:  2.0E6  lb/in**2 
-p:  0  34 

GP03-0624-10-D'ks 


Figure  4-23.  PCB  Configuration 
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Figure  4-24.  Out-of-Plane  Displacements 


PCB  LEVEL 
MODEL 


Figure  4-25.  Three  Dimensional  FEA  Model 
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For  this  detailed  PCB  level  analysis,  the  leads  were  crudely  modeled  with  bending  bar- 
type  elements  which  have  an  approximate  equivalent  stiffness  (keq)  to  the  J-lead  depicted 
in  Figure  4-26.  Typically,  properties  of  finite  element  bar  elements  include  elastic  modu¬ 
lus  (E),  width  (VV),  thickness  (Th)  and  length  (L).  The  axial  stiffness  is  given  by  Eq.  21, 

K«,  (axial)  =  (E)*(W*Th)/(L)  Eq.  21 

The  equivalent  axial  stiffness  of  the  J-lead  shown  in  Figure  4-26  is  given  by  Eq.  22, 

Keq  (axial)  =  [(1/KO  +  (1/K2)  +  (l/Kj)]'1  Eq.  22 

where:  =  SE^i^LO3 

K2  =  EiA2/L2 
Kj  -  ESA,/U 


Therefore,  the  elastic  modulus  used  in  the  three  dimensional  lead  elements  was  deter¬ 
mined  from  Eq.  23, 


Eeq  =  (Keq  *  L)  /  (W*Th) 


Eq.  23 


Node 

-<?- 


Ear  Element 
With 

Properties 
(W,  Th.  E.  L) 


— O — 
Node 


J-Load  Configuration 

GP03-0624- 1 2  0/5U2 


Finite  Element  Representation 
In  PCB  Level  Analysis 


Figure  4-26.  J-Lead  Axial  Stiffness  Representation 
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From  the  three  dimensional  FEA,  axial  and  bend;ng  loads  were  predicted  for  the 
leads.  The  maximum  loads  predicted  were  in  the  corner  lead  in  the  lower  left  hand 
side.  This  is  due  to  the  relative  displacement  of  the  chip  with  respect  to  the  PCB  as 
shown  in  Figure  4-7.7.  These  loads  were  then  imposed  as  enforced  displacements  on 
a  detailed  three  dimensional  finite  element  model  of  the  J-lead  (Figure  4-28).  The 
axial  displacement  and  moment  displacement  were  determined  from  Eqs.  24  and 
25  respectively. 


Faxiai/Keq  =  A*'31  Displacement  (yaxial)  Eq.  24 

XTanG  =  Bending  Moment  Displacement  (ybending)  Eq.  -5 

where:  0  =  M*L/(2EI)  Eq.  26 

and  are  shown  in  Figure  4-29. 

CHIP  PLANE 


PCB  PLANE  -  S  Largest  Relative  Deflection 


Figure  4-27.  Relative  Displacement  -  Chip  to  PCB 
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Figure  4-28.  Detailed  Three-Dimensional  J-Lead  FEM 


Node 


PCB  Level  Representation 


Axial  Load  Moment  Load 


' r  J  If  ' 1 


F/KEg 

Detailed  Lead  Model  Combined  Loading  Scheme 

G  P03-0624- 1  3-D/suZ 


Figure  4-29.  Axial  and  Moment  Displacement 
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The  results  of  the  load  transfer  analysis  are  shown  in  Figures  4-30,  4-31,  and  4-32. 
The  deformation  is  illustrated  in  Figure  4-30.  The  peak  load  stress  is  15,700  psi  in  the 
top  horizontal  arm  of  the  lead  shown  in  Figure  4-31.  The  peak  solder  stress  is  1095  psi 
at  the  lead/solder  interface  shown  in  Figure  4-32. 


Figure  4-30.  Load  Transfer  Analysis  Results  -  Deformation 
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Figure  4-32.  Solder  Stress 
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4.7  Modeling  Consideration  for  Leads  and  Solder  Joints 

Several  modelling  and  analysis  techniques  have  been  used  to  predict  solder  joint 
stresses  (References  4-6,  4-7,  4-8,  4-9,  4-10  &  4-11).  Modelling  techniques  employing 
both  two  and  three  dimensional  assumptions  have  been  demonstrated.  Two  dimensional 
models  require  plane  stress  assumption  in  the  state  of  stresses  in  the  lead.  The  various 
three  dimensional  modelling  techniques  in  the  literature  include  partial  models  of  the 
chip  carrier,  lead,  so'der  joint  and  PCB.  The  chip  carrier  and  PCB  will  be  represented 
based  on  the  type  of  component  connections  (surface  mounted  or  poke  thru).  The  fol¬ 
lowing  discussion  will  focus  on  modelling  techniques  which  illustrate  the  failure  mecha¬ 
nisms  of  components  with  leads  and  without.  The  examples  will  include  rationale  and 
explanation  for  any  modelling  assumptions. 

4.7.1  Plated  Through  Hold  (PTH)  Lead  and  Solder  Joint  Failures 

PTH  failures  normally  occur  during  thermal  cycling  environments  where  expansion  of  the 
PCB  in  the  Z-axis  produces  high  tensile  forces  and  stresses  in  the  copper  barrel  of  the  PTH 
as  shown  in  Figure  4-33.  The  PCB  expands  more  than  the  copper  in  the  PTH  since  the  co¬ 
efficient  of  thermal  expansion  (CTE)  of  the  PCB  is  much  higher  than  the  CTE  of  the  cop¬ 
per.  When  the  number  of  stress  cycles  and  the  corresponding  stress  levels  are  high  enough, 
a  fatigue  fracture  can  occur  in  the  copper  barrel  and  produce  an  electrical  malfunction. 


z 


GP03-0624-8-D/SUZ 


Figure  4-33.  Plated  Thru  Hole  Z-Axis  Expansion 
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Two  different  FEM  techniques  can  be  used  to  determine  the  thermal  expansion 
stresses  induced  in  the  structural  elements;  two  dimensional  (2-D)  and  three  dimensional 
(3-D)  modelling.  Two  dimensional  modelling  can  be  accomplished  using  plate-type  ele¬ 
ments  where  only  the  thermal  expansion  forces  are  acting;  no  externa!  forces  are  acting 
on  the  electrical  leads.  It  is  necessary  to  model  a  segment  of  the  connection  including  the 
lead,  solder  joint,  copper  barrel,  and  PCB  (Figure  4-34).  The  differences  in  the  CTE  of 
the  materials  will  generate  the  internal  stress  field.  Three  dimensional  models  must  be 
constructed  with  isoparametric  (solid)  elements  to  obtain  the  full  3-D  characteristics. 
Generally,  the  3-D  model  represents  the  actual  structure  more  accurately. 

Most  FEM  codes  have  an  eight  node  isoparametric  element  and  a  twenty  node  isopa¬ 
rametric  element  available.  The  eight  node  solid  contains  a  linear  stress  distribution 
within  the  element.  Therefore,  when  the  stress  distribution  through  a  structure  is  known 
to  be  linear,  the  use  of  the  eight  node  element  will  save  a  considerable  amount  of  com¬ 
puter  time  because  there  are  far  fewer  node  points  in  the  model. 


GP03-0624-9D/gms 

Figure  4-34.  Plated  Through  Hole  Section 
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When  the  boundary  conditions  are  properly  defined,  it  is  often  easier  to  slice  a  section 
out  of  the  3-D  model,  as  shown  in  the  Figure  4-35.  To  simplify  the  model  by  reducing 
the  number  of  solid  elements,  a  30°  slice  is  removed  from  the  fall  model.  This  section 
examines  the  stress  distribution  in  the  30°  slice  (cyclic  symmetry  has  been  assumed). 

A  relatively  simple  model  can  be  constructed  if  the  magnitude  of  the  external  lead 
loads  are  known  from  prior  experiments.  Very  often  the  magnitude  of  the  external  loads 
are  not  known,  so  the  model  developed  must  include  the  structural  elements  responsible 
for  producing  the  critical  loads. 


Concentrated  Load.  No  Fillet 


GP03  0183  35  0 

Figure  4-35.  Three  Dimensional  Model  Slice 
4.7.2  Thermal  Expansion  Stress  Analysis  in  DIPs 

Another  type  of  thermal  expansion  problem  that  can  be  solved  using  FEM  techniques 
involves  long  components  such  as  DIPs  and  hybrids.  A  typical  quarter  model  of  a  DIP  on 
a  PCB  is  shown  in  Figure  4-36.  The  DIP  component  has  a  much  lower  CTE  than  the 
PCB  in  the  x-y  plane  of  the  PCB.  At  the  high  temperature  end  of  the  thermal  cycle,  the 
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PCB  expands  more  than  the  component.  This  produces  relative  bending  in  the  electrical 
leads  as  shown  in  the  schematic  in  Figure  4-37.  In  Section  4.8,  examples  show  how  the 
loads  on  the  lead  wires  can  be  determined.  Methods  for  modelling  the  electrical  leads  on 
the  DIP  are  shown  in  Figure  4-38. 

The  leads  generate  a  shear  tearout  stress  in  the  PCB  solder  joint  due  to  the  bending 
moment  induced  by  the  relative  expansion  difference  between  the  component  and  the 
PCB.  The  solder  joint  should  not  be  allowed  to  exceed  a  true  value  of  about  400  psi  to 
insure  a  15  to  20  year  life  in  an  environment  where  extensive  thermal  cycling  is  expected. 


Figure  4-36.  Quarter  Model  of  DIP  on  PCB 
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Figure  4-37.  Thermal  Expansion  in  Plane 


Figure  4-38.  Modelling  DIP  Leads 


Some  typical  methods  for  modelling  the  leads  and  solder  joints  for  a  DIP  are  shown  in 
Figure  4-39.  Many  different  combinations  of  isoparametric  solid  elements,  shell  ele¬ 
ments  and  beam  elements  can  be  combined  to  develop  the  model. 
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4.7.3  Thermal  Stress  Analysis  in  LCCCs 

Leadless  ceramic  chip  carrier  (LCCC)  components  are  being  used  in  great  quantities 
for  surface  mount  PCBs.  The  ability  of  these  components  to  survive  extended  thermal 
cycling  environments  is  determined  by  the  solder  joint  that  supports  these  devices.  The 
solder  point  integrity  is  directly  related  to  the  solder  joint  fatigue  properties  and  the  qual¬ 
ity  of  the  process  that  manufactures  the  solder  joints.  Since  the  process  controls  are  not  a 
design  function,  the  emphasis  here  is  directed  towards  a  better  understanding  of  the  creep 
characteristics  and  fatigue  properties  of  the  solder  joints  at  various  temperatures  and 
thermal  cycling  conditions. 

The  most  common  solder  used  today  is  a  eutectic  solder  with  63%  tin  and  37%  lead 
and  a  melting  point  of  about  184°C.  Its  tensile  strength  and  its  modulus  of  elasticity 
show  an  increase  with  an  increase  in  the  loading  rate.  These  same  factors  show  a  de¬ 
crease  with  increasing  temperatures.  The  stress-strain  curve  is  highly  non-linear,  as  shown 
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in  Figure  4-40,  and  its  modulus  of  elasticity  rapidly  decreases  with  increasing  temperature 
as  shown  in  Figure  4-41.  Solder  is  highly  strain  rate  sensitive  as  shown  in  Figuie  4-42. 


Figure  4-40.  Highly  Nonlinear  Stress-Strain  Curve 


Temperature  -  °C 
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Figure  4-41.  Modulus  Change  With  Temperature 
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Figure  4-42.  Strength  vs.  Strain  Rate 


It  tends  to  exhibit  stress  relaxation  at  high  temperatures,  above  100°C.  Solder  will  also 
rupture  from  creep  under  the  action  of  a  steady  applied  load  for  extended  periods,  as 
shown  in  Figure  4-43.  .Another  unusual  characteristics  of  eutectic  solder  is  that  its  stress 
fatigue  life  is  directly  related  to  frequency  of  the  alternating  load  as  well  as  the  tempera¬ 
tures.  as  shown  in  Figure  4-44. 

Section  4.8  will  illustrate  typical  models  used  for  LCCC  devices  and  leaded  surface 
mount  components. 
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4.8  Thermal  Stress  Solder  Joint  Finite  Element  Analysis 

This  section  describes  thermal  stress  analyses  of  a  J-lead  solder  joint  and  a  leadless 
chip  carrier  solder  joint.  The  J-lead  problem  was  analyzed  using  five  different 
procedures.  The  first  and  second  procedures  consist  of  a  two  dimensional  (2-D)  Nastran 
and  a  2-D  ABAQUS  analysis.  The  third  procedure  used  a  similar  2-D  Probe  analysis. 
The  fourth  procedure  consists  of  a  2-D  Probe  analysis  of  the  joint  and  the  chip  carrier, 
and  finally,  the  fifth  procedure  consists  of  a  2-D  Probe  analysis  of  the  joint,  the  chip  car¬ 
rier  and  the  PCB.  The  fifth  part  examines  the  effect  of  allowing  thermal  expansion  in 
both  the  vertical  and  horizontal  directions. 

The  second  example  is  a  linear  elastic  finite  element  analysis  (FEA)  of  a  leadless  sol¬ 
der  joint  and  compares  four  different  procedures.  The  first  procedure  used  a  2-D 
ABAQUS  model  of  the  solder  joint  subjected  to  thermal  displacements.  The  second 
approach  consists  of  a  2-D  Probe  analysis  of  the  solder  material.  The  third  approach 
consists  of  a  2-D  Probe  analysis  of  the  solder  joint  and  the  chip  carrier.  The  fouth  pan 
consists  of  a  2-D  Probe  analysis  with  the  chip  carrier  and  an  applied  temperature  differen¬ 
tial.  which  examines  the  effect  of  allowing  thermal  expansion  in  both  the  horizontal  and 
vertical  directions. 

These  examples  address  the  accuracy  of  the  stress/strain  analysis  of  the  lead/solder 
joint.  Nastran  and  ABAQUS  are  h-version  finite  element  codes.  Probe  is  a  p-version 
finite  element  code  from  the  MacNeal-Schwendler  Corporation.  In  Nastran  and  ABA¬ 
QUS,  the  discretization  error  is  controlled  by  mesh  refinement.  In  Probe,  the  discretiza¬ 
tion  error  is  controlled  by  increasing  the  polynomial  degree  of  the  interpolation  function 
and/or  mesh  refinement  (Figure  4-45).  When  using  Probe,  the  relative  error  in  energy 
can  easily  be  established  by  performing  successive  analysis  at  three  different  p-levels. 
Calculating  the  relative  error  in  the  strain  energy  and  monitoring  the  convergence  of 
functionals  across  like  element  boundaries  ensures  the  quality  of  the  finite  element 
analysis.  The  functionals  include  stresses,  strains,  and  stress  resultants. 
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An  additional  quality  control  benefit  of  using  p-version  technology  is  the  minimization  of 
round-off  error,  of  concern  because  of  the  very  small  elements  used  to  model  the  joints. 
These  quality  control  procedures  are  integral  to  extracting  the  location  and  magnitude  of 
stresses/strains  for  use  in  fatigue  life  predictions. 


WHAT  IS  MEANT  BY  P-VERSION  AND  H-VERSION 


P-version  of  the  finite  element  method 


P'  is  Ihe  first  letter  ol  Ihe  word 
polynomial. 

In  Ihe  p-version  the  polynomial  order 
of  the  approximating  function  is  increased 
progressively  In  order  to  get  better 
results.  {The  FE  mesh  remains  fixed.) 


H-version  of  the  finite  clement  method 


H'  Indicates  the  size  of  Ihe  element. 
The  mesh  gets  progressively  denser 
(the  elements  gel  smaller)  In  order 
to  get  better  results. 


I 

T 


Figure  4-45.  Discretization  Error  Control 
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4.8.1  J-lead/Solder  Joint  Finite  Element  Model  Geometry 

The  2-D  plane  stress  finite  element  models  (FEM)  of  the  J-lead/solder  joint  are 
shown  in  Figures  4-46  and  4-47.  The  basic  J-lead/solder  FEM  consists  of  the  following: 


Q)  NRSTRRN  MODEL  b)  PROBE  MODEL 


Figure  4-46.  J-Lead  /  Solder  Joint  Geometry 
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Two  Dimensional  Linear  Analysis  of  a  J— Lead 


Figure  4-47.  ABAQUS  FEM 

1)  The  Kovar  J-lead  material  is  0.008  inches  thick. 

2)  The  height  of  the  J-lead/solder  joint  is  0.125  inches. 

3)  The  chip  carrier  is  0.084  inches  thick. 

4)  The  PCB  is  0.067  inches  thick. 

5)  All  finite  element  models  have  a  membrane  thickness  of  0.012  inches  for  the  plane 
stress  analysis. 

6)  Material  properties  for  the  various  materials  used  in  the  J-lead/solder  joint  are 
shown  in  Table  4-3. 

7)  The  60Sn/40Pb  solder  material  properties  are  at  room  temperature  conditions. 

8)  Symmetric  boundary  conditions  are  used  so  that  only  one  half  of  the  chip  carrier 
and  PCB  are  modeled. 

TABLE  4-3.  MATERIAL  PROPERTIES 


Material 

Elastic 

Modulus 

(psi) 

Poisson’s 

Ratio 

Thermal 
Expansion 
(in/ln/deg  C) 

Kovar 

2.E+07 

0.31 

5.5E-06 

60/40  Solder 

2.E  +  06 

0.40 

28.0E-06 

Chip  Carrier 

54.E  +  06 

0.25 

6.4E-06 

PCB 

2.E  +  06 

0.15 

9.5E-06 
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4.8.1. 1  ABAQUS  &  Nastran  Finite  Element  Analysis  of  the  J-lead/Solder 

The  Nastran  FEM  is  shown  in  Figure  4-46a  and  consists  of  78  Quad4  membrane  lead 
elements  and  56  Quad4  membrane  solder  elements.  The  ABAQUS  FEM  is  shown  in 
Figure  4-47  and  consists  of  228  lead  elements  and  216  solder  elements. 

The  top  flange  of  the  J-lead  in  both  models  is  fixed  in  both  the  X  (horizontal)  and  Y 
(vertical)  directions.  The  applied  loading  represents  the  differential  thermal  expansion  in 
the  plus  X  direction,  between  the  chip  carrier  and  PCB,  due  to  a  power  up  cycle.  This 
plus  X  expansion  will  be  called  the  PXTHERM  load  condition.  The  temperature  differ¬ 
ential  (AT  =  Tmax  -  Tmin)  of  the  chip  carrier  was  assumed  to  be  53  degrees  C  and  the 
temperature  differential  of  the  PCB  was  assumed  to  be  46  degrees  C  for  the  power  up 
cycle.  The  loading  was  applied  to  the  FEM  using  imposed  displacements  in  the  X  direc¬ 
tion  along  the  base  of  the  solder  as  shown  in  Figure  4-48a.  At  the  left  comer  of  the  sol¬ 
der  base  the  imposed  displacement  is  +  .319E-04  inches  in  the  X  direction.  At  the  right 
corner  of  the  solder  base  the  imposed  displacement  is  +  .625E-04  inches  in  the  X  direc¬ 
tion.  These  displacements  were  calculated  using  the  following  formula: 

Ax  =  Cl  *  (x  +  C4)  -  C2  *  C3.  Eq.  27 

where:  Cl  =  4.37E-4  in/in,  PCB  thermal  strain 
C2  =  3.392E-4  in/in,  chip  thermal  strain 
C3  =  .483  in.,  one  half  the  chip  length 
C4  =  .448  in.,  center  of  the  chip  to  beginning  of  solder 

The  entire  length  of  the  solder  base  was  restricted  from  any  motion  in  the  Y  direction. 

The  resultant  stresses  at  potential  failure  points  (Figure  4-46a)  are  shown  in  Table  4-4. 
The  stress  on  the  inner  surface  of  the  upper  J-lead  elbow  was  extrapolated  from  stresses 
provided  at  the  element  centers.  As  shown  in  Figure  4-46a,  only  a  crude  representation 
of  the  elbow  was  made  in  the  Nastran  FEM.  This  is  the  reason  for  the  low  stresses  calcu¬ 
lated  in  this  region  of  the  Nastran  FEM.  Stress  contours  for  the  Nastran  FEA  are  shown 
in  Figure  4-49  and  the  deflected  shape  is  shown  in  Figure  4-48a.  Similar  results  were 
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obtained  with  ABAQUS.  The  displacements  were  scaled  by  a  factor  of  1000  to  clearly 
show  the  deformations.  A  graphical  representation  of  the  boundary  conditions  is  included. 
An  important  characteristic  of  the  deformed  shape  is  the  elongation  of  the  J-lead  and 
solder  base. 

SCALE:  tOOOK 


5=  0.3I9E-04  S=  0.625E-04 

a)  NRSTRRN  MODEL  b)  PROBE  MODEL 


c)  PROBE  MODEL  <1)  PROBE 

IU/CHIP  CRRRIER  FULL  20  MODEL 


Figure  4-48.  Loading  and  Imposed  Displacements 


Figure  4-49.  (b)  Max  Shear  Stress  -  Nastran 
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TABLE  4-4.  J-LEAD  SOLDER  JOINT  STRESSES/STRAINS 


FEM  Code  Model 

Load  Condition 

Principal 

Stress 

Lead  Elbow 
(psi) 

Max.  Shear 
Stress 
Solder 
(psi) 

ABAQUS 

J-Lead/Solder 
Section  4. 8. 1.1 

PXTHERM:  Imposed 
Displacements  to 
Simulate  Chip  Carrier/ 
Board  Expansion 

2580 

98 

NASTRAN 

J-Lead/Solder 
Section  4. 8. 1.1 

PXTHERM:  Imposed 
Displacements  to 
Simulate  Chip  Carrier/ 
Board  Expansion 

750 

90 

PROBE 

J-Lead/Solder 
Section  4.8.1 .2 

PXTHERM:  Imposed 
Displacements  to 
Simulate  Chip  Carrier/ 
Board  Expansion 

1840 

94 

PROBE 

J-Lead/Solder 
W/Chip  Carrier 
Section  4. 8. 1.3 

PXTHERM:  Imposed 
Displacements  to 
Simulate  Chip  Carrier/ 
Board  Expansion 

1390 

43 

PROBE 

J-Lead/Solder 
W/Chip  Carrier 
Section  4. 8. 1.3 

PDTEMP:  Applied 
Temperature  & 

PXTHERM:  Imposed 
Displacements  to 
Simulate  Chip  Carrier/ 
Board  Expansion 

1240 

1260 

PROBE 

Full  2-D  Model 
(J-Lead/Solder 
W/Chip  Carrier 
&  PC  Board) 
Section  4.8.1 .4 

PDTEMP:  Applied 
Temperature 

1400 

1240 
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4.8.1. 2  J-Lead/Solder  Joint  Probe  Finite  Element  Analysis 

Figures  4-46b  through  4-46d  are  Probe  FEM’S  used  to  investigate  the  effects  of  using 
different  boundary  conditions  and  load  conditions  within  the  J-lead  and  solder  joint 
finite  element  analyses.  The  FEM  in  Figure  4-46b  was  analyzed  using  the  identical 
boundary  conditions  and  load  conditions  as  the  Nastran  FEM  shown  in  Figure  4-46a.  It 
contains  76  J-lead  elements  and  31  solder  elements.  A  comparison  of  the  stress  results 
presented  in  Table  4-4  shows  higher  stresses  calculated  with  Probe,  using  a  p-level  of  8, 
than  those  calculated  using  NASTRAN.  The  Nastran  analysis  used  linear  interpolation 
within  the  element  boundaries,  which  is  equivalent  to  a  p-level  of  1.  The  stresses  in 
Table  4-4  are  located  at  potential  crack  initiation  points,  shown  in  Figure  4-46. 


The  Probe  FEM  included  mesh  refinements  around  the  area  of  the  J-lead  and  solder 
intersections.  This  was  necessary  to  ensure  convergence  of  both  energy  and  stresses/ 
strains  in  the  area  of  interest.  In  general,  a  finite  element  analysis  ensures  only  displace¬ 
ment  continuity  across  element  boundaries.  The  Probe  FTLA  used  higher  order  interpolation 
functions  and  mesh  refinement  to  ensure  that  the  stresses/strains  converged  across  like 
element  boundaries  within  some  acceptable  tolerance.  In  the  Probe  analysis,  the  toler¬ 
ance  was  less  than  3%  in  stresses  across  J-lead  element  boundaries  and  solder  element 
boundaries.  The  stresses  across  J-lead  element  and  solder  element  boundaries  will  al¬ 
ways  be  discontinuous  because  of  the  differences  in  material  properties  at  the  interface. 
Stress  contours  for  the  Probe  FEA  are  shown  in  Figure  4-50. 


The  Nastran  FEA  provides  stresses  at  element  centers.  In  Nastran,  stresses  calculated 
at  the  same  grid  point  from  different  elements  will  not  be  the  same.  When  grid  point 
stresses  are  calculated  in  Nastran,  a  weighting  procedure  is  used  which  gives  an  average 
stress  at  the  grid  point.  The  grid  point  stress  averaging  procedure  can  produce  misleading 
results  near  the  J-lead/solder  material  interface  if  separate  averaging  procedures  are  not 
used  for  the  different  materials.  A  single  stress  averaging  process  across  different  material 
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Figure  4-50.  (a)  Max  Principal  Stress  -  Probe  FEA 
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Figure  4-50.  (b)  Max  Shear  Stress  -  Probe  FEA 
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boundaries  will  generate  a  single  value  at  the  grid  point  using  an  average  of  stresses  on  ei¬ 
ther  side  of  the  material  interface.  Other  finite  element  codes  use  similar  methods  to  extract 
grid  point  stress  results.  The  methods  of  extracting  grid  point  stresses  is  important  and  is 
one  of  the  reasons  for  the  discrepancies  between  the  Nastran  FEA  and  the  Probe  FEA. 

The  deflected  shape  of  the  Probe  TEM  is  shown  in  Figure  4-48b.  Included  is  a  graphical 
representation  of  the  boundary  conditions,  identical  to  the  Nastran  FEM  of  Figure  4-48a. 
The  deformed  shape  is  very  similar  to  the  deflected  shape  of  the  Nastran  FEM  of  Figure 
4-48a.  Again,  the  J-lead  and  solder  base  elongation  is  noticeable. 

4.8.1.3  J-lead/Solder  Joint  Probe  FEA  with  Chip  Carrier 

The  Probe  FEM  in  Figure  4-46c  is  similar  to  tne  Probe  FEM  in  Figure  4-46b  except 
for  the  addition  of  the  symmetric  portion  of  the  chip  carrier.  This  model  includes  83  J- 
lead  finite  elements,  31  solder  finite  elements  and  60  chip  finite  elements.  Only  half  of 
the  chip  carrier  is  modeled  because  of  symmetry.  The  addition  of  the  chip  carrier  al¬ 
lowed  symmetrical  boundary  conditions  to  be  used  with  the  model.  The  symmetrical  con¬ 
ditions  assume  all  motion  in  the  X  and  Y  direction  is  fixed  at  the  center  of  the  chip 
carrier.  In  the  previous  Nastran  and  Probe  analyses,  all  motion  in  the  Y  direction  was 
fixed  at  the  upper  flange  of  the  J-lead.  Two  load  conditions  were  applied  to  this  2-D 
Probe  FEM.  The  first  load  condition  (PXTHERM)  was  identical  to  that  used  in  the  pre¬ 
vious  Nastran  and  Probe  analysis.  The  second  load  condition  (PDTEMP)  included  the 
thermal  expansion  induced  load  between  the  J-lead  material  and  solder  material.  Load 
condition  1  and  2  represent  the  thermal  expansion  mismatch  between  the  chip  carrier  and 
PCB  in  the  plus  X  direction.  The  applied  temperature  differential  (AT  =  Tmax  -  Tmin) 
for  load  condition  2  (PDTEMP)  was  50  degrees  C.  The  temperature  differential  in  this 
case  effected  only  the  J-lead  and  sol  !er  mateiials.  The  thermal  expansion  of  the  chip 
carrier  and  PCB  were  included  in  the  PXTHERM  imposed  displacements.  Therefore,  the 
CTE  of  the  chip  carrier  was  set  equal  to  0.0  when  the  second  load  condition  was  applied. 
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The  differences  between  the  results  of  the  Probe  FEA  of  Section  4. 8. 1.2  and  this  Probe 
FEA  for  the  PXTHERM  load  condition  are  a  direct  result  of  the  relaxation  of  the  Y  direc¬ 
tion  support  of  the  upper  J-lead  flange.  Removing  the  Y  direction  support  allowed  the  chip 
carrier  to  move  down  towards  the  PCB  when  the  load  was  applied.  As  shown  in  Table  4-4, 
the  J-lead  upper  elbow  stress  was  reduced  from  1840  psi  to  1390  psi  when  the  Y  direction 
support  was  eliminated.  This  is  also  noticeable  in  stress  contours  around  the  upper  flange 
radius  as  shown  in  Figure  4-50  and  4-51.  The  deflected  shape  for  the  PXTHERM  load 
condition  is  shown  in  Figure  4-48c.  Graphical  representation  of  the  boundary  condition  at 
the  base  of  the  solder  is  included.  Again,  there  is  elongation  of  the  solder  base.  There  is  a 
slight  change  in  the  deformed  shape  of  the  upper  J-lead  flange  attributed  to  the  relaxation 
of  the  Y  direction  support.  The  deflected  shape  for  the  second  load  condition  is  similar  to 
that  shown  in  Figure  4-48d  with  the  exception  of  the  flat  solder  base. 
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Figure  4-51 .  (b)  Max  Shear  Stress  -  Probe  w/Chlp  FEA 
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Figure  4-51.  (d)  Max  Shear  Stress  -  Probe  w/Chlp  FEA 
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Figure  4-51 .  (e)  Max  Principal  Strain  -  Probe  w/CNp  FEA 
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4.8.1.4  J-lead/SoIder  Joint  Probe  FEA,  Full  2-D  Model 

The  Probe  FEM  in  Figure  4-46d  includes  the  J-lead/solder  joint,  the  chip  carrier  and 
the  PCB.  This  model  contains  68  J-lead  elements,  37  solder  elements,  13  chip  carrier 
elements,  and  76  PCB  elements.  Symmetrical  conditions  are  imposed  at  the  center  of  the 
chip  carrier  and  at  the  center  of  the  PCB.  Only  one  half  of  the  chip  carrier  and  PCB  are 
modeled  because  of  symmetry.  Motion  in  the  X  direction  is  fixed  and  motion  in  the  Y 
direction  is  free  along  the  line  of  symmetry  and  the  right  edge  of  the  PCB  was  con¬ 
strained  in  the  y  direction.  The  applied  temperature  differential  is  50  degrees  C  and  will 
be  called  the  PDTEMP  load  condition.  Previous  finite  element  analyses  considered  the 
differential  expansion  between  the  chip  carrier  and  PCB  in  the  plus  X  direction  only. 

The  results  in  Table  4-4  show  a  big  difference  in  the  maximum  shear  stress  between 
the  FEA  using  the  PDTEMP  load  conditions  and  the  FEA  using  the  PXTF1ERM  load 
conditions.  The  stress  contours  of  Figure  4-52  shows  the  concentration  of  stress  around 
the  area  where  cracks  are  known  to  propagate.  This  indicates  the  importance  of  including 
the  differential  expansion  between  all  materials  and  the  use  of  realistic  boundary  conditions. 
These  results  compare  well  with  the  second  load  condition  of  the  Probe  model  in  section 
4.8. 1.3,  although  the  maximum  shear  strain  for  the  Probe  FEM  of  section  4.8. 1.3  was 
slightly  higher  because  the  solder  base  remained  flat.  The  increase  in  the  J-lead  elbow 
stress  from  1240  psi  to  1400  psi  was  due  to  the  increased  deflections  of  the  PCB. 
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Figure  4-52.  (c)  Max  Principal  Strain  -  Probe  w/Thermal  Loads 
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Figure  4-52.  (d)  Max  Shear  Strain  -  Probe  w/Thermal  Loads 
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The  deflected  shape  is  shown  in  Figure  4-48d.  There  is  a  clear  difference  between  the 
deformed  shape  of  this  Probe  FEM  and  the  other  FEMs  shown  in  Figures  4-48  a  through 
c.  The  deformations  at  the  J-lead/solder  boundary  result  from  the  thermal  expansion  co¬ 
efficients  mismatch.  The  displacements  at  the  center  of  the  solder  base  and  at  the  upper 
J-lead  flange/chip  carrier  interface  are  listed  in  Table  4-5.  These  grid  locations  are 
shown  in  Figure  4-46.  The  displacements  in  the  X  direction  match  up  well  with  the  ex¬ 
pected  thermal  expansion  at  50  deg  C.  The  displacements  at  the  center  of  the  solder 
base  are  larger  than  those  used  in  the  previous  analysis  because  of  the  larger  temperature 
differential  applied  to  the  PCB.  The  previous  analysis  assumed  a  53  degree  C  tempera¬ 
ture  differential  within  the  chip  carrier  and  a  46  degree  C  temperature  differential  within 
the  PCB.  The  50  degree  C  temperature  represented  an  average  of  these  two  temperatures. 

TABLE  4-5.  THERMAL  DISPLACEMENTS 


Grid 

X 

Y 

sx 

8Y 

Calculated 
(ST=50  C) 

8X  =  8*T*L 

Error 

% 

8 

0.10CE-01 

0.117E  +  00 

0.1 4."  ^-03 

0.2842E-03 

0.1466E-03 

0.00 

9 

0.100E-01 

0.121E  +  00 

0.14bo_-03 

0.2854E-03 

0.1466E-03 

0.00 

381 

0.35GE-01 

0.000E  +  00 

0.2299E-03 

0.241 4E-03 

0.2294E-03 

0.21 

4.8.2  Leadless  Solder  Joint  FEAs 

The  analyses  of  leadless  solder  joints  \v>h  Probe  FEMs  was  accomplished  in  three  steps: 
one  with  the  solder  only,  one  with  the  solder  and  chip  carrier,  and  one  with  the  solder,  chip 
carrier,  PCB  and  temperature  change.  ABAQUS  was  used  to  model  the  solder  joint  only. 

4.8.2.1  Leadless  Solder  Joint  Probe  FEA 

The  leadless  solder  joint  geometry  is  shown  in  Figure  4-53a  and  the  material  proper¬ 
ties  are  shown  in  Table  4-3.  This  georr~try  was  generated  using  information  from  Refer¬ 
ence  4-6.  The  2-D  Probe  FEM  contains  90  solder  elements.  For  boundary  conditions  at 
the  vertical  support,  motion  is  fixed  in  the  X  direction  and  motion  is  free  in  the  Y  direc¬ 
tion.  For  boundary  conditions  at  the  horizontal  support,  motion  is  fixed  in  the  X  direc¬ 
tion  and  Y  direction.  The  applied  loading  represents  the  differential  thermal  expansion 
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Figure  4-53.  '.eadiess  Solder  Joint  Geometry 

in  the  plus  X  direction  between  the  chip  carrier  and  PCB.  This  plus  X  expansion  is  called 
the  PXTHERM  load  condition.  The  temperature  differential  (AT  =  Tmax  -  Tmin)  of 
the  chip  carrier  is  53  degrees  C  and  the  temperature  differential  of  the  PCB  is  46  degrees 
C.  The  loading  was  applied  to  the  finite  element  model  using  imposed  displacements  in 
the  plus  X  direction  along  the  base  of  the  solder  as  shown  in  Figure  4-54.  At  the  left 
corner  of  the  solder  base,  the  imposed  displacement  is  +  .424E-04  inches  in  the  X  direc¬ 
tion.  At  the  right  comer  of  the  solder  base,  the  imposed  displacement  is  +  .691E-04 
inches  in  the  X  direction.  These  displacements  were  calculated  using  the  following  formula: 

8x  =  (Cl  -  C3)  *  (x  +  C2),  when  -0.05  <  x  <  0.00  Eq.  28 
8x  =  Cl  *  (x  +  C2)  -  C3  *  C2,  when 

0.00  <  x  <  0.05  Eq.  29 

where:  Cl  —  4.37E-4  in/in,  PCB  thermal  strain 
C2  =  .483  in,  one  half  chip  carrier  length 
C3  =  3.392E-4  in/in,  chip  carrier  thermal  strain 

The  entire  length  of  the  solder  base  was  restricted  from  any  motion  in  the  Y  direction. 
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Figure  4-54.  Loading  and  Imposed  Displacements 


The  resultant  stresses  at  potential  failure  locations  (shown  in  Figure  4-53)  are  included 
in  Table  4-6.  As  expected,  the  sharp  notch  near  the  chip  carrier  corner  was  the  location 
of  the  peak  stress  although  there  were  aiso  high  stresses  at  the  left  end  of  the  solder  joint. 
The  stress  contours  for  the  2-D  Probe  FEA  are  shown  in  Figure  4-55.  The  deformed 
shape  shown  in  Figure  4-54a  indicates  shea/  deformation  along  the  solder  base  between 
the  chip  carrier  and  PCB.  The  displacements  were  scaled  by  a  factor  of  100  to  clearly 
show  the  deformations.  A  graphical  representation  of  the  boundary  conditions  used  in 
the  FEA  is  included. 
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TABLE  4-6.  LEADLESS  CHIP  CARRIER  SOLDER  JOINT 

Max.  Shear  Stress 


Model 

Load  Condition 

Location  1 : 

Near  the  Chip 
Carrier  Corner 
(PSI) 

Location  2: 

Upper  Left 

End 

(PSI) 

Location  3: 

Lower  Left 
End 
(PSI) 

Solder 
Joint 
Section 
4.8.2. 1 

PXTHERM:  Imposed 
Displacements  to 
Simulate  Chip  Carrier 
/Board  Expansion 

11100 

5970 

2170 

Solder 

Joint 

W /  Chip 

Carrier 

Section 

4. 8.2.2 

PXTHERM:  Imposed 
Displacements  to 
Simulate  Chip  Carrier 
/Board  Expansion 

3120 

1070 

1590 

Solder 

Joint 

W/  Chip 
Carrier 
Section 
4.8.23 

PDTEMP:  Applied 
Temperature  & 
PXTHERM:  Imposed 
Displacements  to 
Simulate  Chip  Carrier 
/Board  Expansion 

6100 

2300 

2460 

4-100 


CG SA  DATA  FILE: 


Figure  4-55.  (b)  Max  Shear  Stress  -  Probe  FEA 
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4.8.2.2  Leadless  Solder  Joint  Probe  FEA  with  Chip  Carrier 

The  ieadless  solder  joint  with  chip  carrier  geometry  is  shown  in  Figure  4-53b  and  the 
material  properties  are  shown  in  Thble  4-3.  The  geometry  was  generated  using  informa¬ 
tion  from  Reference  4-6.  The  2-D  Probe  FEM  contains  97  solder  elements  and  64  chip 
carrier  elements.  Symmetric  boundary  conditions  were  used  at  the  center  of  the  chip  car¬ 
rier.  Motion  in  the  X  direction  is  fixed  and  motion  in  the  Y  direction  is  free.  At  the  sol¬ 
der  base  motion  in  the  Y  direction  is  fixed  and  motion  in  the  X  direction  is  free.  The 
applied  loading  represents  the  differential  thermal  expansion  in  the  plus  X  direction  be¬ 
tween  the  chip  carrier  and  PCB  and  is  identical  to  that  used  previously  in  section  4. 8. 2.1. 
This  load  condition  is  called  PXTHERM. 

The  resultant  stresses  at  potential  failure  locations  are  shown  in  Table  4-6.  As  ex¬ 
pected,  the  sharp  notch  near  the  chip  carrier  corner  was  the  location  of  the  peak  stress. 
When  the  chip  carrier  was  added  to  the  leadless  solder  joint  FEM,  the  magnitude  of  the 
solder  stresses  decreased  because  of  chip  carrier  bending.  The  stress  contours  for  the 
2-D  Probe  FEA  are  shown  in  Figure  4-56. 
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Figure  4-56.  (a)  Max  Principal  Stress  -  Probe  w/Chlp  FEA 
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Figure  4-56.  (b)  Max  Shear  Stress  -  Probe  w/Chip  FEA 
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The  deflected  shape  of  the  Probe  FEM  the  chip  carrier  is  shown  in  Figure  4-54b.  In¬ 
cluded  is  a  graphical  representation  of  the  boundary  conditions  and  an  enlarged  plot  of 
the  solder  joint.  The  deformation  of  the  chip  carrier  is  clearly  shown.  This  deformation 
increased  the  compressive  loading  on  the  left  edge  of  the  solder  joint.  The  compression 
is  shown  in  the  stress  contour  plots  of  Figure  4-56. 

4.8.2.3  Leadless  Solder  Joint  Probe  FEA  with  Chip  carrier  and  Thermal  Loads 

The  leadless  solder  joint  and  chip  carrier  geometry  is  shown  in  Figure  4-53b.  The  ma¬ 
terial  properties  are  shown  in  Table  4-3  and  the  geometry  was  generated  using  informa¬ 
tion  from  Reference  4-6.  The  2-D  Probe  FEM  contains  97  solder  elements  and  64  chip 
carrier  elements.  Symmetric  boundary  conditions  were  used  at  the  center  of  the  chip  car¬ 
rier.  Motion  in  the  X  direction  is  fixed  and  motion  in  the  Y  direction  is  free.  The 
applied  loading  includes  imposed  displacements  which  represent  the  thermal  expansion  of 
the  PCB  in  the  plus  X  direction  for  a  temperature  differential  (AT  =  Tmax  -  Tmjn)  of  50 
degrees  C  and  an  applied  temperature  differential  of  50  degrees  C.  The  imposed  dis¬ 
placements  load  is  called  PXTHERM  and  the  thermal  load  is  called  PDTEMP. 

The  resultant  stresses  and  strains  for  the  2-D  Probe  FEA  with  thermal  loads  are  pres¬ 
ented  in  Table  4-6.  When  the  temperature  loads  are  applied,  radically  different  results 
are  obtained  because  of  the  thermal  expansion  coefficient  at  the  chip  carrier/solder  inter¬ 
face.  This  difference  is  shown  in  the  stress  contour  plots  presented  in  Figure  4-57. 

The  deflected  snap*  oi  u»  s  t-EM  A  shown  iu  i;g-*e  4-54c.  A  graphical  representation 
of  the  boundary  conditions  and  an  enlarged  plot  of  the  solder  joint  is  included.  The  de¬ 
formation  of  the  chip  carrier  is  clearly  shown.  The  increased  bending  of  the  chip  carrier 
was  the  result  of  using  an  average  temperature  differential  of  50  degrees  C.  The  imposed 
displacements  used  in  the  previous  leadless  solder  joint  FEA  assumed  a  53  degree  C  tem- 
peratur  ',<fferential  within  the  chip  carrier  and  a  46  degree  temperature  differential  with¬ 
in  the  PCB. 
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Figure  4-57.  (a)  Max  Principal  Stress  -  Probe  w/Thermal  Loads 
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Figure  4-57.  (c)  Max  Principal  Strain  -  Probe  w/Thermal  Loads 
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Figure  4-57.  (d)  Max  Shear  Strain  -  Probe  w/Thermal  Loads 
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4.8.3  Correlation  of  FEA  Results  and  Known  Failure  Locations 

Solder  joint  cracks  in  the  J-lead/solder  joint  originate  either  near  the  outside 
J-lead/solder  intersection  within  the  solder  material  or  from  voids  within  the  solder 
material  near  the  J-lead  interface.  The  cracks  propagate  around  and  near  the  J-lead 
interface  within  the  solder  material.  In  the  models  discussed  above,  the  peak  stresses  do 
not  occur  near  the  outer  surface  of  the  solder  at  the  J-lead/solder  intersection  (See  Fig¬ 
ures  4-49  through  Figure  4-52).  This  indicates  that  either  voids  in  the  solder  material 
near  the  J-lead  interface  or  surface  flaws  near  the  J-lead/solder  interface  play  a  part  in 
the  crack  initiation  process.  A  circular  shaped  void  will  magnify  the  through  stress  by  a 
factor  of  2  to  4,  depending  upon  the  stress  distribution  around  the  void.  For  examples  of 
J-lead/solder  joint  cracking,  see  References  4-8  and  4-9.  In  the  2-D  Probe  J-lead/sold- 
er  FEM  that  allowed  thermal  expansion  in  both  X  and  V  directions  (Figure  4-52),  the 
maximum  total  principal  strains  and  maximum  shear  strains  were  located  in  areas  where 
cracking  occurs.  However,  as  previously  shown,  the  strains  at  the  outside  intersection  of 
the  J-lead/solder  materials  were  lower  than  the  maximum  principal  strains.  None  of  the 
FEMs  included  flaws  on  the  solder  surface  or  interior  to  the  solder  material. 

Solder  joint  cracks  in  leadless  joints  have  been  shown  to  originate  either  near  the  in¬ 
side  solder/chip  carrier  interface  within  the  solder  material  or  from  voids  within  the  sol¬ 
der  material  (see  References  4-7  and  4-12  for  examples  of  leadless  soider  joint  cracking). 
These  reports  show  the  crack  propagating  from  the  inside  solder/chip  carrier  interface, 
along  the  interface  until  reaching  the  corner  and  then  proceeds  normal  to  the  outside  sol- 
dei  surface  until  fracture  occurs.  This  process  is  fully  described  in  Reference  4-7.  The 
peak  stresses/strains  will  occur  at  the  corner  notch  as  shown  in  Figure  4-55  through 
Figure  4-57  and  there  is  the  possibility  that  when  subjected  to  thermal  loads,  cracking 
may  originate  near  this  corner  notch. 
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Chapter  5 

RELIABILITY  PREDICTIONS 


5.0  Introduction 

The  reliability  of  a  comporent  is  described  as  its  ability  to  function  without  failure. 

This  chapter  addresses  the  theories  and  procedures  needed  to  analyze  the  failure  mecha¬ 
nisms  covered  in  Chapter  3.  These  failure  theories  can  be  used  to  determine  if  fracture 
of  the  material  will  occur  under  operational  conditions.  This  chapter  also  covers  proce¬ 
dures  for  predicting  fatigue  life  and  creep  rupture  time  to  failure. 

5.1  Deformation 

Finite  element  methods  for  computing  deformation  of  electronics  during  vibration  and 
temperature  changes  were  covered  in  Chapter  4.  Deformation  results  in  stresses  which 
can  exceed  the  strength  of  the  material  and  cause  short  fatigue  lives.  This  chapter  focuses 
on  using  these  calculated  stresses  to  predict  the  reliability  of  electronics. 

5.2  Fatigue 

5.2.1  Finite  Element  Stresses  for  Fatigue  Analysis 

The  fatigue  life  of  components  subjected  to  bending  or  axial  loads  is  dictated  by  the 
principal  stress,  Oi,  obtained  from  FEA  This  principal  stress  represents  the  maximum 
normal  stress  acting  on  an  element.  In  general,  cracks  propagate  in  a  direction  perpen¬ 
dicular  to  the  direction  of  the  principal  stress  (Figure  5-1).  The  principal  stresses  should 
be  used  in  predicting  the  fatigue  life  of  leads  connecting  the  components  to  the  circuit 
board.  Once  vibration  and  thermal  stresses  are  determined  from  the  finite  element  anal¬ 
ysis,  a  curve  of  axial  stress  amplitude  (Ao/2)  vs.  cycles  to  failure  (N)  for  the  appropriate 
material  should  be  used  to  predict  fatigue  life.  Stress  vs.  life  curves  for  various  materials 
are  included  in  Chapter  2.  A  different  characteristic  stress  is  needed  for  conditions  in 
which  the  cracks  are  forced  to  propagate  under  shear  deformation.  This  type  of  crack  is 
found,  for  example,  along  the  boundary  between  a  solder  joint  and  a  leadless  chip  carrier 
(LCC).  As  the  board  expands  during  an  increase  in  temperature,  the  joint  is  subjected  to 
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a  shear  deformation  as  illustrated  in  Figure  5-2.  The  maximum  shear  stress,  rmax,  should 
be  used  in  this  case  for  the  prediction  of  fatigue  life.  If  the  component  analyzed  is 
sheared  primarily  in  one  direction,  the  maximum  shear  stress  is  one  of  the  FEA  output 
shear  stresses:  Txy,  or  ryz.  This  occurs  in  the  case  of  the  LCC  solder  joint  under  tem¬ 
perature  changes.  If  a  more  complex  load  condition  occurs,  the  maximum  shear  stress 
can  be  computed  from  the  principal  stresses  oi,  02,  and  03.  The  maximum  shear  stress, 
t max'  is  the  greater  of: 


[  er,  -  cr2 1  |o,-a,j  |a2-a3|  , 

2  Eq‘  1 

I  he  fatigue  life  is  then  obtained  from  a  plot  of  shear  stress  amplitude  (At/2)  versus 
cycles  to  failure  (N).  Under  multiaxial  load  conditions,  the  use  of  a,  or  Tmax  to  predict 
fatigue  life  may  be  invalid.  In  these  cases,  the  use  of  the  Von  Mises  stress  is  recommended 
because  it  contains  the  three  principal  stresses.  The  Von  Mises  stress  computed  by  FEA 
can  then  be  used  to  predict  life  with  material  data  in  terms  of  Von  Mises  stress  vs.  N. 
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Expansion  Differential  J 
Under  Temperature  Increase  — 


GP03-0624-29-DSUZ 


Figure  5-2.  Shear  Deformation 


Fatigue  data  can  also  be  presented  in  terms  of  strain  amplitude  (Ae/2)  vs.  cycles  (N)  to 
failure.  This  format  is  typically  used  when  cyclic  stresses  are  high  enough  to  cause  yield¬ 
ing  and  plastic  deformation  in  the  critical  area.  The  resulting  condition  is  known  as  “low 
cycle  fatigue”  because  failure  under  high  stress  occurs  in  relatively  few  cycles  (less  than 
105  cycles).  Strain-life  data  is  obtained  by  controlling  the  strain  as  the  load  is  applied. 
The  applied  strain  is  repeated  in  each  cycle  until  failure  of  the  test  specimen  occurs.  In 
low  cycle  fatigue,  strain  is  the  controlling  parameter  because  the  material  is  in  the  plastic 
region  of  the  stress-strain  curve.  In  this  region,  small  deviations  in  stress  can  result  in 
large  changes  in  strain.  Therefore,  it  is  difficult  to  reach  a  desired  point  in  the  plastic  re¬ 
gion  of  the  stress-strain  curve  without  monitoring  strain. 

The  cycles  to  failure  can  be  predicted  from  strain-life  data  and  the  strains  computed  in 
the  finite  element  analysis.  Life  predictions  for  parts  under  axial  or  bending  loads  should  be 
based  on  the  principal  strain,  €1.  The  life  of  materials  undergoing  shear  deformation  should 
be  predicted  by  using  the  maximum  shear  strain.  ymax.  In  a  two  dimensional  analysis,  the 
maximum  shear  strain  can  be  determined  from  the  principal  strains  ei  and  e?- 

ymax  =  O  -  e2  Eq.  2 

When  using  strains  from  finite  element  analyses,  the  correct  strain  needs  to  be  used  to 
predict  fatigue  life.  In  Chapter  4,  stresses  and  strains  in  a  J-lead  exposed  to  temperature 
changes  were  obtained  with  PROBE  finite  element  analyses.  In  several  of  the  analyses, 
deflections  were  imposed  on  the  J-lead  finite  element  model  to  simulate  thermal 
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expansion  of  the  circuit  board  and  the  chip  carrier.  The  principal  strains  obtained  in 
these  analyses  can  be  used  directly  with  strain  amplitude  vs.  life  data  to  predict  the  fa¬ 
tigue  life  of  the  material. 

In  another  analysis  of  the  J-lead,  the  temperature  change  was  imposed  on  the  chip 
carrier-lead-circuit  board  model.  This  allowed  the  FEA  to  determine  the  thermal  expan¬ 
sion  of  the  materials  without  externally  imposed  displacements.  The  principal  strain  out¬ 
put  from  this  PROBE  example  consists  of  the  mechanical  strain  and  the  thermal  strain: 

£l  =  ^mech  "h  ^thermal  EQ-  3 

The  thermal  strain  is  the  product  of  the  material’s  coefficient  of  thermal  expansion 
(CTE)  and  the  temperature  differential: 


^thermal  =  CTE  X  AT  Eq.  4 

The  mechanical  strain  is  caused  by  the  interaction  of  adjacent  materials  having  differ¬ 
ent  thermal  expansion  coefficients.  Materials  with  lower  CTE  prevent  the  expansion  of 
materials  with  higher  CTE.  This  results  in  compressive  mechanical  strains  in  the  high 
CTE  materials  and  tensile  mechanical  strains  in  the  low  CTE  materials.  Mechanical 
strains  cause  stresses  which  lead  to  fatigue  damage.  Therefore,  when  predicting  fatigue 
life,  only  the  mechanical  component  of  ei  should  be  used.  The  mechanical  strain  can  be 
extracted  from  €1  a'  follows: 


fmech  -  €i-  (CTE  x  AT)  Eq.  5 

This  mechanical  strain  is  then  used  with  strain  amplitude  vs.  life  curves  to  predict  the 
number  of  cycles  to  failure.  Fatigue  life  can  also  be  determined  by  using  the  FEA  princi¬ 
pal  stress,  Oi,  and  stress  vs.  life  data.  Both  procedures  will  result  in  equivalent  fatigue  life 
predictions.  Table  5-1  summarizes  the  application  of  the  various  FEA  stresses  and  strains 
in  fatigue  analysis. 
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TABLE  5-1.  APPLICATION  OF  FEA  OUTPUT  IN  FATIGUE  LIFE  PREDICTION 


Parameter 

Symbol 

Application 

Fatigue  Data  Needed 
for  Life  Prediction 

Principal 

Stress 

cr, 

Materials  under  axial  or 
bending  loads 

Act/2  vs.  N 

Max.  Shear 
Stress 

Tmax 

Materials  under  constrained 
shear  deformation 

At/2  vs.  N 

Principal 

Strain 

e, 

Materials  under  axial  or 
bending  loads 

A  e/2  vs.  N 
(low  cycle  fatigue) 

Mechanical 

Principal 

Strain 

^mech 

Materials  exposed  to 
temperature  changes  and 
analyzed  with  FEA  under 
applied  AT 

A  e/2  vs.  N 
(low  cycle  fatigue) 

Max.  Shear 

Strain 

y  max 

Materials  under  constrained 
shear  deformation 

Ay/2  vs.  N 
(low  cycle  fatigue) 

Stress  or  strain  output  from  finite  element  analyses  can  be  used  in  fatigue  life  predic¬ 
tion  with  equal  accuracy.  The  choice  depends  on  whether  stress-life  or  strain-life  curves 
are  available  for  the  material.  Even  in  the  case  where  only  stress-life  curves  are  avail¬ 
able,  these  data  can  be  converted  to  strain  vs.  life  by  using  the  stress-strain  curve  for  the 
material. 

Table  5-1  assumes  that  the  applied  stresses  (oi  or  Tmax)  and  strains  (eL  or  ymax)  are  ful¬ 
ly  reversed  during  the  cycle.  For  a  material  under  axial  loads,  the  stress  at  a  point  should 
vary  from  Oi  to  -oq  during  each  cycle.  This  condition  occurs  when  the  temperature  in¬ 
crease  above  a  baseline  condition  is  balanced  by  a  temperature  drop  of  the  same  magni¬ 
tude  below  the  baseline.  As  an  example,  consider  a  temperature  cycle  defined  by: 

Maximum  temperature  =  100°C 
Minimum  temperature  =  -50°C 
Baseline  temperature  =  25°C  (Room  temperature) 

This  temperature  cycle  will  result  in  fully  reversed  stresses  because  the  75°C  increase 
above  the  baseline  equals  the  75°C  drop  below  the  room  temperature  condition. 
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In  the  case  of  cycles  where  the  stresses  are  not  fully  reversed,  the  above  procedures  need 
to  be  modified.  Methods  for  analyzing  mean  stress  effects  and  combinations  of  vibration 
and  thermal  stresses  are  covered  later  in  this  chapter.  Load  cycles  where  stresses  exceed  the 
yield  strength  of  the  material  and  plastic  deformation  occurs  are  also  discussed. 

5.2.2  Fatigue  Life  Prediction 

Once  the  appropriate  stress  information  is  obtained  from  the  finite  element  analysis, 
the  number  of  stress  cycles  which  can  be  sustained  by  a  material  can  be  determined  with 
a  fatigue  analysis.  If  the  same  stress  level  is  attained  in  each  cycle,  the  stress  profile  is 
known  as  a  constant  amplitude  stress  history.  The  number  of  cycles  to  failure  can  then  be 
determined  from  a  material  curve  of  stress  versus  cycles  as  described  in  Section  3.3.  If 
the  stress  levels  vary  from  cycle  to  cycle,  the  stress  profile  is  a  variable  amplitude  stress 
history.  Variations  in  stress  levels  can  be  caused  by  variations  in  maximum  temperature, 
and  random  vibration,  for  example.  The  fatigue  damage  from  the  various  stress  cycles 
must  be  accounted  for  when  predicting  the  fatigue  life  of  the  material. 

The  first  step  in  the  analysis  requires  grouping  stress  cycles  with  similar  amplitudes 
together.  This  results  in  a  series  of  constant  amplitude  blocks  which  represent  the  stress 
history.  The  damage  caused  by  each  block  of  cycles  is  defined  by: 


where:  ns  =  number  of  cycles  at  Oj 

N,  =  number  of  cycles  to  failure  at  Oj 
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This  is  illustrated  in  Figure  5-3.  Once  the  damage  from  each  block  is  determined. 
Miner’s  cumulative  damage  formulation  (Reference  5-1)  can  be  used  to  sum  the  damage 
and  determine  if  failure  will  occur.  According  to  Miner,  failure  occurs  if: 

2—  >  1.0  Eq.  7 

Nj 

Reference  5-2  recommends  that  0.7  should  be  used  for  electronics  instead  of  1.0  in 
Equation  7.  This  results  in  more  conservative  life  predictions. 


GP03-0624-30-D/qg 

Figure  5-3.  Miner’s  Rule  Applied  to  Counted  Effective  Strain  Ranges 
5.2.3  Fatigue  Life  Predictions  Under  Combined  Stresses 

High  frequency  vibration  stresses  can  occur  in  combination  with  slower  stress  cycles 
caused  by  temperature  changes  or  high  G  aircraft  maneuver  loads  (Figure  5-4).  The  first 
step  in  a  fatigue  life  analysis  is  to  determine  the  stresses  caused  by  these  conditions.  A 
particular  temperature,  or  an  RMS  vibration  level,  is  selected  as  the  reference  condition. 
The  stresses  in  the  electronic  material  are  obtained  with  a  finite  element  analysis  for  the 
reference  condition.  Each  time  the  component  is  exposed  to  this  reference  condition  rep¬ 
resents  a  “load  cycle”.  The  magnitude  of  the  load  cycles  can  vary  depending  on  the  maxi¬ 
mum  temperature  of  each  cycle  or  the  level  of  vibration.  Higher  temperatures  cause 
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Temperature 


Thermal 


o 


I — 1 

Time 


Figure  5-4.  Slow  Thermal  Stresses  and  High  Frequency  Vibration  Stresses 


higher  thermal  stresses,  and  higher  vibration  accelerations  lead  to  greater  dynamic 
stresses.  In  general,  stress  is  proportional  to  temperature  and  vibration  level  as  long  as 
the  stresses  are  elastic  (below  the  yield  strength  of  the  material).  If  stresses  exceed  the 
yield  strength,  then  the  plastic  behavior  of  the  material  needs  to  be  considered.  The  plas¬ 
tic  stress  analysis  was  described  in  Section  4.3.4,  and  models  for  analyzing  cyclic  stresses 
in  the  plastic  range  are  covered  in  Section  5.2.4.  Once  the  stresses  are  determined  for 
the  different  conditions,  they  need  to  be  combined  into  a  load  spectrum  representing  the 
load  history  of  the  component.  The  fatigue  damage  caused  by  each  load  cycle  can  then 
be  summed  to  predict  how  many  cycles  can  be  experienced  by  the  material  until  failure 
occurs. 


5.2,3. 1  Independent  Vibration  and  Thermal  Stresses 

Steinberg  (Reference  5-2)  has  proposed  one  of  the  simplest  models  for  combining  vi¬ 
bration  and  thermal  stresses.  In  this  model,  vibration  is  assumed  to  be  independent  from 
thermal  stresses.  The  fatigue  damage  due  to  vibration  and  thermal  stress  cycles  are  de¬ 
fined  as  follows: 


Vibration  Fatigue  Damage  = 


_n v 
Nv 


Eq.  8 
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and 

Thermal  Fatique  Damage  =  Eq.  9 

where:  nv  =  number  of  vibration  cycles  in  a  year 

Nv  =  number  of  vibration  cycles  which  will  cause  failure 

(obtained  from  stress  vs.  life  S-N  data  for  the  material) 

nt  —  number  of  thermal  cycles  in  a  year 

Nt  =  number  of  thermal  cycles  which  will  cause  failure  (obtained  from  S-N 
curve  for  the  material) 

The  values  of  Nv  and  Nt  are  obtained  from  material  S-N  curves  similar  to  Figure  5-5. 
Constant  amplitude  stresses  on  these  curves  are  typically  “fully  reversed  stresses”  in 
which  the  "R”  ratio  of  minimum/maximum  stress  equals  -1.  Random  vibration  stresses 
do  not  contain  mean  stress  or  preload.  To  predict  the  number  of  years  of  operation  (Ny). 
the  vibration  and  thermal  fatigue  damage  are  then  combined  as  follows: 


Eq.  10 


Vibration  Cycles  Low  Frequency  Cycles 


Cycles  to  Failure  Cycles  to  Failure 

GPC3-0624-32  D 


Figure  5-5.  Fatigue  Data  Obtained  Under  Fully  Reversed  Stresses 
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Solving  for  Ny  gives  the  fatigue  life.  Although  this  model  will  give  an  approximate  fa¬ 
tigue  life  for  the  material,  it  does  not  account  for  simultaneous  action  of  the  vibration 
and  thermal  cycles  illustrated  in  Figure  5-6. 


Maneuver 

Stress 


Combined 

Stress 


Vibration 

Stress 


GP03-0624-33-DSUZ 


Figure  5-6.  Combined  Manuever  and  Vibration  Stresses 

Without  Relaxation  or  Creep 


5.2.3.2  Generating  Effective  Strain  vs.  Life  Curves 

When  vibration  and  temperature  cycling  occur  simultaneously,  the  thermal  stress  acts 
as  a  static  mean  stress  with  vibration  stresses  superimposed  (Figure  5-6).  The  determina¬ 
tion  of  vibration  fatigue  damage  under  this  condition  is  complicated  by  the  limitations  of 
the  available  S-N  data.  Most  S-N  curves  for  electronic  materials  have  been  generated 
for  tests  without  mean  stress.  These  data  need  to  be  modified  for  use  with  actual  load 
histories  which  contain  thermal  preloads.  This  data  modification  is  accomplished  by  first 
selecting  data  points  from  a  stress  vs.  life  curve  for  the  particular  material  (Figure  5-7). 
The  stress  values  (omax)  read  from  the  curve  are  used  to  estimate  the  corresponding  val¬ 
ues  of  strain  amplitude  (Ae/2): 


Ae  _  Omax 
2  ~  E 


Eq.  11 
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Cycles  to  Failure 

G  PC 3-0624- 34-D/cjg 


Figure  5-7.  Assembly  of  Effective  Strain  vs  Life  Curve 

where  E  is  the  material  modulus  of  elasticity.  Equation  11  is  not  valid  if  the  maximum 
stress  exceeds  the  yield  strength  of  the  material.  If  yielding  occurs,  the  strains  (Ag/2) 
need  to  be  estimated  from  the  cyclic  stress  vs.  strain  curve  for  the  material  as  illustrated 
in  Figure  5-8.  The  maximum  stress  (omaxX  the  strain  (Ae/2),  and  E  are  then  substituted 
into  the  'following  effective  strain  equation  developed  by  Smith,  Watson  and  Topper 
(Reference  5-3): 


Eq.12 


An  effective  strain  is  computed  for  each  point  selected  in  the  original  o-N  curve  to 
form  a  new  eeff  vs.  life  curve.  In  the  final  step  in  the  fatigue  life  prediction  procedure,  the 
stresses  listed  in  the  stress  history  are  converted  into  effective  strains.  The  fatigue  dam¬ 
age  caused  by  each  stress  cycle  can  then  be  determined  from  the  eeff  vs.  life  curve. 
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Figure  5-8.  The  Strain  Amplitude  Ae/2  Can  Be  Obtained 
From  the  Stress-Strain  Curve 


5.2J.3  Generating  Effective  Strain  Histories 

Vibration  load  histories  are  composed  of  high  frequency  cycles  with  random  ampli¬ 
tude.  Reference  5-2  suggests  that  random  amplitude  cycles  can  be  modeled  with  a  Gaus¬ 
sian  distribution  of  maximum  and  minimum  stress  levels  (Figure  5-9).  The  amplitudes  in 
a  Gaussian  distribution  are  arranged  as  follows: 


68.3%  of  cycles  are  between  0  and  1  RMS 
27.1%  of  cycles  are  between  1  and  2  RMS 
4.3%  of  cycles  are  between  2  and  3  RMS 


G  P03-0624-36-  D/dp! 


Figure  5-9.  Gaussian  Distribution  of  Stress  Amplitudes 
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A  Rayleigh  distribution  an  also  be  used  during  vibration  analysis  (Figure  5-10)  to  de¬ 
termine  the  peak  stress  distributions  (Reference  5-4).  The  stress  peaks  in  a  Rayleigh  dis¬ 
tribution  are  broken  down  as  follows: 

39.3%  of  peaks  are  between  0  and  1  RMS 
47.2%  of  peaks  are  between  1  and  2  RMS 
12.3%  of  peaks  a*e  between  2  and  3  RMS 
1.2%  of  peaks  are  between  3  and  4  RMS 

Rayleigh  Distribution 


Figure  5-10.  Rayleigh  Distribution  of  Stress  Peaks 


With  these  assumed  distributions,  the  random  vibration  cycles  can  be  represented  by  a 
much  simpler  mode!  composed  of  four  blocks  of  constant  amplitude  cycles  (Figure  5-11). 
The  maximum  and  minimum  stresses  are  defined  ior  each  block  as: 


(^max)!  — 
(0)nin)l  — 

0)  herbal  +  CTrms 
^thermal  “  <?RMS 

Eq.  13 

(0)nax^2 

^thermal  -  Orms 

Eq.  14 

>  'rruii)2  = 

^thermal  —  -  Or  ms 

(Omax)3  = 
(0)71111)3  = 

Otnermal  +  3  *  O^MS 

thermal  "  3  *  CTrm.S 

Eq.  15 

(ormax)4  — 

Orthermal  +  4  *  <7rms 

Fq.  IP 

(0)1110)4  ~ ■ 

O'thcrmal  “  4  *  °rRMS 
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3*0  RMS 


Figure  5-11.  Model  of  Random  Vibration  With  Blocks  of  Constant  Amplitude  Cycles 

where:  othermai  =  thermal  mean  stress 

oRMS  =  root-mean-square  of  the  vibration  stress 
The  strain  amplitude  of  each  block  can  be  determined  with  the  following  procedure: 


(^max)i 
(fmin)i  — 


(tTrnax)i 

E 

(°min)i 


Eq.  17 
Eq.  18 


where  i  -  1,  2.  3,  or  4  RMS  levels.  The  strain  amplitude  is  then: 

(Af)i  ==  (fmax)i  —  (fmin)i  ^ 


The  effective  strain  characterizing  each  block  of  constant  amplitude  cycles  is  computed 
with  Equation  12: 


yfoVnax)' 

E  2 


(Ae)j 


Eq.  20 


Finally,  the  fatigue  damage  due  to  each  block  of  cycles  is: 

Fatigue  Damage,  =  — r  Eq.  21 

^  I 

where:  n,  =  number  of  vibration  cycles  of  amplitude  i*RMS  which  occur  in  a  year 

N,  =  number  of  cycles  of  amplitude  i*RMS  which  will  cause  failure 
(obtained  from  e,.tf  vs.  N  curve) 
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As  an  example,  consider  a  component  mounted  on  a  PCB  which  will  vibrate  100  hours 
in  a  year.  If  the  first  mode  natural  frequency  (of  the  PCB)  is  200  Hz,  the  number  of 
cycles  with  1*Orms  amplitude  is  defined  by  the  Rayleigh  distribution: 

n  [  =  0.393  *  100  Hrs  *  3600  sec/Hr  *  200  cycles/sec 
=  2.83  *  107  cycles  in  a  year 


The  value  of  Ni  is  obtained  from  the  eeff  vs.  life  curve  at  an  effective  strain  (eeff)i  as 
illustrated  in  Figure  5-12.  Similar  computations  can  be  repeated  for  the  other  blocks  of 
amplitude  2*Orms>  3*cjrms,  and  4*orms-  To  predict  the  number  of  years  of  operation 
(Nv),  the  damage  due  to  each  block  is  added: 


Solving  for  Ny  gives  the  fatigue  life  of  the  material  in  terms  of  years  of  operation. 


Effective 

Strain 


Figure  5-12.  Effective  Strain  Is  Used  to  Predict  Fatique  Life 
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The  way  in  which  the  stresses  are  combined  depends  on  the  material  and  load  condi¬ 
tion.  When  the  analysis  focuses  on  materials  which  do  not  creep  the  thermal  stresses  are 
sustained  during  the  entire  temperature  cycle.  The  procedures  described  above  are  then 
valid. 


5.2.4  Cyclic  Plastic  Stresses 

As  long  as  the  applied  stresses  do  not  exceed  the  material  yield  strength,  the  cyclic 
stress-strain  behavior  will  remain  linear  (Figure  5-13).  However,  if  the  stress  exceeds  the 
yield  strength,  plastic  deformation  occurs,  and  the  material  will  follow  its  stress-strain 
curve  to  the  peak  of  the  cycle  (Figure  5-14). 


Stress 


Stress 


Figure  5-13.  Linear  Elastic  Cyclic  Figure  5-14.  Plastic  Stress-Strain 

Stress-Strain  Behavior  Condition  at  the  Peak  of  the  Cycle 

Section  4.3.4  described  procedures  for  estimating  plastic  stresses  from  linear  elastic 
finite  element  analyses.  One  procedure  assumes  that  the  actual  plastic  strain  equals  the 
strain  computed  with  finite  element  analysis.  The  plastic  stress  is  read  from  the  stress- 
strain  curve  at  the  given  strain.  A  second  more  conservative  estimate  can  be  obtained  by 
using  Neuber's  rule  which  assumes  that: 


(oinax^  max)actual  —  (ormax^max)FEA 


Eq.  24 
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The  actual  plastic  stress  and  strain  combination  must  fall  on  the  curve  of  a  vs.  e. 

These  procedures  are  used  to  find  the  stress  and  strain  at  the  peak  of  the  load  cycle.  A 
similar  procedure  can  be  used  to  determine  the  stress  and  strain  state  under  reversed 
loads.  Just  as  the  material  follows  its  stress-strain  curve  during  loading,  the  material  fol¬ 
lows  a  set  path  in  the  unloading  part  of  the  cycle  (Figure  5-15).  This  path,  known  as  th^ 
hysteresis  curve,  describes  the  cyclic  stress-strain  behavior  in  the  plastic  range  (Figure 
5-16).  The  hysteresis  curve  for  a  material  can  be  determined  from  cyclic  stress-strain 
characterization  tests.  If  data  is  not  available,  the  hysteresis  curve  can  be  approximated 


Figure  5-15.  Stress-Strain  Behavior  Figure  5-16.  Hysteresis  Stress-Strain 

During  a  Reversed  Load  Curves  Under  Fully  Reversed  Cycles 

trom  the  stress-strain  curve.  One  approach  is  to  assume  that  the  hysteresis  curve  is  simi¬ 
lar  to  the  stress-strain  curve,  but  twice  as  large  (Reference  5-5).  To  generate  the  curve,  a 
and  e  values  of  selected  data  points  on  the  .>tress-strain  curve  are  multiplied  by  two: 


ahys  =  2  x  a 

fhys 


Eq.  25 


1  hese  data  points  then  describe  the  hysteresis  curve  as  illustrated  in  Figure  5-17. 

Once  the  hysteresis  curve  is  established,  it  can  be  used  to  determine  the  change  in 
stress  (Ao)  and  strain  (Ae)  during  the  cycle  illustrated  in  Figure  5-18.  In  this  figure,  the 
origin  of  the  “inverted”  hysteresis  curve  is  at  the  point  of  load  reversal.  If  plasticity  oc- 
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curs  during  the  unloading  part  of  the  cycle,  Neuber’s  rule  can  be  used  once  again  to  ob¬ 
tain  the  actual  state  of  stress  and  strain: 


(A a  x  Atr)actuai  =  (A a  x  A6)fea 


Eq.  26 


where:  AopEA  —  (<7max  ~  £Fnin)FEA 


Eq.  27 


Acfea  =  Actfea  /  E 


GPC3-0624-49-D/qg 


Figure  5-17.  Hysteresis  Curve  Generated 
From  the  Material  Stress-Strain  Curve 


Figure  5-18.  Cyclic  Stress-Strain 
Range 


The  actual  change  in  stress  Ao  and  the  change  in  strain  Ae  are  obtained  from  the  hys¬ 
teresis  curve  at  the  point  of  the  curve  which  satisfies  Equation  26: 


ACTactual  —  ^hys 


Eq.  29 


Afactual  —  fhys  Eq.  30 

Subsequent  constant  amplitude  load  cycles  will  reach  the  actual  maximum  stress  and 
strain  aud  will  have  the  same  stress  and  strain  range  as  illustrated  in  Figure  5-19. 
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GP03-0624-50-D/qg 

Figure  5-19.  Stress-Strain  Loop  During  Constant 
Amplitude  Cycles 

5.3  Creep  and  Stress  Relaxation 

5.3.1  Creep  Failure  Under  Constant  Load 

Solder  is  an  example  of  a  material  which  undergoes  time  dependent  deformation,  or 
creep,  under  sustained  stress.  This  type  of  material  will  continue  to  deform  until  rupture 
occurs.  In  general,  the  higher  the  stress  the  shorter  is  the  time  to  failure.  Figure  5-20 
illustrates  creep  rupture  data  for  60-40  solder  at  various  temperatures  (Reference  5-6). 
To  generate  the  data  in  Figure  5-20,  the  stress  in  each  test  was  held  constant  until  failure 
of  the  specimen  occurred. 

After  completing  a  finite  element  analysis,  the  time  to  failure  can  be  determined  from 
Figure  5-20  at  the  calculated  stress  acting  in  the  material.  This  type  of  analysis  is  appli¬ 
cable  to  cases  where  the  stress  level  is  maintained  continuously.  This  includes  electronics 
which  remain  on  for  significant  periods  of  time. 
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TIME  (HOURS) 

A'023-CJ 

Figure  5-20.  Creep  Rupture  Stress  of  Solder  (60SN  -  40  Pb) 


5.3.2  Combination  of  Creep/Stress  Relaxation  and  Fatigue 

When  predicting  the  fatigue  life  of  solder,  the  load  history  model  must  account  for 
creep  and  stress  relaxation  (Figure  5-21).  After  power-on,  the  heat  generated  by  the 
component  increases  the  temperature  of  the  case.  As  the  heat  dissipates  into  the  board, 
board  temperatures  increase,  but  at  a  lower  rate  than  the  case  temperature  increase  (Ref¬ 
erence  5-7).  The  chip  carrier,  therefore,  expands  first  as  illustrated  by  Condition  2  in 
Figure  5-21.  This  results  in  solder  thermal  stresses.  As  the  maximum  temperature  of  the 
case  is  reached,  the  board  catches  up  in  temperature.  Because  of  the  higher  coefficient 
of  thermal  expansion,  the  board  expands  more  than  the  chip  carrier,  which  results  in  the 
reversed  stresses  at  Condition  3.  The  board  and  the  component  remain  expanded  as  long 
as  the  maximum  temperatures  are  maintained.  Under  this  constant  solder  displacement, 
the  stresses  decrease  or  relax  as  discussed  in  Section  3.4.  The  stresses  decrease  very  rap¬ 
idly  to  Condition  4  in  Figure  5-21.  Upon  power-off.  the  temperature  in  the  case  de¬ 
creases  rapidly  as  the  remaining  heat  dissipates  into  the  board.  This  causes  a  contraction 
of  the  case  before  the  board  begins  to  contract,  and  creates  stresses  in  the  solder  (Condi¬ 
tion  5).  The  components  and  the  board  eventually  return  to  their  original  dimensions  as 
the  temperatures  decrease  to  the  original  levels  before  power  on.  This  results  in  reversed 
solder  stresses  (Condition  6),  which  will  also  relax  over  time  to  Condition  7. 
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Figure  5-21 .  Relaxation  of  Thermal  Stress  in  Solder 


Low  thermal  stress  cycles  can  be  ignored  to  simplify  the  above  thermal  stress  profile. 
This  truncation  procedure  is  used  to  filter  out  cycles  which  cause  negligible  fatigue  dam¬ 
age.  Examples  of  these  type  in  Figure  5-21  include  the  stress  cycles  caused  by  transient 
temperature  changes  during  power  on  and  power  off.  Although  the  component  case 
heats  up  first  during  power  on,  the  difference  in  temperature  between  the  case  and  the 
PCB  is  generally  low.  Initial  expansion  of  the  case  is  also  small  because  of  its  low  coeffi¬ 
cient  of  thermal  expansion.  Deleting  the  stress  cycle  between  Conditions  1  and  2  and  the 
cycle  between  Conditions  4  and  5  results  in  the  simpler  thermal  stress  profile  in  Figure 
5-22. 
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Stress 


It'  vibration  occurs  while  the  electronics  are  on,  the  majority  of  vibration  will  occur 
with  the  thermal  stresses  close  to  zero  (Figure  5-23).  Therefore,  the  vibration  and  ther¬ 
mal  stresses  are  essentially  independent  as  in  Figure  5-24,  and  their  effect  on  fatigue  life 
can  be  analyzed  with  the  model  discussed  in  Section  5.2.3. 1. 


Stress 


Figure  5-23.  Combined  Vibration  Stresses  and  Thermal  Stresses  in  Solder 
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Stress 


5.4  Fracture 

5.4.1  Brittle  Fracture 

Brittle  materials,  such  as  ceramics,  undergo  little  or  no  plastic  permanent  deformation 
under  large  stress  (Figure  5-25).  Fracture  (physical  separation  of  a  component  into  two 
or  more  parts)  occurs  when  the  stress  in  the  critical  area  reaches  the  ultimate  strength  of 
the  material.  To  check  if  failure  will  occur  under  the  operational  conditions,  the  stress  (o) 
computed  with  the  finite  element  analysis  should  be  compared  with  the  ultimate  strength 
(oUit).  The  material  will  fail  if: 


o  >  ou  K 


Eq.  31 


Brittle  Material 


Load  >  Oyit 

t 

□  Fractured 
Specimen 

U  0 

1 


Figure  5-25.  Brittle  Fracture 
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5.4.2  Ductile  Fracture 

As  discussed  in  Chapter  3,  ductile  materials  undergo  localized  plastic  deformation 
when  the  stress  in  the  critical  area  exceeds  the  yield  strength.  As  the  load  increases,  larg¬ 
er  volume  of  plastic  zone  increases  until  the  entire  cross  section  has  yielded.  Failure  oc¬ 
curs  when  the  ultimate  strength  is  exceeded  across  the  entire  section.  These  ductile 
failures  are  characterized  by  highly  deformed  broken  parts  (Figure  5-26).  Equation  3 1 
can  be  used  once  again  to  determine  if  ductile  failure  will  occur  under  the  operational 
environment.  However,  the  stresses  along  the  entire  cross  section  must  exceed  the  ulti¬ 
mate  strength.  If  only  a  localized  region  experiences  large  stress,  plastic  yielding  will  pre¬ 
vent  fracture. 

If  stresses  are  lower  than  the  ultimate  strength  of  the  material,  fracture  can  still  occur 
under  fatigue  or  creep.  Failure  from  repetitive  load  cycles  occurs  when  fatigue  cracks 
grow  to  a  critical  size  which  results  in  sudden  fracture  of  the  part.  A  sustained  load  act¬ 
ing  on  a  creep  sensitive  material  can  eventually  result  in  fracture  of  the  component.  Fa¬ 
tigue  and  creep  failure  theories  and  analysis  procedures  were  covered  in  Sections  5.2  and 
5.3,  respectively. 


Ductile  Material 

Load 

t 

p  Fractured 
II  Specimen 

J  0 

▼ 

GP03  C624-52-D  gms 

Figure  5-26.  Ductile  Fracture 

Leads  and  solder  joints  in  which  fracture  is  caused  by  fatigue  show  little  permanent 
plastic  deformation.  None  of  the  large  scale  deformation,  which  ductile  materials  experi¬ 
ence  during  a  single  large  load,  occurs  in  fatigue  failures.  Because  of  this,  fatigue  failures 
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are  sometimes  called  “brittle  failures”  although  the  material  may  be  ductile.  The  term 
brittle  failure  in  reference  to  fatigue  is  inaccurate,  since  localized  plastic  zones  do  occur 
around  fatigue  crack  tips  (Figure  5-27). 


Figure  5-27.  Plastic  Zone  Surrounding  a  Crack  Tip 


5.5  Buckling 

Buckling  can  occur  during  out  of  plane  deformation  of  a  printed  circuit  board.  The 
leads  connecting  components  to  the  board  act  as  columns  which  are  loaded  in  tension 
and  compression  (Figure  5-28).  To  determine  if  buckling  will  occur,  the  loads  acting  on 
the  leads  are  first  obtained  with  finite  element  analysis.  The  critical  buckling  load  can  be 
computed  by  using  the  following  formula  for  a  fixed-fixed  beam  (Reference  5-8): 


47T2  E  I 
L2 


Eq.  32 


where:  E  =  Young's  Modulus 

I  =  Moment  of  Inertia  -  Beam  Cross  Section 
L  =  Length  of  Beam 
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if  the  load  acting  on  the  lead  exceeds  the  critical  load,  the  lead  will  buckle.  If  some 
of  the  leads  buckle,  loads  get  redistributed  to  the  other  leads  on  the  component.  This  is 
a  very  complicated  problem  which  may  require  a  buckling  finite  element  analysis  of  the 
multiple  lead  system  to  obtain  actual  deformation.  Since  the  lead  can  return  to  its  origi¬ 
nal  shape  upon  unloading,  buckling  of  a  lead  may  not  necessarily  imply  failure  of  a  lead. 
However,  a  buckled  lead  may  impose  large  deformations  on  the  solder  joint,  and  may  re 
suit  in  solder  fatigue  damage.  A  finite  element  analysis  of  the  deflected  solder  joint  can 
then  be  used  to  determine  stresses  in  the  solder. 


mwmmwmw1 


Figure  5-28.  Out  of  Plane  Board  Deflections  Can  Cause  Buckling  of  Leads 
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Chapter  6 

RELIABILITY  APPLICATIONS 


6.0  Introduction 

This  chapter  describes  how  the  results  of  finite  element  analyses  are  used  to  predict 
fatigue  life  and  reliability  using  the  procedures  from  Chapter  5.  Several  examples  are  in¬ 
cluded  to  illustrate  the  procedures.  The  first  two  examples,  J-lead/soIder  joint  and  the 
leadless  chip  carrier  solder  joint,  include  the  finite  element  results  documented  in  Chap¬ 
ter  4.  The  third  example  describes  a  life  prediction  of  a  dual-in-line  package  lead/solder 
connection.  The  stresses  used  in  this  example  were  obtained  using  non-finite  element 
methods  to  illustrate  alternate  analysis  procedures.  The  fourth  example  covers  a  driver 
amplifier  using  MMIC  technology  (gallium  arsenide  chip  and  gold  metallization  with 
monolithic  circuitry  containing  six  transistor  cells).  This  component  is  part  of  a  Transmit/ 
Receive  Radar  Module.  The  life  of  the  gold  material  exposed  to  thermal  stresses  is  pre¬ 
dicted  in  this  example.  The  fifth  example  covers  procedures  for  checking  the  likelihood 
of  buckling  in  leads  subjected  to  compression  loads  during  vibration. 


6.1  J-Lead  Fatigue  Analysis 

The  following  is  a  list  of  the  basic  parameters  used  in  this  J-lead  analysis: 


Failure  Mechanism 

Failure  Theory 

Device 

Component 

Material 

Fnvironment 


-  Fatigue,  Creep  and  Stress  Relaxation 

-  Miner’s  Rule  on  Cumulative  Fatigue  Damage 

-  Flight  Control  Computer  Printed  Circuit  Board 

-  68  J-Lead  Surface  Mounted  Chip  Carrier 

-  Kovar  and  Solder 

-  On-Off  Temperature  Cycles  and  Aircraft  Vibration 


6.1.1  Thermal  Fatigue  Analysts 

A  J-lead/solder  joint  finite  element  analysis  was  completed  to  dotermine  the  stress  in 
critical  locations.  Different  finite  element  codes  and  boundary  conditions  were  used  to 
compare  various  ways  of  analyzing  the  same  problem.  These  finite  element  analyses. 
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described  in  Chapter  4,  simulated  a  temperature  increase  of  approximately  50  degrees  C. 
From  the  various  J-lead  finite  element  analyses,  one  was  chosen  to  illustrate  the  fatigue 
life  prediction  procedures.  This  analysis  consisted  of  a  PROBE  model  of  the  chip  carrier 
J-lead.  solder  joint,  and  circuit  board  under  an  imposed  temperature  change  (Figure  6-1 ). 
The  upper  lead  elbow  and  the  solder/lead  interface  were  selected  as  critical  locations 
where  fatigue  failure  could  occur  (Figure  6-2),  since  the  lead  elbow  has  the  highest  ther¬ 
mal  stress  in  the  lead  material  and  the  solder  experiences  large  thermal  stress  in  the  sol¬ 
der/lead  interface,  a  known  failure  location. 


Figure  6-1.  J-Lead/Solder  Joint  Geometry 


Critical  Location 
Solder 


Figure  6-2.  Critical  Areas  Analyzed  in  the  Thermal  Cycle  Reliability  Assessment 
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6. 1.1.1  Lead  Material  Fatigue 

The  upper  lead  bend  experiences  a  principal  stress,  Oi,  of  1400  psi  because  of  the  50°C 
temperature  increase.  The  maximum  temperature  is  maintained  long  enough  for  com¬ 
plete  stress  relaxation  and  creep  in  the  solder,  allowing  the  lead  to  return  to  its  original 
shape.  As  a  result,  the  stress  in  the  lead  material  decreases  to  zero. 


.As  the  temperature  decreases  to  its  initial  level,  the  lead  is  deformed  in  the  opposite 
direction  and  reversed  stresses  equal  to  -1400  psi  occur.  The  original  temperature  is 
maintained  long  enough  for  complete  stress  relaxation  and  creep  in  the  solder,  once 
again  resulting  in  zero  stress  in  the  lead  after  a  period  of  time  has  elapsed. 


This  stress  cycle  consists  of  fully  reversed  stress  (1400  to  -1400  psi)  and  a  stress  ratio 
of  -1.  Figure  6-  3  (Reference  6-1)  illustrates  the  lead  material  fatigue  data  used  to  pre¬ 
dict  the  life.  According  to  this  figure,  the  lead  material  can  sustain  more  than  109  cycles 
without  failure. 

io 5 
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102  103  104  105  10s  107  108  109  1010  1011 


Cycle  to  Failure  -  N 


GP03-0624-55-D  dot 


Figure  6-3.  Constant  Amplitude  S-N  Fatigue  Curve  for  Electrical  Lead  Wire, 
Type  D,  Type  K,  Reversed  Bending 

6. 1.1.2  Solder  Fatigue 


The  solder  experiences  a  maximum  shear  stress  of  1240  psi  near  the  solder/lead  inter¬ 
face.  The  maximum  temperature  is  maintained  long  enough  for  complete  stress  relax¬ 
ation  and  creep  in  the  solder  (Figure  6-4).  As  the  temperature  decreases  to  its  initial 
level,  a  reversed  stress  equal  to  -1240  psi  occurs  in  the  solder.  Stress  relaxation  and 
creep  causes  a  return  to  the  original  zero  stress  and  strain  state.  The  stress-strain 
oehavior  illustrated  in  Figure  6-4  was  described  by  Engelmaier  (Reference  6-2). 
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Figure  6-4.  Solder  Stress-Strain  Cycle  In  a  J-Lead  Connection 

The  fully  reversed  stress  cycles  with  amplitude  of  1240  psi  will  result  in  a  fatigue  life  of 
700,000  cycles,  as  read  from  the  S-N  curve  in  Figure  6-5.  The  stress  vs.  life  curve  in 
Figure  6-5  was  generated  from  strain  vs.  life  data  in  References  6-3  and  6-4. 


Figure  6-5.  S-N  Curve  for  63-37  Solder  at  Room  Temperature 


6.1.2  Vibration  Fatigue  Analysis 

A  J-lead/solder  joint  finite  element  analysis  was  completed  to  determine  the  stress  in 
critical  locations.  The  ABAQUS  finite  element  code  was  used  to  determine  deformations 
and  stresses  in  the  lead  and  solder  joint  material.  These  3-dimensional  finite  element 
analyses,  described  in  Chapter  4,  simulated  a  vibration  environment. 
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Based  on  the  finite  element  results,  two  critical  locations  were  selected  for  fatigue  life 
assessment.  Figure  6-6  illustrates  a  side  view  of  the  3-D  finite  element  model  and  the 
critical  locations.  The  first  critical  area  is  located  in  the  Kovar  material  where  the  lead 
first  meets  the  chip  carrier.  Large  vibration  stresses  in  the  Kovar  lead  material  occur  in  this 
region.  The  second  possible  failure  area  occurs  in  the  solder  near  the  solder/lead  interface. 


Figure  6-6.  Critical  Areas  Analyzed  In  the 
Vibration  Reliability  Assessment 


6. 1.2.1  Lead  Material  Fatigue 

The  principal  stress,  aj,  in  the  lead  critical  region  was  computed  with  ABAQUS  FEA 
to  be  14,390  psi  RMS.  Because  of  stress  relaxation,  thermal  stresses  are  approximately 
zero  by  the  time  vibration  starts.  As  a  result,  the  vibration  stresses  are  fully  reversed. 
Figure  6-7  (from  Reference  6.1)  presents  random  vibration  RMS  stress  vs.  cycles  to  fail¬ 
ure  data  for  lead  material.  The  fatigue  life  corresponding  to  a  14,390  psi  RMS  stress  is: 

N  =  4  x10s  cycles 
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Figure  6-7.  Random  Amplitude  S-N  Fatigue  Curve  for  Electrical  Lead  Wire, 
Type  D,  Type  K,  Reversed  Bending 


If  no  random  vibration  fatigue  data  is  available  for  the  material,  a  life  prediction  based 
on  constant  amplitude  fatigue  data  can  still  be  performed.  This  procedure,  described  in 
Section  5  3.2,  consists  of  assuming  a  Gaussian  distribution  of  random  vibration  stresses 
and  modeling  this  distribution  with  the  following  constant  amplitude  blocks  of  stress: 


68.3%  of  cycles  with  ampitude  =  o-rmc 
27.1%  of  cycles  with  amplitude  =  2  x  ctrms 
4.3%  of  cycles  with  amplitude  =  3  x  ctrMS 

This  results  in  the  following  peak  stresses  for  this  example: 


(CTmax)l  =  crRMS  =  14.390  psi 
(CTmax)2  =  2  X  ctrms  =  28,780  psi 
(CTmax)3  =  3  X  ctrms  =  43,170  psi 


The  number  of  cycles  to  failure  at  each  stress  level  can  be  read  from  Figure  6-3 
(Reference  6-1): 


At  (CTmax)i,  N[  =  1.2  x  10 13  cycles 
At  (CTmax)2,  N2  =  6.0  x  108  cycles 
At  (am ax)3,  N3  =  1.8  X  106  cycles 


Failure  is  detined  to  occur  when: 


ni  n2  n3 

—  +  — -  +  —  > 
Ni  N2  n3 


1.0 


Eq.  1 


Eq.2 
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where  rij  (i  =  1,  2,  3)  is  the  number  of  cycles  at  each  (omax),.  The  total  number  of  cycles 


to  failure  is: 


N  =  ni  4-  n2  +  n3 


Eq.  3 


From  the  assumed  Gaussian  distribution: 


ni  =  0.683  x  N 


n2  =  0.271  x  N 
n3  =  0.043  x  N 


Eq.  4 
Eq.  5 
Eq.  5 


Substituting  into  Equation  (2): 


xr  /  0.683  0.271  0.043 

I  N!  N2  N3 


=  1.0 


Eq.  7 


and  solving  for  N  gives  the  number  of  cycles  to  failure: 


xr  0.683  0.271  0.043 

N=  - + - + - 

Nt  N2  N3 


Eq.  8 


Substituting  the  values  of  Ni  from  Equation  1  into  Equation  8  results  in  the  following 
life  prediction  for  this  example: 


0.683  0.271  0.043 

1.2  x  1013  6.0  x  108  1.8  x  106 


Eq.  9 


=  4.1  x  107 


This  life  prediction  using  constant  amplitude  fatigue  data  was  shorter  than  the  predic¬ 
tion  of  4.0  x  108  cycles  based  on  random  vibration  data.  The  constant  amplitude  model 
was  designed  to  be  conservative  to  cover  inaccuracies  inherent  with  the  assumed  Gaus¬ 
sian  distribution.  This  was  achieved  by  assuming  that  any  cycle  with  amplitude  between  0 
and  orms,  has  an  amplitude  of  orms-  Any  between  Orms  and  2  x  Orms*  has  an  amplitude 
of  2  x  orms-  Similar  assumptions  apply  to  3  x  Orms-  resulting  in  a  conservative  model. 
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6.1. 2.2  Solder  Fatigue 

From  the  same  FEA  analysis,  the  maximum  shear  stress,  Tmax,  in  the  solder  critical  re¬ 
gion  was  computed  to  be  600  psi  RMS.  Once  again,  thermal  stress  relaxation  results  in 
fully  reversed  vibration  stresses.  Figure  6-8  (from  Reference  6-1)  presents  random  vi¬ 
bration  RMS  stress  vs.  cycles  to  failure  data  for  solder  material.  The  fatigue  life  corre¬ 
sponding  to  a  600  psi  RMS  stress  is: 

N  =  1.0  x  1010  cycles 


Figure  6-8.  Random  Amplitude  S-N  Fatigue  Curve  for  Soft  Solder 
(37%  Lead  -  63%  Tin),  Reversed  Shear 

The  constant  amplitude  procedure  is  also  used  in  this  example  to  compare  with  results 
using  the  random  stress  fatigue  data.  The  maximum  shear  stresses  for  the  Gaussian  in 
this  example  are: 

(^max)!  =  TRMS  ~  600  psi 
(Tnax)2  =  2  X  TrMS  =  1200  psi 
(7max)3  =  3  xrRMS  =  1800  psi 

The  number  of  cycles  to  failure  at  each  stress  level  can  be  read  from  Figure  6-9 
( Reference  6-1 ): 


At  (rmax)i,  Ni  =  6.0  x  1013  cycles 
At  (rmax)z  N2  =  3.5  x  1010  cycles 
At  (rmax)3,  N3  =  4.4  x  108 


Eq.  10 
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Figure  6  9.  Constant  Amplitude  S-N  Fatigue  Curve  for  Soft  Solder 
(37%  Lead  -  63%  Tin),  Reversed  Shear 


Substituting  into  Equation  8: 


=  /  0.683  0.271  +  0.043  \ 

~  \  6.0  x  10lJ  +  3.5  x  1010  +  4.4  x  108  ) 

=  9.5  x  109  cycles 


Once  again,  the  life  prediction  using  the  constant  amplitude  method  is  shorter  than  the 
prediction  based  on  random  stress  fatigue  data. 


6.2  LCC  Solder  Joint  Thermal  Fatigue  Analysis 

The  following  is  a  list  of  the  basic  parameters  used  in  this  leadless  chip  carrier  analysis. 


Failure  Mechanism 

Faiiuie  Theory 

Device 

Component 

Material 

Envi  ronment 


-  Low  Cycle  Fatigue  Cracking 

-  Strain  vs.  Cycles  to  Failure 

-  Flight  Control  Computer  Printed  Circuit  Board 

-  Leadless  Chip  Carrier  with  68  Solder  Joints 

-  Solder  63-37 

-  On-Off  Temperature  Cycles 


6.2.1  Thermal  Stress  Analysis 

A  leadless  chip  carrier  solder  joint  finite  element  analysis  was  completed  to  determine 
the  stress  in  critical  locations.  Different  finite  element  codes  and  boundary  conditions 
were  used  to  compare  various  ways  of  analyzing  the  same  problem.  These  finite  element 
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analyses,  described  in  Chapter  4,  simulated  a  temperature  increase  of  approximately  50 
degrees  C.  From  the  various  solder  joint  finite  element  analyses,  one  was  chosen  to  illus¬ 
trate  the  fatigue  life  prediction  procedures.  This  analysis  consisted  of  a  PROBE  model  of 
the  chip  carrier  and  solder  joint  under  an  imposed  temperature  change  (Figure  6-10). 

The  critical  location  in  the  solder/chip  carrier  interface  was  located  near  the  corner  of  the 
chip  carrier.  The  highest  shear  stresses  and  strains  were  found  in  this  area.  At  the  criti¬ 
cal  location,  the  elastic  finite  element  analysis  gave: 


Maximum  shear  stress  =  blOO  psi 
Maximum  shear  strain  =  0.00855  in/in. 


Figure  6-10.  2-D  Model  of  a  Leadless  Chip  Carrier  and  Solder  Joint 

The  shear  stress-strain  curve  in  Figure  6-11  was  constructed  with  data  from 
Reference  6-3  and  assuming  a  proportional  limit  of  2000  psi.  The  6100  psi  stress  obtained 
with  finite  element  analysis  indicates  that  plastic  deformation  occurs  in  the  solder  at  the 
critical  location. 
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Figure  6-11.  Shear  Stress-Strain  Curve  for  Solder  63-37 

The  procedures  outlined  in  Section  4.3.4  were  used  to  determine  the  actual  nonlinear 
stress  and  strain  state  from  the  linear  finite  element  analysis.  The  first  approximation  as¬ 
sumes  that  the  actual  shear  strain  is  the  same  as  the  calculated  strain.  The  actual  stress  is 
then  obtained  from  the  shear  stress-strain  curve  in  Figure  6-1 1.  This  procedure,  known 
as  the  "linear  rule",  gives  the  following  actual  conditions: 

Actual  maximum  shear  stress  =  3100  psi 
Actual  maximum  shear  strain  =  0.00855  in/in. 

T  he  second  procedure  used  the  following  "Neuber's  rule"  calculation: 

Elastic  FEA  stress  x  strain  =  6100  psi  x  0.00855  in/in  Eq.  12 

=  52.155  psi  in/in 

The  actual  stress  and  strain  are  obtained  by  finding  a  point  on  the  stress-strain  curve 
(Figure  6-1 1 )  that  gives  this  product.  Neuber's  rule  gives: 

Actual  maximum  sheai  stress  =  3600  psi 
Actual  maximum  shear  strain  =  0.0145  in/in. 

The  solder  experiences  these  peak  stresses  when  the  maximum  temperature  is  reached. 
I  hese  stresses  are  caused  when  the  PCB  expands  more  than  the  chip  carrier  whi<'h  induces 
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solder  joint  deflections.  The  temperature  and  deflection  is  maintained  until  stress  rci.iv 
ation  occurs  (Figure  6-12).  After  power-off.  the  temperature  returns  to  its  original  lew; 
and  the  stresses  are  reversed  as  the  circuit  board  contracts.  Stress  relaxation  returns  the 
solder  to  the  initial  state  of  stress. 


Figure  6-12.  Solder  Stress-Strain  Cycle 
in  a  Leadless  Solder 


Once  the  stress  and  strain  cycle  is  defined,  the  number  of  cycles  to  failure  can  be  pre¬ 
dicted  from  available  fatigue  data.  Two  alternate  procedures  are  used  in  this  example  te 
predict  the  life.  The  first  uses  the  calculated  stress  to  predict  life  from  a  plot  of  stress  \s. 
cycles.  The  second  method  is  based  on  effective  strain  to  predict  life. 

6.2.2  Fatigue  Life  Prediction  Based  on  Stress 

The  maximum  shear  stresses  estimated  from  the  finite  element  analysis  results  were 
used  to  predict  the  number  of  cycles  to  failure.  The  S-N  data  used  is  illustrated  in  Figure 
6-5.  The  following  fatigue  lives  were  obtained: 

"Linear  rule"  shear  stress  =  3100  psi 
Predicted  Life:  15.000  cycles 

“Neuber's  rule"  shear  stress  =  3600  psi 
Predicted  Life:  3000  cycles 
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The  range  of  fatigue  lives  obtained  with  this  method  is  conservative.  Although  the 
stresses  are  fully  reversed  in  Figure  6-12,  the  strain  varies  from  0  to  a  maximum  value. 
The  data  in  the  S-N  curve  in  Figure  6-5  corresponds  to  fully  reversed  stresses  and  strains, 
anu  therefore,  gives  shorter  lives  than  expected  for  th;s  example  problem. 

6.2.3  Fatigue  Life  Predictions  Based  on  Effective  Strain 

The  effective  strain  parameter  described  in  Chapter  5  was  used  to  address  the  effect  of 
non-reversed  strains  on  life.  The  strain  vs.  life  data  in  Figure  6-13  (References  6-3  and 
6-4)  was  converted  to  effective  strain  vs.  life.  Data  points  from  Figure  6-13  were 
selected  and  converted  by  using  the  following  relation: 


yeff  = 


Eq.  13 


where: 


G  =  Shear  Modulus  =  7.14  x  10s  psi 
Ay  =  Shear  strain  range  =  ymax  -  ymin 
rmax  =  Maximum  shear  stress  at  the  peak  of  the  cycle 


Strain 

Amplitude 

Ay/ 2 


7 


Cycle  to  Failure  -  N 

GP03-0621-63-D/dpt 


Figure  6-13.  Strain-Life  Curve  for  63-37  Solder  at  Room  Temperature 
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The  following  points  selected  from  Figure  6-13  illustrate  the  conversion: 


Ay/2 

N 

0.01 

104 

0.025 

103 

0.05 

102 

The  maximum  stress  of  the  cycle  can  be  read  from  the  stress-strain  curve  in  Figure 
6-11.  These  values  are  then  used  with  Equation  13  to  compute  the  following  effective 

strain  for  each  data  point: 

Ay/2 

I’max  (PSI) 

yeff 

0.010 

3250 

0.0067 

0.025 

4150 

0.0121 

0.050 

4900 

0.01S5 

Values  of  yeff  and  N  were  then  plotted  in  Figure  6-14 


Figure  6-14.  Effective  Strain  vs  Life  for  63-37  Solder 

The  next  step  in  the  analysis  is  to  determine  the  effective  strain  for  the  thermal  cycle 
used  in  this  example.  The  “linear  rule"  approximation  resulted  in  the  following  values  of 
maximum  stress  and  strain  range: 

Tnax  =  3100  iroi 
Ay  =  0.00855 
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Substituting  these  values  into  Equation  13  gives: 

■yeff  =  0.0043 


From  Figure  6-14: 


N  =  48,000  cycles 

Neuber’s  rule  resulted  in: 

Tmax  =  3600  psi 
Ay  =  0.0145 

Substituting  into  Equation  13  gives: 

yeff  =  0.00605 


and  from  Figure  6-14: 


N  =  15,000  cycles 

This  last  value  of  life  represents  a  reasonable  reliability  prediction  for  solder  under  the 
50  degree  C  thermal  cycles. 


6.3  Dual-Inline-Package  (DIP)  Lead/Solder  Fatigue  Analysis 

The  following  is  a  list  of  parameters  used  in  the  DIP  analysis. 


Failure  Mechanism  - 

Failure  Theory 

Device 

Component 

Material 

Environment 


Fatigue  Cracking 

Stress  vs.  Cycles  to  Failure 

Inertial  Sensor  Unit 

40  Pin  Dual-Inline-Package 

Kovar  and  Solder 

On-Off  Temperature  Cycles 
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This  example  illustrates  non-finite  element  procedures  for  analyzing  the  reliability  of 
lead  and  solder  material  connecting  a  dual-inline-package  to  a  printed  circuit  board 
(PCB).  The  example,  extracted  from  Reference  6-4,  is  included  as  an  alternate  proce¬ 
dure  for  rapidly  assessing  electronics  life  or  for  checking  finite  element  results. 

Consider  a  40  pin  DIP  soldered  to  a  PCB  that  must  operate  in  an  environment  where 
the  temperature  range  is  expected  to  vary  from  -55°C  to  105°C.  An  examination  of  the 
thermal  expansion  shows  that  the  PCB  will  expand  and  contract  more  than  the  ceramic  DrP 
component.  The  expansion  differences  will  force  the  electrical  lead  wires  to  bend  as  shown 
in  Figure  6-15.  A  deformation  equation  can  be  derived  to  describe  the  displacements  af¬ 
fecting  the  lead.  The  subscripts  P  and  C  represent  the  PCB  and  the  component  respectively: 


x: 


■Large  TCE 

Figure  6-15.  Thermal  Expansion  in  Plane 


GPO3  0183-37  D 


X  =  CTEp  Lp  AT  -  CTEc  Lc  AT 


Eq.  14 


where:  CTEp  =  coefficient  of  thermal  expansion  of  an  epoxy  fiberglass  circuit  board 
=  15  x  10'6  in/in/C 
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CTEc  =  coefficient  of  thermal  expansion  of  ceramic  DIP  body 
=  6  x  10~6  in/in/C 

Lc  =  length  from  center  to  end  of  DIP 
=  1.0  inch 

Lp  =  effective  length  of  PCB,  same  as  DIP 
=  1.0  inch 

AT  =  temperature  amplitude 
=  [105  -(-55)]/2  =  80°C 

X  =  relative  displacement 

Note  that  half  the  maximum  to  minimum  temperature  range  is  used  to  find  the  displace¬ 
ment  difference  between  the  PCB  and  DIP.  This  assumes  that  the  expansion  from  the 
center  of  the  component  to  one  end  of  the  component  is  the  same  as  the  contraction. 
This  also  assumes  that  the  cycle  duration  does  not  permit  stress  relaxation  and  creep  in 
the  solder. 

Substituting  into  Equation  14  gives  the  expansion  difference  between  the  PCB  and 
DIP: 


X  =  (15  -  6)  x  10-6  (1.0)(80) 

=  0.00072  inches  Eq.  15 

The  horizontal  force  acting  on  the  electrical  leads  can  be  determined  from  the  spring 
rate  of  the  lead  and  the  displacement  calculated  below.  Therefore,  the  lead  spring  rate 
must  be  obtained  in  order  to  obtain  the  force  in  the  lead.  The  standard  DIP  lead  has  a 
small  horizontal  leg  that  extends  outward  from  the  DIP  body  before  it  makes  a  90  degree 
bend  down,  as  shown  in  Figure  6-16.  Therefore,  when  the  DIP  lead  is  forced  to  bend 
due  to  the  thermal  expansion  difference,  the  vertical  part  of  the  lead  will  bend  and  the 
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Figure  6-16.  Typical  DIP  Package 

horizontal  part  of  the  lead  will  twist,  as  shown  in  Figure  6-17.  The  problem  can  be  sim¬ 
plified  by  considering  only  the  end  lead  on  the  DIP  for  the  bending  stiffness.  Since  the 
stiffness  of  the  DIP  and  the  PCB  along  the  X  axis  is  much  greater  than  the  bending  stif¬ 
fness  of  all  the  leads,  there  will  be  very  little  change  in  the  relative  displacement  when 
only  the  end  DIP  leads  are  used  for  the  analysis. 

Superposition  can  be  used  to  find  the  total  displacement  due  to  bending  and  twisting. 
The  lead  can  be  considered  as  a  beam  fixed  at  both  ends,  with  a  lateral  displacement  as 
follows: 

p  ^  3 

Xi  =  — —  =  Lead  bending  displacement  Eq.  16 

12  E  I[ 


where: 


P  =  load 

E  =  Young’s  Modulus 
Ii  =  Moment  of  Inertia  of  lead 
Li  =  Lead  length,  vertical 
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TOTAL  DEFLECTION  BENDING  TWIST 


Figure  6-17.  DIP  Lead  Deformation 


Considering  the  torsion  in  the  standard  DIP  lead  next,  the  horizontal  part  of  the  lead 
will  twist  through  the  angle  theta  (9)  resulting  in  a  displacement  for  this  segment  of: 

X2  =  Li  0  Displacement  due  to  wire  rotation  Eq.  17 

The  angle  of  twist,  theta,  will  be  related  to  the  torque  and  the  torsional  stiffness  of  the 
lead  as  follows: 

e  -  7TT  Et>- 18 

where:  T  =  Lead  Torque 

L2  =  Lead  length,  horizontal 
J  =  Torsion  form  factor 
G  =  Shear  Modulus 

This  torque  is  equal  to: 

T  =  P  Li  Eq.  19 
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Substituting  Equations  18  and  19  into  Equation  17  gives: 


X2  - 


p  U  2  L2 
G  J 


Eq.  20 


The  total  lead  displacement  will  be  the  sum  of  the  two  displacements: 

X  =  Xi  +  X2  Eq.  21 

Substituting  Equations  16  and  20  into  Equation  21: 


X  =  P  U 2 


12  E  Ii 


Eq.  22 


The  load  P  on  the  lead  is: 


Eq.  23 


where:  Li  =  Length  of  vertical  wire  leg,  which  includes  part  of  the  thinner  wire  that 

extends  into  solder  joint  in  the  plated  through  hole,  one  wire  width  of 
0.018  inch. 

=  Lu  +  Llb  =  0.156  +  0.018  =  0.174  inch 
L2  =  Length  of  horizontal  wire  leg  =  0.038  inch 
d  =  constant  wire  thickness  =  0.010  inch 


h  =  width  at  wide  part  of  wire  =  0.050  inch 
h  ib  =  width  at  narrow  part  of  wire  =  0.018  inch 
I ia  =  d  h3/12  =  (0.010)  (0.050)3/12  =  1.042  x  10'7  in.4 
lib  =  (0.010)  (0.018)3/12  =  4.86  x  10'9  in.4 
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E  =  modulus  of  elasticity  for  Kovar  =  20  x  106  psi 

J  =  wire  torsion  form  factor  =  h  d3/3 
=  (0.050)  (0.0 10)3/3  =  1.67  x  10'8  in.4 

G  =  Shear  modulus  for  Kovar  =  8.27  x  106  psi 

A  weighted  average  lead  moment  of  inertia  was  used  to  compensate  for  the  sharp  re¬ 
duction  in  the  wire  cross  section  at  the  narrow  segment  of  the  DIP  lead. 


Iav  - 


Lia  I  la  +  Lib  fib 

Lla  +  Lih 


Average  moment  of  inertia 


Eq.  24 


Substituting  the  various  values  results  in: 


Iav 
Eq.  25 


(0.156)  (1,042  x  IQ"7  )  +  (0.018)  (4,89  x  IQ'9  ) 
0.156  +  0.018 


=  9.39  x  10'8  in4 


Substituting  the  above  values  into  Equation  23  results  in  the  load  on  the  lead: 

0.00072 


P  = 


Eq.  26 


(0.174)2 


0.174 


+ 


0.038 


(12)  (20  x  106)  (9.39  x  10'8)  (8.27  x  106)  (1.67  x  10'8 


=  0.084  lb 


The  bending  moment  acting  on  the  lead  at  the  solder  joint  can  be  obtained  with  the 
following  relation: 

M  _  PL,  _  (O084HO174)  Eq  27 

2  2 

=  0.0073  in -lb 

The  bending  stress  at  the  narrow  portion  of  the  lead  can  be  determined  from  standard 
bending  stress  equation: 
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M  C  (0.0073)  (0.009) 

I  4.89  x  10-9 

=  13,450  psi 


Eq.28 


Since  the  endurance  limit  for  Kovar  is  about  40,000  psi,  the  design  is  safe  for  the  lead 
material. 


The  shear  tearout  stress  in  the  solder  joint  is  considered  next.  The  stress  state  is 
caused  by  forces  which  pull  the  lead  out  of  the  plated  through  hole  as  the  board  expands. 
The  expression  for  this  relation  is: 

M 

Sst  =  r -  Solder  shear  tearout  stress  Eq.  29 

hp  As 

where:  M  =  Bending  Moment  =  0.0073  in-lb 

hp  =  PCB  thickness  =  0.080  in. 

A§  =  Area  of  solder  section,  using  an  average  radii  of  the 
solder  joint  equal  to  0.014  in. 

=  v  (0.014)2  =  0.000616  in2 


Substituting  results  in: 


Sst 


0.0073 

(0.080)  (0.000616) 
148  psi 


Eq.30 


An  examination  of  the  solder  fatigue  curve  in  Figure  6-5  shows  that  the  fatigue  life  for 
solder  will  exceed  107  cycles. 

6.4  Case  Study  of  MMIC  GaAs  Component  for  Thermal  Reliability 

The  reliability  of  a  driver  amplifier  using  Microwave/Millimeter- Wave  Monolithic  In¬ 
tegrated  Circuits  (MMIC)  technology  is  examined  in  this  example.  This  component  is  part 
of  a  Transmit/Receive  Radar  Module  as  shown  in  Figure  6-18.  The  highest  heat  producing 
component  in  the  module  is  the  driver  amplifier  integrated  circuit  (IC)  (Figure  6-19)  which 
p.oduces  a  peak  value  of  9  watts  of  heat  and  an  average  value  of  2.31  watts.  This  critical 
component  uses  MMIC  technology  consisting  of  a  gallium  arsenide  (GaAs)  chip  and  gold 
metallization  with  monolithic  circuitry  containing  six  field  effect  transistor  (FET)  cells.  The 
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Figure  6-18, 


Figure  6-19.  Driver  Amplifier 
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module  is  attached  to  a  controlled  heat  sink,  and  the  driver  amplifier  duty  cycle  consists  of  1 
millisecond  (msec)  on/3  msec  off  pulses.  The  pulsing  of  the  amplifier  causes  temperatures  in 
the  component  to  vary  with  time.  These  temperatures  vary  above  and  below  the  average 
values  that  would  be  obtained  from  steady  state,  time  averaged,  heat  dissipation  thermal 
simulation.  Thus,  the  stresses  vary  with  time  and  have  a  non-zero  mean  stress.  Hence, 
fatigue  data  (S-N)  must  be  modified  to  account  for  non-zero  mean  stresses  (or  R>  1). 

6.4.1  Finite  Element  Analysis  Results 

Thermal  stress  FEA  was  performed  on  the  component  (see  Reference  6-5  for  further 
detail  on  the  thermal  finite  element  modeling  and  analysis).  Due  to  the  high  concentra¬ 
tion  of  heat  being  generated  in  a  small  area,  calculated  stresses  in  some  regions  of  the 
device  are  beyond  the  material's  elastic  limit.  Thus,  the  linear  FEA  results  were  modi¬ 
fied  to  estimate  true  plastic  stresses  and  strains  using  Neuber's  Rule. 

The  six  FET  cells  are  represented  by  the  FEM  showm  in  Figure  6-20.  Only  three  cells 
on  1/2  of  the  chip  needed  to  be  modeled  because  of  thermal  symmetry  along  the  full  chip 
center  line.  Also  shown  are  the  adiabatic  lines  of  symmetry  that  define  the  boundaries  of 
the  follow-on  1/4  FET  cell  model. 

The  temperature  distribution  shown  in  Figure  6-21  is  from  the  thermal  transient  analy¬ 
sis  of  the  chip  and  surrounding  surface  metal  at  one  time  step.  At  this  time  step.  Figure 
6-21  represents  temperatures  at  average  conditions  while  Figure  6-22  shows  the  stresses 
in  the  gold  metallization  at  the  same  instant.  Figure  6-23  shows  that  temperatures  vary 
above  and  below  this  condition.  Likewise,  the  stresses  wall  also  vary  above  and  below 
these  average  values  shown  in  Figure  6-22.  For  reference,  the  stresses  are  in  units  of 
newtons  per  square  micrometer  (N/pm2).  To  acquire  maximum  and  minimum  thermal 
stresses,  all  nodal  temperatures  at  the  maximum  and  minimum  temperature  times  (either 
1  or  5  milliseconds  and  4  milliseconds  respectively,  on  Figure  6-23)  are  read  into  the  ther¬ 
mal  stress  analysis  files.  The  zero  stress  state  is  assumed  to  be  at  room  temperature.  20°C. 
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6.4.2  Gold  Metalization  Circuitry  Reliability  Prediction 

The  following  is  a  list  of  the  basic  parameters  and  assumptions  used  in  this  analysis. 


Failure  Mechanism 
Failure  Theory 
Device 
Component 


Material 

Environment 


-  Fatigue  Cracking 

-  Coffin-Manson  Plastic  Strain  vs.  Life  Analysis 

-  Transmit/Receive  Radar  Module 

-  Driver  Amplifier  using  MMIC  technology 
(gallium  arsenide  chip  and  gold  metallization  with 
monolithic  circuitry  containing  6  transistor  cells) 

-  Commercial  Gold 

-  Module  attached  to  a  controlled  heat  sink. 

Driver  amplifier  pulsating  with  a  1  millisecond 
on/3  millisecond  off  pulse 


6.4.2.1  Linear  Stresses  and  Neuber  Rule  Approximations 

Stress  and  strain  within  certain  areas  of  the  gold  metallization  are  at  stresses  that  ex¬ 
ceed  the  elastic  limit.  It  is  at  these  locations  where  the  following  stresses  were  acquired 
using  the  elastic  FEA: 

Maximum  Tensile  Stress  -  11,387  psi 
Average  Tensile  Stress  -  8,540  psi 
Minimum  Tensile  Stress  -  7,246  psi 

Note  that  for  commercial  gold,  the  material’s  yield  point  is  approximately  1200  psi. 
and  therefore,  the  gold  metallization  is  undergoing  plastic  deformations.  In  order  to 
modify  the  linear  stress  and  strain  to  account  for  plastic  deformation,  Neuber’s  relation¬ 
ship  is  used.  Two  items  are  necessary  for  this  calculation:  the  material's  stress-strain 
curve  in  Figure  6-24  and  the  equation: 
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Stress 


(a  X  factual  —  ip  X  £)fea 


Eq,  31 


where:  o  =  stress  (psi) 

e  =  strain  (in/in) 

The  linear  stress  and  strain  are: 


ctfea  =  11,387  psi 

epEA  =  o/E  =  11,387  /  10.8  x  106  =  0.001054 


Eq.  32 


where:  E  =  modulus  of  elasticity 

Finally: 


(o  x  €)fea  =  (11,387)  (0.001054)  =  12.0  Eq.  33 

The  point  on  the  stress-strain  curve  in  Figure  6-24  which  gives  this  product  is: 

^actual  =  3,703  psi 
^actual  =  0.00324 

The  above  values  are  the  stress  and  strain  at  maximum  temperature: 

crmax  =  3,703  psi 
€max  =  0,00324 
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A  similar  procedure  can  be  followed  to  find  the  stress  and  strain  at  minimum  tempera¬ 
ture.  First  the  change  in  the  linear  stress  due  to  the  drop  in  temperature  is  computed  as 
follows: 


Aofea  =  11,387  -  7,246  =  4,141  psi  Eq.  34 

The  change  in  strain  is  given  by: 

AepEA  =  4,141  /  10.8  x  106  =  0.0003834  Eq.  35 

The  actual  change  in  stress  and  strain  can  be  determined  by  using  the  hysteresis  curve 
in  Figure  6-25  and  Neuber’s  relation  one  more  time.  The  hysteresis  curve  was  generated 
by  multiplying  o  and  e  values  of  selected  data  points  on  the  stress-strain  curve  (Figure  6-24) 
by  a  factor  of  two: 


Ao  =  2xo 
Ae  =  2  x  € 


Ao 


Figure  6-25.  Gold  Hysteresis  Curve 


These  points  were  then  plotted  as  Ao  vs.  Ae.  The  hysteresis  curve  is  basically  a  re¬ 
scaled  stress-strain  curve  which  represents  the  assumed  behavior  of  the  material  during 
unloading.  Neuber’s  rule  in  this  case  becomes: 


(Ao  x  Ae)acnia,  =  (Ao  x  Ae)reA 


Eq.  36 


where:  (Ao  x  A€)feA  =  (4,141)  (0.0003834)  =  1.588 


Eq.  37 
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The  actual  Ao  and  Ae  values  are  then  read  from  Figure  6-25  at  a  point  on  the 
hysteresis  curve  which  gives  the  above  product.  The  actual  values  are: 

A  (Tactual  =  2700  psi 
Aeactuai  =  0.00059 

The  stress  and  strain  at  the  minimum  temperature  are: 

CTmin  =  3,703  -  2,700  =  1003  psi  Eq.  38 

€min  =  0.00324  -  0.00059  =  0.00265 

In  summary,  the  stress  and  strain  values  at  the  maximum  an  minimum  temperatures 
are: 


FEA  Stress  Actual  Stress 


Actual  Strain 


Max.  -  1 1 ,  387  psi  3,703  psi  0.00324  in/in 

Min.  -  7,246  psi  1 ,003  psi  0.00265  in/in 


6.4.2.2  Thermal  Cycles  to  Failure 

The  Coffin-Manson  formula  that  relates  strain  to  number  of  cycles  to  failure  will  be 
used.  The  formula  is  modified  such  that  only  plastic  strain  will  be  considered.  The  strain 
values  acquired  from  the  Neuber  analysis  include  both  the  elastic  and  plastic  portions. 
However,  the  strain  value  will  be  assumed  to  be  all  plastic.  Figure  6-24  shows  this  sim¬ 
plification  to  be  justified.  The  Coffin-Manson  formula  (References  6-6  and  6-7)  is: 


N 


Eq.  39 


where:  D 

c 

Ae/2 


material  ductility  (in/in) 

ductility  exponent 

strain  amplitude 

one  half  of  strain  range  (in/in) 
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This  formula  applies  to  fully  reversed  cyclic  loading  with  a  zero  mean  stress.  Our 
loading  is  not  fully  reversed  (Figure  6-26).  In  order  to  account  for  the  change  in  fatigue 
life  caused  by  a  non-zero  mean  stress  cyclic  loading,  effective  strains,  eeff,  are  computed. 
The  number  of  cycles  to  failure  are  then  calculated  using  eeff  vs.  number  of  cycles  to  fail¬ 
ure  data.  The  effective  strain  is  given  by: 


Eq.  40 


where:  eeff  =  the  effective  strain  (in/in) 

E  =  modulus  of  elasticity  =  10.8  x  106  psi 
Ae  =  strain  range  (in/in) 
omax  =  the  maximum  stress  in  the  cycle 

Stress 


o 


GP03-0624-67-D/qg 


Figure  6-26.  Stress  History 


The  temperature  excursions  experienced  by  the  component  results  in  the  following 
stress  and  strains  in  the  gold: 


anax  =  3,703  iroi,  e^ax  =  0.00324  iv/iv 
<Vv  =  1,003  7TCT1,  e^lv  =  0.00265  iv/iv 

where:  De  =  0.00324  -  0.00265  =  0.00059  in/in  Eq.  41 

and:  De/2  =  0.00059/2  =  0.000295  in/in  Eq.  42 
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From  the  above  equation: 


feff  = 


3703 

10.8E6 


0.000295  =  .00032 


Eq.  43 


Data  relating  eeff  to  the  number  of  cycles  to  failure  can  be  determined  with  the  Coffin- 
Manson  equation.  Using  a  value  of  ductility  =  0.3  and  a  ductility  exponent  =  0.5  for 
commercial  gold,  a  plot  of  ecff  vs.  Life  (N)  is  determined  as  follows: 

1)  Using  various  values  of  maximum  stress,  determine  plastic  strain  amplitude  (Ae/2) 
from  the  stress-strain  curve  (Figure  6-24). 

2)  Calculate  the  number  of  cycles  to  failure  (N)  using  Equation  39. 

3)  Calculate  the  effective  strain  (€eff)  using  Equation  40  where  o  =  maximum  stress 
for  each  Ae/2  and  E  =  10.8xl06  psi. 

Sample  data  points  for  the  eeff  vs.  N  curve  are: 


Maximum  Stress 

Strain  Amplitude 

N 

eeff 

v  (psi) 

De/2  (in/in) 

(cycles) 

(in/in) 

3,703 

0.00324 

4,287 

0.0011 

3,162 

0.00212 

10,012 

0.0008 

2,898 

0.00168 

15,944 

0.0007 

2,500 

0.00120 

31,250 

0.0005 

2,200 

0.00080 

70,312 

0.0004 

2,000 

0.00066 

103,660 

0.00035 

1,900 

0.00058 

132,270 

0.00032 

For  this  loading,  eeff  =  0.00032,  and  N  =  132,270  cycles  to  failure  (Figure  6-27). 


Effective 

Strain 


Eeff 


0.00032 


Figure  6-27.  Fatigue  Life  From  an  Effective  Strain  vs  Life  Curve 


6-38 


Computer  Aided  Assessment  of  Reliability 
_ Using  Finite  Element  Methods 


6.4.3  Gallium  Arsenide  Reliability  Prediction 


Failure  Mechanism  - 
Failure  Theory 
Device 
Component 


Material 

Environment 


Brittle  Fracture 

Maximum  Stress  Failure  Criteria 
Transmit/Receive  Radar  Module 
Driver  Amplifier  using  MMIC  technology  (gallium 
arsenide  chip  and  gold  metallization  with 
monolithic  circuitry  containing  6  transistor  cells) 
Gallium  Arsenide 

Module  attached  to  a  controlled  heat  sink. 

Driver  amplifier  pulsating  with  a  1  millisecond 
on/3  millisecond  off  pulse 


Gallium  arsenide  is  a  brittle  material  with  an  ultimate  strength  of  6100  psi.  The  fail¬ 
ure  mechanism  is  brittle  fracture,  and  the  component  will  break  if  the  stresses  reach  or 
exceed  the  strength  of  the  material.  The  finite  element  analysis  indicates  that  the  thermal 
stress  in  the  gallium  arsenide  can  reach  12,035  psi.  This  results  in: 

°max  ->  CTu]t  Eq.  44 

where:  omax  =  12,035  psi 
ouit  =  6,100  psi 

Therefore,  the  material  is  expected  to  fail  at  the  operational  temperatures.  These  tem¬ 
peratures  need  to  be  lowered  to  bring  the  thermal  stress  below  the  ultimate  strength  of 
gallium  arsenide. 


6.5  Lead  Buckling  Analysis 


Failure  Mechanism 

Failure  Theory 

Device 

Component 

Material 

Environment 


-  Buckling 

-  Euler’s  Buckling  Formula 

-  Flight  Control  Computer  Printed  Circuit  Board 

-  68  J-Lead  Surface  Mounted  Chip  Carrier 

-  Kovar 

-  Aircraft  Vibration 
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A  finite  element  analysis  simulating  a  vibration  environment  is  described  in  Chapter  4. 
This  analysis  indicates  that  some  of  the  leads  connecting  the  surface  mount  component 
experience  compression  while  the  board  vibrates.  The  load  (Frms)  acting  on  the  critical 
lead  was  computed  with  the  finite  element  analysis  to  be  0.324  lb.  During  random  vibra¬ 
tion,  peak  loads  can  reach  levels  four  times  as  high  as  Frms: 

Maximum  Applied  Load,  Fapplied  =  4  x  Frms  Eq.  45 

=  4  x  0.324  lb 
=  1.296  lb 


The  critical  buckling  load  is  given  by: 


F cr  = 


4?rE  I 
L2 


c.q.  46 


where:  E  =  modulus  of  elasticity  =  20  x  106  psi 

L  =  length  of  lead  =  0.1  inch 
I  =  moment  of  inertia  of  lead  cross  section 
=  1/12  (width)  (thickness)3 
=  1/12  (0.012)  (0.008)3 
=  5.12  x  10-10 

Substituting  into  Equation  (46): 


^  4 jf_  (20  x  IQ6)  (5.12  x  10~10) 

(0.1)2 

=  40.4  lbs 


Eq.  47 


Because  the  applied  load  (1.296  lb)  is  much  smaller  than  the  critical  load  (40.4  lbs), 
the  lead  will  not  buckle. 
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